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Abstract 

As  part  of  the  Shallow  Water  Acoustics  in  a  Random  Medium  (SWARM)  experiment 
[1],  a  sixteen  element  WHOI  vertical  line  array  (WVLA)  was  moored  in  70  meters 
of  water  off  the  New  Jersey  coast.  This  array  was  sampled  at  1395  Hz  or  higher  for 
the  seven  days  it  was  deployed.  Tomography  sources  with  carrier  frequencies  of  224 
and  400  Hz  were  moored  about  32  km  shoreward,  such  that  the  acoustic  path  was 
anti-parallel  to  the  primary  propagation  direction  for  shelf  generated  internal  wave 
solitons.  Two  models  for  the  propagation  of  normal  modes  through  a  2-D  waveg- 
uide with  solitary  internal  wave  (soliton)  scattering  included  are  developed  to  help 
in  understanding  the  very  complicated  mode  arrivals  seen  at  the  WVLA.  The  sim- 
plest model  uses  the  Preisig  and  Duda  [2]  sharp  interface  approximation  for  solitons, 
allowing  for  rapid  analysis  of  the  effects  of  various  numbers  of  solitons  on  mode  ar- 
rival statistics.  The  second  model,  using  SWARM  thermistor  string  data  to  simulate 
the  actual  SWARM  waveguides,  is  more  realistic,  but  much  slower.  The  analysis  of 
the  actual  WVLA  data  yields  spread,  bias,  wander,  and  intensity  fluctuation  signals 
that  are  modulated  at  tidal  frequencies.  The  signals  are  consistent  with  predicted 
relationships  to  the  internal  wave  distributions  in  the  waveguides. 
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1-1  A  simple  two-layer  model  of  an  internal  soliton,  with  the  upper  layer, 
hi  =  15  m,  lower  layer,  h2  =  55  m,  and  reduced  gravity,  gAp/p  =  .01 
m/s2.  The  20  m  depression  propagates  at  a  phase  speed,  V  =  .8  m/s.        28 

1-2  Time  and  space  interpolation  of  four  thermistor  records.  The  thermis- 
tors were  attached  to  the  WVLA  at  water  depths  of  12.5,  20.5,  30.5, 
and  40.5  meters  and  sampled  at  30  second  intervals 29 

1-3  Space  transformation  of  the  time  series  for  the  prominent  soliton  in  the 
lower  frame  of  Figure  1-2.  All  six  available  thermistors  are  used,  so  the 
depth  spans  12.5  to  60.5  meters.  The  dimensions  shown  are  properly 
scaled  for  soliton  with  a  combined  phase  and  advection  speed  of  .8  m/s.  29 

2-1  Source  positions  are  noted  by  "*".  The  400  Hz  source  is  the  shallower 
of  the  two.  WVLA  hydrophone  positions  are  noted  by  "o"  and  the 
thermistors  on  the  array  are  marked  as  "+" 34 

2-2  Circular  convolution  of  an  ideal  delayed  400  Hz  transmitted  sequence 
with  an  ideal  non-delayed  transmitted  sequence.  Both  sequences  have 
been  zero-padded  from  length  7129  to  length  8192.  The  lower  frame  is 
a  window  of  the  full  output  above 38 

2-3     Typical  400  Hz  arrivals  at  the  WVLA  hydrophones.    The  sets  of  16 

arrivals  are  separated  by  6  minutes 39 
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3-1  Mode  rejection  levels  from  mode-filter  outputs,  Ai(*),  A2(x),  ^3(0), 
and  A4(+),  as  a  function  of  array  tilt.  The  four  plots  represent  tilted 
inputs  of  pure  mode  1,  mode  2,  mode  3,  and  mode  4  arrivals.  These 
plots  use  mode  filters  and  mode  arrivals  generated  from  the  same  typ- 
ical SWARM  SSP 46 

3-2  Mode  rejection  levels  from  mode-filter  outputs,  Ai(*),  A2(x),  A3(o), 
and  A4(+),  as  a  function  of  transmission  number.  The  four  plots  are 
based  on  a  series  of  untilted  inputs  of  pure  mode  1,  mode  2,  mode 
3,  and  mode  4  arrivals,  where  the  mode  filters  and  mode  arrivals  are 
generated  from  the  same  series  of  SSPs 48 

3-3  Mode  rejection  levels  from  mode-filter  outputs,  A}(*),  A2(x),  A3(o), 
and  Ai(-l-),  as  a  function  of  transmission  number.  The  four  plots  are 
based  on  1  degree  tilted  inputs  of  pure  mode  1,  mode  2,  mode  3,  and 
mode  4  arrivals,  using  mode  filters  and  mode  arrivals  generated  from 
the  SSPs  separated  by  30  seconds. 50 

3-4  Mode  Coefficient  outputs,  A^solid),  A2(dashed),  A3(dot-dash),  and 
A4 (dotted),  as  a  function  of  time.  The  ten  plots  depict  the  arrival 
structure  for  the  first  of  22  captured  sequences  out  of  ten  sequential 
transmissions 52 

3-5  Mode  coefficient  outputs,  A\ (solid),  A2 (dashed),  A3(dot-dash),  and 
A4 (dotted),  as  a  function  of  time.  The  ten  plots  depict  the  arrival 
structure  for  every  other  sequence  from  the  first  to  the  19th  of  22 
captured  sequences  from  the  time  zero  transmission  in  Figure  3-4.    .  .       54 

3-6     Average  correlation  coefficients  for  outputs,  A\ (solid),  A2(dashed),  ^(dot- 
dash),  and  A4(dotted),  as  a  function  of  elapsed  time.  The  average  is 
taken  over  the  ten  transmissions  of  Figure  3-4 55 
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4-1  An  illustration  of  the  SI  A  approximation  for  solitons  (left),  and  the 
idealized  SSPs  (right)  used  to  model  the  interior  and  exterior  regions 
of  the  SIA  soliton.  The  dotted  curve  represents  the  exterior  or  back- 
ground region,  and  the  solid  curve  represents  the  interior 59 

4-2  Propagation  of  four  modes  through  32  km  of  an  idealized  wave  guide 
with  one  200  meter  soliton  at  the  20  km  point.  The  solid  curves  are 
the  absolute  value  of  coherent  sums  of  the  J2  =  16  dotted  individually 
scattered  arrivals.  Starting  amplitudes  for  mode  1  through  4  were  1.0, 
0.8,  0.6,  and  0.4  respectively 60 

4-3  Propagation  of  four  modes  through  32  km  of  an  idealized  wave  guide 
with  200  meter  solitons  at  the  10  and  20  km  points.  The  solid  curves  are 
the  coherent  sums  of  the  N2  x  N2  =  256  dotted  individually  scattered 
arrivals.  Initial  mode  1  through  4  amplitudes  were  1.0,  0.8,  0.6,  and 
0.4  respectively 62 

4-4  Simulated  mode  coefficient  outputs,  A\ (solid),  A2(dashed),  ^(dot- 
dash),  and  Aj(dotted),  for  propagation  through  seven  idealized  soli- 
tons. The  solitons  propagated  toward  the  receiver  in  lock  step  at  0.8 
m/s.  The  ten  plots  depict  the  arrival  structures  for  sequential  realiza- 
tions at  six  minute  intervals 64 

4-5     Time-series  of  temperature  at  a  depth  of  12.5  meters  on  the  WVLA.  .       65 

4-6  Field  measured  SSPs  for  use  in  the  pulse-propagation  model.  The 
dashed  line  is  representative  of  a  standard  SSP  with  no  soliton.  The 
solid  line  represents  the  soliton  SSP 66 

4-7  Mode  shapes  for  field  measured  (SWARM)  SSPs.  The  dashed  lines  are 
the  355  Hz  (thin  lines)  and  445  Hz  (thick  lines)  mode  shapes  corre- 
sponding to  the  standard  SSP  with  no  soliton.  The  solid  lines  are  the 
same  modes  for  the  soliton  SSP.  The  dotted  line  at  68  m  indicates  the 
sediment  interface 68 
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4-8  Phase  and  group  speeds  for  the  standard  field  measured  SSP.  Mode  1 
has  the  lowest  phase  speed  and  the  highest  group  speed.  The  remaining 
phase  speeds  increase  monotonically  with  mode  number,  and  the  group 
speeds  decrease  monotonically  with  mode  number 69 

4-9  Effects  of  mode  and  frequency  dispersion  on  single  mode  transmissions 
over  32  km.  The  transmissions  include  no  soliton  scattering  or  intrinsic 
mode  attenuation.  The  attenuation  in  the  higher  modes  is  due  solely 
to  frequency  dispersion 71 

4-10  Effects  of  mode  and  frequency  dispersion  combined  with  intrinsic  mode 
attenuation  on  single  mode  transmissions  over  10  km.  The  transmis- 
sions include  no  soliton  scattering 72 

4-11  Mode  coupling  effects  of  the  bottom  topography  shown  in  Figure  2-1 
when  the  SSP  is  a  modified  version  the  "non-soliton"  SSP  of  Figure 
4-6.  The  "No  Ducting"  modification  consists  of  imposing  the  50.5  m 
temperature  from  50.5  m  to  the  bottom  depth,  and  the  "Ducting" 
modification  consists  of  imposing  the  60.5  m  temperature  from  50.5 
m  to  the  bottom  depth.  The  thinner  lines  in  this  figure  represent 
predicted  arrival  shapes  for  modes  that  propagate  adiabatically.    ...       74 

4-12  Comparison  of  fine  scale  coupled  mode  arrivals  using  24  m  (solid)  and 
48  m  (dotted)  step  sizes.  The  32  km  wave  guide  under  evaluation  was 
generated  from  range  propagated  WVLA  thermistor  records 78 

5-1  The  mean  (-),  median  (•••),  standard  deviation  (-•),  and  interquartile 
range  (-  -)  of  an  unscattered  mode  1  arrival  plotted  as  a  function  of 
noise  cutoff  level.  The  median  reaches  its  correct  value  of  zero  at  a 
noise  cutoff  level  of  -38  dB  relative  to  the  peak  of  the  arrival 82 
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5-2  The  solid  line  shows  the  mode  number  space  history  of  one  of  the  620 
subarrivals  that  would  make  up  the  mode  1  arrival  from  a  six  mode 
transmission  through  the  10  SI  A  soliton  waveguide  that  is  indicated 
by  the  dotted  lines 83 

5-3  Random  realizations  of  test  waveguides  with  from  10  to  35  200  m  SI  A 
solitons  over  32  km.  The  first  soliton  is  placed  randomly  about  2.5  km 
from  the  source,  and  the  remainder  are  added  randomly  from  there, 
with  an  average  of  450  m  between  leading  edges 86 

5-4  Single  realizations  of  mode  1  (solid),  mode  2  (dashed),  mode  3  (dot- 
dash),  and  mode  4  (dotted)  arrivals  for  0  to  35  200  m  SI  A  solitons  over 
32  km.  The  simulated  waveguides  used  are  similar  to  those  shown  in 
Figure  5-3.  The  simulated  source  depth  is  27  m  in  68  m  of  water,  and 
the  compressional  wave  attenuation  in  the  bottom  is  .046  dB/m  at  400 
Hz.  Note:  Mode  4  arrivals  are  not  visible  at  this  level  of  attenuation.        87 

5-5  The  same  scenario  as  Figure  5-4,  except  the  compressional  wave  atten- 
uation in  the  bottom  has  been  reduced  from  .046  to  0.0115  dB/m.  .  .       88 

5-6  Histograms  for  peak  position  of  mode  1  arrivals  for  5  through  40  ran- 
dom solitons.  The  solitons  start  about  2.5  km  from  the  source  with  an 
average  separation  of  450  m  between  soliton  leading  edges.  Time  zero 
corresponds  to  the  PAM1  arrival 89 

5-7  The  average  mean  (-),  median  (-  -),  and  peak  (•  •  •)  lag  with  respect  to 
the  pseudo-adiabatic  arrival,  and  standard  deviation  (-•)  of  scattered 
mode  1  arrivals  plotted  as  a  function  of  the  number  of  solitons  placed 
between  the  source  and  the  receiver,  where  the  solitons  are  added  start- 
ing from  the  source  end.  The  upper  three  lines  are  all  bias  measures, 
and  the  lower  line  is  a  measure  of  spread 90 
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5-8  The  average  mean  (-),  median  (-  -),  and  peak  (•  •  •)  lag  with  respect  to 
the  pseudo-adiabatic  arrival,  and  standard  deviation  (-•)  of  scattered 
mode  1  arrivals  plotted  as  a  function  of  the  number  of  solitons  placed 
between  the  source  and  the  receiver.  The  upper  curves,  averages  over 
50  realizations,  represent  the  case  where  the  solitons  are  added  starting 
at  the  receiver  end.  The  lower  set  of  curves  is  replotted  from  figure 
5-5.  60  solitons  is  essentially  saturation,  and  the  same  points  are  used 
in  both  sets  of  curves 91 

5-9  Time  dependent  waveguide  with  two  sets  of  six  SIA  solitons.  The 
upper  frame  shows  the  starting  positions,  the  middle  frame  shows  the 
positions  five  hours  later,  and  the  lower  frame  is  plot  of  peak  position 
bias  (re  PAM1)  as  a  function  of  time.  The  solitons  propagate  in  lock- 
step,  traversing  48  m  per  minute,  the  time  between  each  peak  position 
data  point  (o).  The  mean  arrival  data  is  also  in  one  minute  increments, 
but  the  points  are  connected  by  a  solid  line 96 

5-10  Waterfall  of  mode  1  arrival  patterns  from  the  time  dependent  waveg- 
uide with  two  sets  of  six  SIA  solitons.  The  300  peak  positions  in  the 
lower  frame  of  Figure  5-9  were  picked  from  these  arrival  patterns.     .  .       97 

5-11  Snapshots  of  the  arrival  shapes  for  all  eight  modes  propagated  in  the 
Figure  5-10  simulation.  The  left  hand  column  shows  the  minute  1  ar- 
rival shapes  and  the  right  column  shows  the  minute  300  arrival  shapes. 
The  thicker  lines  are  the  scattered  arrival  shapes,  and  the  thinner  lines 
show  the  corresponding  unscattered  mode  arrivals.  The  eighth  mode  is 
only  minimally  excited  by  the  source  at  27  m  depth,  and  it  undergoes 
the  most  attenuation,  so  it  is  not  even  visible  as  an  unscattered  arrival.    98 
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5-12  A  comparison  of  mean  (-)  and  peak  (o)  arrival  time  as  function  of 
the  number  of  evenly  spaced  solitons  in  a  32  km  waveguide.  The  first 
soliton  is  always  200  m  from  the  source  and  the  last  soliton  is  roughly 
32  km  divided  by  the  number  of  solitons  away  from  the  receiver.  ...     100 

5-13  The  top  graph  shows  12.5  m  depth  temperature  readings  from  the  2 
hrs  of  thermistor  data  used  to  simulate  modal  excitation  levels.  The 
middle  graph  shows  the  resulting  modal  excitation  levels  for  mode  1 
(solid),  mode  2  (dashed),  mode  3  (dot-dash),  and  mode  4  (dotted). 
Mode  5  is  similar  in  magnitude  and  used  as  an  input  to  the  model,  but 
it  is  omitted  from  this  figure  to  avoid  clutter.  The  lower  graph  shows 
the  output  bias  levels  corresponding  to  the  excitation  inputs  shown  in 
the  middle  graph.  The  solid  line  is  the  mode  1  arrival  envelope  peak 
arrival  time  in  excess  of  the  PAM1  arrival  time  and  the  dotted  line  is 
the  mode  1  arrival  envelope  mean  arrival  time  in  excess  of  the  PAM1 
arrival  time 101 

5-14  Pycnocline  displacements  from  an  Orr  backscatter  record  compared  to 
a  thermistor  string  record  using  a  propagation  speed  of  .62  m/s  over 
the  distance  between  the  back  scatter  instrument  and  the  thermistor 
string  at  each  backscatter  reading 103 

5-15  Distribution  of  peaks  for  the  mode  1  arrivals  of  simulated  transmissions 
using  sequential  WVLA  thermistor  string  generated  waveguides.  The 
sequential  32160  by  68  m  waveguides  are  separated  by  six  minutes. 
The  size  of  each  arrival  dot  is  scaled  by  the  amplitude  of  the  peak 
arrival.  The  upper  frame  includes  a  plot  of  the  PAM1  arrival  time  in 
conjunction  with  the  individual  peak  arrival  times;  the  zero  reference 
on  the  arrival  time  axis  corresponds  to  21.6  seconds.  The  lower  frame 
plots  the  bias  or  offset  of  the  peak  arrivals  relative  to  the  PAM1  arrival 
time  for  each  realization 105 
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5-16  Distribution  of  peaks  for  the  mode  2  arrivals  corresponding  to  the  mode 
1  arrivals  of  Figure  5-15.  The  upper  frame  includes  a  plot  of  both  the 
PAM1  arrival  time  (dotted  line)  and  PAM2  arrival  time  (solid  line),  in 
conjunction  with  the  individual  peak  arrival  times;  the  zero  reference 
on  the  arrival  time  axis  corresponds  to  21.6  seconds.  The  lower  frame 
plots  the  bias  or  offset  of  the  mode  2  peak  arrivals  relative  to  the  PAM2 
arrival  time  for  each  realization 106 

5-17  30  minute  bin  averaged  IQR  measurements  and  the  WVLA  tempera- 
ture profile  for  the  simulation  test  period.  The  profiled  temperatures 
range  from  7.2  to  25.9  degrees  C.  The  top  frame  is  for  mode  1  and  the 
bottom  is  for  mode2 108 

5-18  A  comparison  of  1  hour  bin  averaged  IQR  levels  for  the  simulated 
mode  2-4  data  (solid  lines)  and  a  1  hr  bin  average  of  the  Standard 
Deviation  of  the  Temperature  at  22.5  m  depth  over  the  previous  4.3 
hrs  (dashed  lines).  The  standard  deviation  of  the  temperature  is  in 
tenths  of  a  degree  centigrade  to  facilitate  comparison  with  the  IQR 
levels  in  milliseconds 110 

5-19  Cross- correlations  of  1.4  hr  averaged  peak  height  fluctuations  for  mode 
1  (solid),  2  (dashed),  3  (dot-dashed),  and  4  (dotted)  with  the  corre- 
sponding 1.4  hr  averaged  mode  1  IQR  (spread)  levels Ill 

6-1  Distribution  of  SWARM  Mode  1  Peak  arrival  times.  The  zero  refer- 
ence is  at  an  arrival  time  of  21.985  seconds.  The  plotted  transmissions 
occurred  every  six  minutes,  with  each  transmission  producing  22  se- 
quence arrivals  over  a  two  minute  period.  The  location  of  the  highest 
peak  in  each  sequence  is  plotted.  Though  not  evident  in  most  cases, 
there  are  gaps  in  the  record  that  range  from  single  lost  transmissions 
to  occasional  larger  gaps  of  up  to  two  hours  in  duration 116 
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6-2  Predicted  current  effects  on  PAM1  arrival  time  plotted  with  the  dis- 
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based  on  the  mean  current  between  39  and  59  m  depth  and  the  dotted 
curve  is  based  on  the  average  current  between  11  and  67  m  depth.   .  .     118 

6-3  Predicted  temperature  fluctuation  effects  on  PAM1  arrival  time  (solid 
curve)  plotted  with  the  distribution  of  "current- corrected"  peak  arrival 
times.  The  solid  curve  is  zero  mean  (by  definition)  and  offset  from  the 
data  for  convenience 119 

6-4  Mode  1  peak  arrival  time  bias  histograms  for  the  data  shown  in  Figure 
6-1  and  the  model  distribution  shown  Figure  5-15.  The  bias  measure- 
ments are  relative  to  estimated  PAM1  arrival  times  (i.e.  the  leading 
edge)  for  the  data  and  calculated  PAM1  arrival  times  for  the  model.  .     120 

6-5  Mode  1  arrival  spread,  as  measured  by  the  IQR  of  the  arrival  enve- 
lope (unscattered  IQR  is  about  7.2  ms).  The  upper  frame  is  plot  of 
the  spread  when  averaged  over  each  transmission  (22  sequences  per 
transmission).  The  lower  frame  shows  bin  averaged  spreads  for  five 
bin  widths  ranging  from  1.4  to  4.1  hrs 122 

6-6     Power  Spectral  Density  of  the  bin  averaged  spreads  shown  in  the  lower 

frame  of  Figure  6-5 124 

6-7  WVLA  temperature  as  function  of  time  and  depth  overlaid  by  the  1.4 
hr  bin  averaged  IQR  of  mode  1  arrivals  at  the  WVLA.  The  tempera- 
tures in  this  figure  range  from  7.2  to  26.5  degrees  C 125 

6-8  A  comparison  of  1  hour  bin  averaged  IQR  levels  for  the  SWARM 
data  (solid  line),  the  corresponding  propagated  thermistor  string  model 
(dotted  line),  and  the  Standard  Deviation  of  the  Temperature  at  22.5  m 
depth  over  the  previous  4.3  hrs  (dashed  line).  The  standard  deviation 
of  the  temperature  is  in  tenths  of  a  degree  C  to  facilitate  comparison 
with  the  IQR  levels  in  milliseconds 127 
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6-9  A  comparison  of  1.4  hr  averaged  peak  intensity  fluctuations  and  the 
corresponding  1.4  hr  averaged  IQR  levels  for  the  SWARM  data.  The 
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6-10  A  comparison  of  1.4  hr  averaged  peak-height  fluctuations  for  the  first 
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mode  3  is  dashed,  and  mode  4  is  dotted 130 

6-11  Cross-correlations  of  1.4  hr  averaged  peak  height  fluctuations  for  mode 
1  (solid),  2  (dashed),  3  (dot-dashed),  and  4  (dotted)  with  the  corre- 
sponding 1.4  hr  averaged  mode  1  IQR  (spread)  levels  for  the  SWARM 
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6-12  Mode  1  (solid),  mode  2  (dashed),  mode  3  (dot-dashed),  and  Mode  4 
(dotted)  peak  arrival  time  leading  edge  (connects  earliest  peak  arrivals 
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6-13  Distribution  of  SWARM  224  Hz  Mode  1  Peak  arrival  times.  The  zero 
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The  later  zero  reference  time,  relative  to  Figure  6-1,  stems  from  the 
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Chapter  1 


Background 


The  Shallow  Water  Acoustics  in  a  Random  Medium  (SWARM)  experiment  [1],  con- 
ducted off  the  coast  of  New  Jersey  near  the  continental  shelfbreak  in  July  and  August 
of  1995,  was  dedicated  to  furthering  our  understanding  of  acoustic  mode  scattering 
effects  of  both  linear  internal  waves  and  nonlinear  solitary  internal  waves  (solitons) 
in  shallow  water. 


1.1      Solitons 

A  stratified  ocean  will  support  a  variety  of  internal  waves,  the  most  common  being  the 
internal  equivalent  of  a  linear  surface  gravity  wave.  In  the  deep  ocean,  these  waves  are 
generally  isotropic  in  direction  and  have  frequency/vertical  mode  power  spectra  that 
are  red  and  dominated  by  the  first  few  modes.  The  canonical  Garrett-Munk  (GM) 
spectrum  [3]  normally  provides  a  good  model  for  the  observed  power  spectra  in  deep 
water.  The  story  is  quite  different  along  the  United  States  East  Coast  continental 
shelf  and  other  shallow  regions,  where  the  linear  waves  are  non-isotropic  and  the 
internal  wave  field  can  also  be  heavily  influenced  by  nonlinear  solitary  waves. 

Internal  solitons  appear  to  be  generated  in  groups  or  packets  at  the  continental 
shelfbreak  when  the  tide  shifts  from  an  ebb  to  a  flood.  Each  "developed"  packet  usu- 
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ally  contains  six  to  twelve  solitons  that  are  large  enough  to  have  a  surface  expression 
on  radar,  so  the  term  "solitary  wave"  is  something  of  a  misnomer  in  this  case.  An 
account  of  the  first  observation  of  a  surface  solitary  wave,  made  by  John  Scott  Russell 
in  1834,  is  provided  by  Strang  [4]: 

I  was  observing  the  motion  of  a  boat  which  was  rapidly  drawn  along  a 
narrow  channel  by  a  pair  of  horses,  when  the  boat  suddenly  stopped. . .  not 
so  the  mass  of  water  in  the  channel  which  it  had  put  in  motion;  it  accu- 
mulated round  the  prow  of  the  vessel  in  a  state  of  violent  agitation,  then 
suddenly  leaving  it  behind,  rolled  forward  with  great  velocity,  assuming 
the  form  of  a  large  solitary  elevation,  a  rounded,  smooth  and  well-defined 
heap  of  water,  which  continued  its  course  along  the  channel  apparently 
without  change  of  form  or  diminution  of  speed.  1  followed  it  on  horseback, 
and  overtook  it  still  rolling  at  a  rate  of  some  eight  to  ten  miles  an  hour, 
preserving  its  original  figure  some  thirty  feet  long  and  a  foot  to  a  foot  and 
a  half  in  height.  Its  height  gradually  diminished,  and  after  a  chase  of  one 
or  two  miles  I  lost  it  in  the  windings  of  the  channel. 

This  surface  soliton  can  be  modeled  by  the  Korteweg-de  Vries  (KdV)  equation, 

Vt  +  6Wx  -r  rjxxx  =  0,  (1.1) 

where  the  wave  solution, 

r,(x  -  ct)  =  ^sech2  (- yfc{x  -  d))  ,  (1.2) 

has  a  single  hump  whose  height,  |c,  is  proportional  to  the  speed  at  which  it  travels 
to  the  right  [4]. 

In  contrast  to  the  surface  soliton's  hump,  internal  solitons  normally  appear  as  a 
sech2  shaped  depression  in  the  pycnocline.  In  covering  some  theoretical  aspects  of 
internal  solitons,  Apel  et  al.    [5]  provide  a  solution  for  the  following  internal  wave 
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KdV  equation, 

Vt  +  cr]x  +  ocr}T)x  +  /3r]xxx  =  0,  (1.3) 

namely, 

rj(x,t)  =  r/0sech2  (-  -      j,  (1.4) 

where  the  soliton  phase  speed,  V,  is, 

V  =  c  +  ar70/3  (1.5) 

and  the  characteristic  width,  A,  is  given  by, 

A  2       12/3  ,      x 

A  =  ™  ■  (L6^ 

For  a  simple  two-layer  ocean  density  model,  the  constants  are 


c  = 


1i 
2 

_        3c       ( P2h2,-pihl 


(P2hj-Pihl\  Q   -jA 

2/ll/l2     \P2/ll+Pl/l2/  V     '     ' 

a  _  c/i]/i2   /  pi  /11  -1-^2/12^ 
^  6        \P2hi+p\h2) 

Where  h\  and  /12  are  the  thicknesses  of  the  upper  and  lower  layers,  and  p\  and  p<i  are 
the  respective  densities.  From  the  characteristic  width  formulation,  we  see  that  770 
will  have  the  same  sign  as  a,  so  for  the  most  common  case  of  a  shallow  pycnocline  over 
a  deeper  layer  (i.e.  p2h\  <  pi/^),  V  will  be  negative  (i.e.  downgoing).  Employing  this 
model  with  reasonable  2-layer  approximations  for  a  large  SWARM  soliton,  770  =  20 
m,  hi  =  15  m,  h2  =  55  m,  and  gAp/p  =  .01  m/s2,  we  predict  a  rather  typical  soliton 
phase  speed  of  V  =  .8  m/s  and  the  wave  shape  shown  in  Figure  1-1.  The  linear  phase 
speed,  c,  for  this  case  is  0.34  m/s  for  comparison. 

The  internal  solitary  wave  packets  are  often  so  numerous  and  closely  spaced  on  the 
New  Jersey  shelf  in  summer  that  a  nearly  continuous  train  of  solitons  can  pass  a  fixed 
instrument  over  the  course  of  several  hours.  This  is  illustrated  in  Figure  1-2,  where  an 
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Figure  1-1:  A  simple  two-layer  model  of  an  internal  soliton,  with  the  upper  layer, 
hi  =  15  m,  lower  layer,  &2  =  55  m,  and  reduced  gravity,  gAp/p  =  .01  m/s2.  The  20 
m  depression  propagates  at  a  phase  speed,  V  =  .8  m/s. 

array  of  thermistors  attached  to  the  WHOI  Vertical  Line  Array  (WVLA)  records  the 
passage  of  several  packets  over  a  six  hour  period.  The  mixture  of  time  and  space  in 
this  figure  creates  a  distorted  view  of  the  soliton  shapes,  but  this  can  be  approximately 
corrected  by  applying  a  constant  phase  speed  to  a  selected  soliton,  thus  transforming 
the  time  coordinate  into  a  space  coordinate.  This  has  been  accomplished  in  Figure 
1-3,  where  the  time  series  for  the  prominent  soliton  in  the  lower  frame  of  Figure  1-2 
has  been  transformed  into  a  properly  scaled  space  representation. 

1.2      Acoustic  Impact  of  Solitons 

The  first  investigation  of  the  acoustic  impact  of  solitons  was  reported  by  Zhou  et  al. 
[6]  in  1991.  This  paper  generated  interest  in  the  subject,  and  further  strides  have 
been  made  in  the  several  investigations  that  have  followed  [1,  2,  7,  8,  9]. 


1.2.1      Time  Spreading 

Internal  solitons  often  have  first  mode  displacements  at  the  pycnocline  ranging  from 
10  to  20  meters.  These  large  displacements  occur  over  relatively  short  wavelengths 
(Figure  1-3),  creating  strong  horizontal  gradients  in  soundspeed,  which  thus  produces 
coupling  between  acoustic  normal  modes.  Differences  in  the  acoustic  modal  group 
velocities  cause  time  spreading  of  pulses  as  energy  is  transferred  between  different 
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Figure  1-2:  Time  and  space  interpolation  of  four  thermistor  records.  The  thermistors 
were  attached  to  the  WVLA  at  water  depths  of  12.5,  20.5,  30.5,  and  40.5  meters  and 
sampled  at  30  second  intervals. 
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Figure  1-3:  Space  transformation  of  the  time  series  for  the  prominent  soliton  in  the 
lower  frame  of  Figure  1-2.  All  six  available  thermistors  are  used,  so  the  depth  spans 
12.5  to  60.5  meters.  The  dimensions  shown  are  properly  scaled  for  soliton  with  a 
combined  phase  and  advection  speed  of  .8  m/s. 
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modes.  In  the  SWARM  waveguide,  mode  1  has  the  highest  group  speed,  so  the 
energy  that  transits  in  other  modes  for  a  portion  of  the  trip  is  delayed  with  respect  to 
energy  that  traveled  exclusively  in  mode  1 .  This  causes  the  spread  of  mode  1  energy 
to  be  all  in  one  direction  (i.e.  there  is  a  positive  bias  in  the  mean  and  median  travel- 
time  of  mode  1  energy  relative  to  an  unscattered  mode  1  arrival  time).  Higher  modes 
may  have  some  energy  coming  in  earlier  (from  traveling  partially  in  lower  modes) 
and  some  coming  in  later  (from  traveling  partially  in  higher  modes),  so  there  may 
be  negative,  positive,  or  even  zero  bias  in  travel-time  for  these  modes.  As  shown  by 
Colosi  and  Flatte  [10]  the  coupling  caused  by  linear  internal  waves  is  similar,  but 
reversed,  in  deep  water  where  mode  1  generally  has  the  lowest  group  speed. 

Repeated  coupling  processes  create  mode  arrivals  that  are  sums  of  numerous 
"mode  multipaths"  that  propagate  through  the  range  and  wavenumber  (mode  num- 
ber) space  between  the  source  and  the  receiver.  Constructive  and  destructive  interfer- 
ence between  "mode  multipaths"  creates  peaks  and  nulls  in  the  arrival  pattern.  The 
detailed  interference  patterns  are  very  sensitive  to  the  spacing  between  the  coupling 
events.  Thus,  given  the  tidal  period  variations  in  the  size  and  numbers  of  solitons 
between  the  source  and  receiver  and  the  continuous  variations  in  soliton  position  and 
spacing  (caused  by  their  differing  phase  speeds),  significant  temporal  variations  in 
arrival-time  patterns  are  to  be  expected. 


1.2.2     Attenuation 

Another  effect  of  the  coupling  of  mode  1  into  higher  modes  is  the  apparent  increased 
attenuation  of  mode  1  discussed  by  Zhou  et  al.  [6].  In  general  modal  attenuation 
is  greater  for  higher  mode  numbers,  so  coupling  is  an  energy  sink  for  low  modes 
and  a  source  for  higher  modes.  This  mode-stripping  along  with  long-term  variations 
in  the  amount  of  coupling  and  short-term  variations  in  interference  patterns  creates 
significant  variability  in  the  apparent  attenuation  of  mode  1  energy. 
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1.3      Outline  of  Thesis 

Chapter  2  discusses  SWARM  data  that  are  relevant  to  this  thesis,  along  with  other 
details  such  as  experiment  geometry,  signal  characteristics,  data  storage  requirements, 
and  pulse- compression  schemes.  Several  issues  associated  with  mode-filtering  are 
then  addressed  in  Chapter  3:  a  basic  theory  of  normal  mode  propagation;  analysis  of 
possible  errors  in  the  mode  filters  used;  and  the  accuracy  of  the  filter  outputs  relative 
to  the  actual  mode  arrivals. 

A  model  for  the  propagation  of  normal  modes  through  a  2-D  waveguide  with 
solitons  is  developed  in  Chapter  4.  Chapter  5  discusses  pulse  propagation  simulations 
using  two  different  forms  of  the  model.  The  simplest  form  of  the  model,  which 
uses  Preisig  and  Duda's  [2]  sharp  interface  approximation  (an  inverted  square-wave) 
for  solitons,  was  tested  by  varying  the  numbers  and  locations  of  solitons  and  the 
amount  of  bottom  attenuation.  In  observing  the  sensitivity  of  mode  arrival  statistics 
to  changes  in  the  model  parameters,  candidate  statistics  for  quantifying  the  spread 
and  bias  of  mode  arrivals  were  evaluated.  The  second,  more  complicated,  form  of  the 
model  is  based  on  range  propagated  thermistor  string  records,  so  it  models  the  effects 
of  both  linear  and  nonlinear  internal  waves.  It  is  used  to  try  to  model  the  temporal 
variations  in  the  SWARM  waveguide  during  the  time  of  the  actual  experiment. 

Chapter  6  discusses  the  analysis  of  mode  1-4  arrivals  for  the  400  Hz  source  and 
the  mode  1-3  arrivals  for  the  224  Hz  source.  The  narrow  pulse  width  of  the  400  Hz 
source  provided  much  better  data,  and  it  receives  a  greater  share  of  our  attention. 
Chapter  7  lists  the  conclusions  reached  and  proposes  outlines  of  future  investigation. 
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Chapter  2 


Data  Collection 


2.1      Experiment  Geometry 

The  sixteen  element  WVLA,  moored  in  70.5  meters  of  water  off  the  New  Jersey  coast, 
operated  continuously  for  seven  days.  Tomography  sources  with  carrier  frequencies 
of  224  and  400  Hz  were  moored  about  32  km  shoreward  on  a  compass  heading  of  300° 
and  transmitted  to  the  WVLA.  The  resulting  acoustic  path  was  anti-parallel  to  the 
primary  propagation  direction  for  shelf  generated  internal  waves.  The  temperature 
profile  in  the  vicinity  of  the  WVLA  was  monitored  by  six  thermistors  attached  to  the 
array  at  various  depths.  Figure  2-1  illustrates  the  2-D  geometry  of  the  experiment. 


2.2      Source  Signals 

The  400  Hz  (100  Hz  bandwidth)  tomography  source  produced  a  511  digit,  phase- 
modulated  signal  sequence  with  a  digit  length  of  4  cycles,  for  a  total  of  5.11  seconds 
per  sequence.  The  sequence  was  repeated  23  times  during  each  transmission,  for  a 
transmission  time  of  117.53  seconds.  The  400  Hz  transmissions  started  on  the  hour 
and  were  repeated  every  six  minutes.  They  were  received  at  the  WVLA  about  22 
seconds  later. 
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Figure  2-1:  Source  positions  are  noted  by  "*".  The  400  Hz  source  is  the  shallower  of 
the  two.  WVLA  hydrophone  positions  are  noted  by  "o"  and  the  thermistors  on  the 
array  are  marked  as  "+". 

The  224  Hz  (16  Hz  bandwidth)  tomography  source  produced  a  63  digit  sequence, 
phase-modulated  signal  with  a  digit  length  of  14  cycles,  for  a  total  of  3.9375  seconds 
per  sequence.  The  sequence  was  repeated  30  times  during  each  transmission,  for 
a  transmission  time  of  118.125  seconds.  The  transmissions  were  repeated  every  five 
minutes,  starting  on  the  hour,  five  after,  etc.  They  were  received  at  the  WVLA  about 
22.5  seconds  later. 


2.3     Receiving  Array 

The  WVLA  had  hydrophones  spanning  water  depths  from  14.9  to  67.4  m.  The 
hydrophones  could  not  be  placed  any  higher  in  the  water  column  without  experiencing 
surface  wave  effects.  The  sixteen  hydrophones  were  sampled  at  1395.089286  Hz  six  of 
the  seven  days  the  array  was  deployed.  The  sample  rate  for  the  remaining  period  was 
1953.1250  Hz.  The  data  were  normally  telemetered  back  to  a  moored  research  vessel 
(R/V  Oceanus)  stationed  several  kilometers  away,  but  a  limited  amount  of  internal 
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Tape 

Number 

First 

Last 

Number 

of  Files 

File 

File 

1 

16 

210-19:54 

211-10:48 

2 

7 

211-16:43 

211-22:58 

3 

2 

212-01:09 

212-01:11 

4 

1 

212-01:50 

212-01:50 

5 

14 

212-03:27 

212-17:01 

6 

32 

213-02:02 

215-06:01 

7 

29 

215-06:27 

216-11:41 

8 

18 

216-13:04 

217-06:49 

9 

34 

210-19:50 

214-05:13 

10 

30 

214-00:34 

217-17:50 

Table  2.1:  Data  set  outline  for  the  WVLA.  Tapes  1  through  8  contain  telemetered 
data  that  is  sometimes  broken  up  when  the  ship  was  out  of  range.  Tapes  9  and  10 
contain  internally  recorded  data  that  fills  in  most  of  the  gaps  in  the  telemetered  data. 


storage  was  available  and  automatically  utilized  when  the  ship  was  out  of  range.  The 
non-integer  sample  rates  were  dictated  by  the  system's  analog  to  digital  converter. 
The  equipment  provided  a  very  accurate  sample  rate,  but  the  sample  frequencies  were 
limited  to 


Sample  Rate 


5MHz 


(2.1) 


256  *  ni  *  2n2 
In  our  case,  rii  =  7  or  5,  and  n2  =  1. 

The  data  were  stored  on  tapes,  with  about  27  files  per  tape.  Each  file  contains 
up  to  160  records  which  are  approximately  1  megabyte  each.  The  records  consist  of 
32768  two  byte  digitized  voltage  samples  for  each  of  the  sixteen  hydrophones.  The 
sample  rate  of  1395.089286  Hz,  yields  about  23  seconds  of  data  per  record  and  62 
minutes  of  data  per  file.  Seven  days  of  data  required  over  25  gigabytes  of  storage, 
and  an  outline  of  the  entire  data  set  is  given  in  Table  2-1. 

The  combination  of  the  1395  Hz  sample  rate  and  sequence  length  for  the  400 
Hz  source  yields  7128.9063  samples  per  sequence,  and  the  224  Hz  arrivals  require 
5493.1641  samples  per  sequence.   Both  of  these  numbers  are  rather  cumbersome  to 
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deal  with.  An  ideal  combination  of  sequence  length  and  sample  rate  would  produce  a 
radix-2  integer  number  of  samples  per  sequence.  During  the  period  from  about  1600Z 
day  213  to  1600Z  on  day  214,  the  sample  rate  was  shifted  from  1395  to  1953  Hz,  and 
the  resulting  files  only  cover  44  minutes.  Seven  transmissions  are  removed  from  each 
160  Mb  file  during  this  period. 

There  are  two  further  data  considerations.  First  is  the  location  of  the  header  on 
each  file.  The  telemetered  data  files  have  headers  at  the  beginning  of  the  files.  The 
internally  recorded  data  have  their  headers  are  at  the  end.  Second,  as  a  consequence 
of  the  shifts  between  telemetered  and  internally  recorded  data,  some  files  with  less 
than  160  Mb  were  created  during  both  the  1395  and  1953  Hz  sampling  periods.  These 
files  obviously  contain  fewer  transmissions. 

2.4      Signal  Isolation 
2.4.1      Sequence  Capture 

The  400  Hz  arrivals  at  the  WVLA  can  be  analyzed  one  sequence  at  a  time,  but  of  the 
23  sequences  in  a  transmission,  only  22  full  sequences  can  be  captured  for  analysis. 
The  loss  of  one  sequence  is  required  to  compensate  for  the  uncertainty  in  the  actual 
arrival  time  of  the  beginning  of  the  first  sequence  in  a  transmission.  The  transmitted 
sequences  have  an  ideal  circular  autocorrelation  property  that  is  unaffected  by  the 
starting  point  in  a  captured  sequence. 

In  our  case,  sequence  capture  is  accomplished  by  starting  the  first  capture  at  a 
fixed  time  (26.85  sees)  after  the  beginning  the  transmission.  The  fixed  time  for  the 
first  sample  was  chosen  through  a  series  of  trials  to  be  about  400  samples  (depending 
on  the  particular  transmission)  or  280  msec  prior  to  the  start  of  the  arrival  of  the 
second  sequence.  The  earlier  parts  of  the  first  sequence  are  discarded,  and  the  last 
400  samples  are  clipped  from  the  second  sequence  to  provide  the  data  buffer  for  the 
beginning  of  the  third  sequence.  This  process  continues  through  the  sequences  until 
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the  last  400  samples  of  the  23rd  sequence  are  clipped  and  discarded.  The  22  sequences 
obtained  in  this  manner  are  rounded  to  7129  samples  long,  and  have  basically  been 
circularly  shifted  to  the  right  (i.e.  delayed)  by  about  400  samples  with  respect  to  the 
transmitted  sequence.  In  a  similar  manner,  29  sequences  of  the  30  transmitted  can 
be  captured  from  each  224  Hz  transmission.  The  224  Hz  sequences  are  rounded  to 
5493  samples  for  analysis. 

Changes  in  travel-time  are  estimated  by  capturing  sequences  using  a  fixed  time 
delay  after  transmission  and  detecting  the  resulting  circular  shifts  in  the  captured 
sequences.  There  are  a  non-integer  number  of  samples  between  transmissions,  so 
choosing  the  first  sample  to  capture  in  a  sequence  creates  an  error  between  the  desired 
fractional  sample  number  and  the  allowable  integer  sample  numbers.  A  correction  for 
the  time  and  phase  error  introduced  by  the  sample  error  can  be  made  by  multiplying 
the  base-band  signal  in  the  frequency  domain  by  et2^^+^At,  where  /  is  the  base-band 
frequency,  fc  is  the  carrier  frequency,  and  At  is  the  sample  error  divided  by  the  sample 
rate.  This  correction  step  is  accomplished  in  conjunction  with  pulse-compression. 

2.4.2      Pulse  Compression 

The  captured  400  Hz  sequences  are  demodulated  to  baseband  and  then  low-pass 
filtered  by  windowing  in  the  frequency  domain  with  a  90  Hz  pass-band  Blackman 
window.  Pulse  compression  is  also  accomplished  in  the  frequency  domain  through 
multiplication  with  the  conjugate  spectrum  of  the  transmitted  sequence.  This  is 
equivalent  to  a  matched-filter  circular  convolution  in  the  time  domain  [11]. 

If  sequence  lengths  of  5493  and  7129  are  maintained,  mixed-radix  DFTs  must  be 
used  in  the  transforms  between  the  time  and  frequency  domains.  These  mixed-radix 
DFTs  are  too  slow  for  practical  application,  so  the  data  was  zero-padded  to  a  length 
of  8192  in  order  to  employ  a  radix-2  FFT. 

The  circular  convolutions  of  delayed  and  zero-padded  signals  with  zero-padded 
transmitted  sequences  are  not  ideal.  The  resulting  sidelobes  (Figure  2-2)  stem  from 
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Figure  2-2:  Circular  convolution  of  an  ideal  delayed  400  Hz  transmitted  sequence  with 
an  ideal  non-delayed  transmitted  sequence.  Both  sequences  have  been  zero-padded 
from  length  7129  to  length  8192.  The  lower  frame  is  a  window  of  the  full  output 
above. 

offsets  in  the  location  of  the  zeros.  The  sidelobe  effects  were  eliminated  from  all 
analyses  presented  through  SNR  discrimination. 

It  is  not  necessary  to  save  the  entire  output  of  the  pulse  compressed  data.  The 
relevant  portions  of  the  data  are  limited  to  an  area  around  the  peak  arrival  that 
encompasses  the  spread  of  energy  caused  by  mode  dispersion  and  coupling.  A  200 
sample  (143  msec)  observation  window  like  that  shown  in  Figure  2-2  is  sufficient  to 
capture  the  spread  expected  from  modal  dispersion  over  32  km,  but  an  additional 
70  msec  is  required  to  accommodate  the  wander.  Figure  3-3  shows  some  typical 
saved  outputs  of  400  Hz  pulse-compressed  arrival  patterns  for  the  hydrophones  at  the 
WVLA.  Due  to  its  wider  pulse  width,  358  msec  are  saved  for  the  224  Hz  arrivals. 
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Figure  2-3:  Typical  400  Hz  arrivals  at  the  WVLA  hydrophones.  The  sets  of  16  arrivals 
are  separated  by  6  minutes. 
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Chapter  3 


Mode  Filtering 


3.1      Normal  Modes 

The  two  most  frequently  used  representations  of  sound  propagation  in  ocean  acoustic 
waveguides  are  rays  and  normal-modes.  Ray  theory  is  perhaps  more  intuitive,  but 
it  is  a  high  frequency  approximation  that  is  better  suited  to  situations  where  the 
acoustic  wavelength  is  small  compared  to  the  dimensions  of  the  sound  channel.  Our 
low  frequencies  and  shallow  waveguide  are  better  suited  to  a  modal  analysis  than  ray 
theory. 

The  modal  analysis  of  400  Hz  acoustic  arrivals  at  the  WVLA  begins  by  treating 
the  acoustic  signal  at  each  hydrophone  as  a  sum  of  vertical  modes 

N 
Ph{t)  =  Y,An{t)<l>n{Zh),  (3.1) 

71=1 

where  <bn{z)  is  the  mode  shape,  An(t)  is  the  mode  coefficient,  and  z^  is  the  depth  of  the 
hydrophone  in  question.  The  number  of  propagating  normal  modes,  N,  as  well  as  the 
mode  shapes  themselves  are  frequency  dependent.  The  mode  shapes  are  analogous 
to  the  eigenmodes  of  a  vibrating  string,  and  are  formed  by  a  constructive  vertical 
interference  pattern  caused  by  a  pair  of  up  and  down  going  plane-waves  with  the 

41 


same  horizontal  wave  number.  The  highest  possible  normal  mode  number  is  roughly 
equivalent  to  the  sound  channel  depth  divided  by  the  acoustic  wavelength,  where  the 
up  and  down  going  pair  of  plane- waves  would  be  propagating  with  a  horizontal  wave 
number  approaching  zero. 

The  six  thermistors  on  the  WVLA,  which  are  sampled  every  30  seconds,  provided 
the  soundspeed  profiles  needed  for  computing  a  time-dependent  set  of  mode  shapes 
using  the  KRAKEN  [12]  normal-mode  code.  The  basis  for  this  code,  a  numerical 
normal-mode  solution  of  the  Helmholtz  equation  at  frequency  u>, 

[v2  +  k2{x,y,z,u)\P{x,y,z,u;)  =  0,  (3.2) 

where  p(r,  uj)  is  the  Fourier  component  of  the  pressure  field,  can  be  found  in  Jensen 
et  al.  [13]. 

The  KRAKEN  [12]  mode  shapes  are  orthonormal  in  the  sense  that, 

/*-(yu^,  0.3) 

J  p{z) 

so  given  p(z,t),  An(t)  can  be  readily  computed  by  applying  the  operator 


/(•)^W  (3.4) 


Since  the  pressure  field,  p(z,t),  is  only  sampled  at  the  sixteen  hydrophone  locations, 
the  integral  computation  of  the  mode  coefficients,  An(t),  must  be  approximated  by 

An{t)c,Y:p{z^\Zh)Az,  (3.5) 

where  Az  is  the  spacing  between  hydrophones.  This  discrete  sum  can  be  easily  im- 
plemented through  matrix  multiplication  of  mode  shape  vectors  with  16  channel  hy- 
drophone data  vectors.  This  process  of  converting  the  16  channels  of  hydrophone  data 
into  N  normal  mode  arrivals  is  referred  to  as  narrowband  direct-projection  (NBDP) 
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mode-filtering. 

An  alternative,  broadband  direct-projection  (BBDP)  mode-filtering  operates  in 
the  frequency  domain,  approximating  the  Fourier  components  of  the  mode  arrivals  as 

AMl)  *  E  E^m^dj^l,  (3.6) 

where 

P(*ft  Jj)  =  EK^gc-^»At.  (3.7) 

m 

The  time  domain  mode  coefficients  can  then  be  recovered  through 

An(tm)  =  £  AntfJtpiV't-Af.  (3.8) 

i 

The  benefits  of  this  method  over  NBDP  depend  upon  the  frequency  dependence 
of  the  mode  shapes  over  the  bandwidth  of  the  arriving  signal.  The  cost  of  BBDP  is 
in  computation  time;  over  4000  different  sets  of  narrowband  mode  filters  are  required 
for  each  frequency  to  cover  the  time-dependence  seen  in  six  days  of  data,  so  in  the 
ideal  extreme  of  a  different  set  of  <f>n(fj)s  for  each  of  the  1000  frequency  bins  that 
cover  the  bandwidth  of  the  signal,  4  million  sets  of  mode  filters  would  be  have  to  be 
computed!  The  extreme  case  is  not  feasible,  but  there  are  compromises  that  can  be 
explored. 

Either  DP  method  should  work  very  well  for  the  low  modes  that  are  adequately 
sampled  by  the  WVLA.  A  typical  SWARM  mode  1  shape  will  yield, 

f  *(**)*(«*)*«  =  0.9989,  (3.9) 

S        pw 

very  close  to  the  ideal  output  of  1.0.  Higher  modes  have  more  energy  outside  the 
depths  sampled  by  the  WVLA,  so  the  above  "auto-filtration"  degrades,  yielding  out- 
puts of  0.9949,  0.9881,  0.9803,  0.9681,  0.9269,  and  0.8243  for  modes  2  through  7, 
respectively.  The  higher  modes  suffer  more  from  partial  coverage  than  they  do  from 
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undersampling,  so  a  simple  output  gain  factor  (like  dividing  the  filter  output  by  its 
auto-filtration  level)  may  be  useful  if  level  comparisons  between  mode  arrivals  are  to 
be  made.  However,  an  output  gain  factor  will  not  provide  relief  from  modal-crosstalk, 
which  is  also  inherent  with  undersampling  and  partial  coverage. 

One  must  keep  in  mind  that  mode  coefficient  outputs  of  the  mode-filtration  are 
dependent  on  the  quality  of  both  the  computed  mode  shapes  and  the  measured  pres- 
sure field.  The  most  significant  error  associated  with  pressure  field  measurements  is 
uncompensated  array  tilt.  Deviations  from  a  perfectly  vertical  array  will  cause  phase 
differences  between  hydrophones  that  degrade  the  mode  coefficient  computations. 


3.2      Array  Tilt 

The  designed  navigation  system  for  the  WVLA  was  extensive  and  complete,  with 
three  transponder  travel-times  recorded  at  four  hydrophones  on  the  array.  Unfor- 
tunately, the  travel-times  from  one  of  the  transponders  to  three  of  the  hydrophones 
did  not  record  properly,  leaving  just  one  channel  (the  shallowest  hydrophone)  fully 
navigated.  Moreover,  all  of  the  navigation  travel-times  appear  to  be  heavily  influ- 
enced by  small-scale  temperature  fluctuations  caused  by  internal  waves.  The  waves 
caused  navigation  travel- time  changes  that  were  three  to  five  times  larger  than  those 
expected  from  mooring  motion.  The  time  scale  of  these  travel-time  changes  is  on  the 
order  of  the  navigation  sample  rate  (once  every  four  minutes),  consistent  with  soliton 
passage.  The  longer  time-scale  navigation  travel-time  changes  are  similarly  out  of 
proportion  with  expected  motions  of  the  array.  They  are  probably  due  to  path  de- 
pendent temperature  fluctuations,  as  well;  thus  these  navigation  data  are  considered 
largely  unusable. 

We  can  try  to  estimate  the  array  tilt  by  hydrodynamics.  The  length  of  the  of 
the  mooring  is  approximately  60  meters  and  it  has  net  buoyancy  of  3230  Newtons. 
Assuming  laminar  flow  in  calculating  the  form  drag  for  the  cable  and  flotation  sphere, 
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a  barotropic  current  will  cause  upper  hydrophone  displacements  that  are  roughly 
proportional  to  the  square  of  the  current  speed,  giving  1  meter  at  30  cm/s;  2  meters 
at  40  cm/s;  and  3  meters  at  50  cm/s.  The  current  records  obtained  from  a  bottom- 
mounted  acoustic  Doppler  current  meter  (ADCP),  moored  about  3  km  from  the 
WVLA,  showed  barotropic  currents  that  seldom  exceeded  30  cm/s,  with  occasional 
baroclinic  currents  from  large  solitons  of  up  to  50  cm/s. 

The  effects  of  baroclinic  currents  on  mooring  motion  are  hard  to  predict,  but  the 
movement  at  the  top  hydrophone  is  certainly  less  than  what  would  be  produced  by 
a  barotropic  current  of  the  same  magnitude.  This  line  of  reasoning  predicts  mooring 
motions  that  are  generally  less  than  a  meter  in  amplitude  at  the  top  hydrophone 
(equivalent  to  1  degree  of  tilt)  with  occasional  excursions  that  might  approach  2  me- 
ters. A  more  detailed  analysis  of  the  WVLA  navigation  data  is  provided  in  Appendix 
A. 


3.3      Mode  Filtering 

3.3.1      Errors  from  Array  Tilt 

As  mentioned  earlier,  phase  differences  caused  by  array  tilt  will  cause  errors  when 
the  16  channels  of  acoustic  data  from  the  WVLA  are  converted  to  normal  mode 
arrivals  through  mode-filtering.  Figure  3-1  illustrates  the  effects  of  array  tilt  on  some 
typical  400  Hz  mode-filter  outputs;  the  inputs  are  pure  mode  shapes  from  the  same 
sound  speed  profile  (SSP)  as  the  mode-filter.  With  no  tilt,  all  four  mode  filters  do 
an  exceptional  job  of  passing  their  respective  mode  and  rejecting  all  others.  The 
mode-filters  continue  to  perform  acceptably  for  small  deflections  (less  than  1  meter) 
of  the  upper  hydrophone,  but  most  become  significantly  degraded  for  deflections 
approaching  2  meters.  With  the  exception  of  mode  2  rejection,  the  mode  1  filter 
performs  robustly  all  the  way  to  2  meters  of  tilt. 

45 


Typical  Tilt  Dependence 
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Figure  3-1:  Mode  rejection  levels  from  mode-filter  outputs,  Ai(*),  A2(x),  ^3(0),  and 
A4(+),  as  a  function  of  array  tilt.  The  four  plots  represent  tilted  inputs  of  pure  mode 
1,  mode  2,  mode  3,  and  mode  4  arrivals.  These  plots  use  mode  filters  and  mode 
arrivals  generated  from  the  same  typical  SWARM  SSP. 
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Additional  Error  from  Mode  Shape  Differences 

Another  source  of  error  for  the  mode  filter  is  the  difference  between  the  sound  speed 
profile  used  in  creating  the  mode  filter  and  the  actual  sound  speed  profile  during 
the  acoustic  arrival.  The  22  sequences  of  a  transmission  are  split  into  two  groups  of 
eight  and  one  group  of  six  sequences  for  further  analysis.  The  mode  filter  for  each 
group  is  calculated  by  using  the  WVLA  thermistor  record  samples  whose  sample  time 
is  closest  to  the  arrival  time  of  the  middle  sequence  of  the  group.  The  thermistors 
on  the  WVLA  are  sampled  every  30  seconds,  and  an  eight  sequence  group  covers 
about  41  seconds.  In  only  three  of  the  22  sequences  is  the  arrival  of  a  particular 
sequence  offset  from  the  sample  time  of  the  mode-filter  by  more  than  30  seconds,  and 
the  worst  of  these  three  is  42  seconds.  Under  most  circumstances  this  error  will  be 
insignificant,  but  during  the  passage  of  a  soliton,  offsets  of  30  seconds  can  produce 
noticeable  effects. 


3.3.2     Time  Dependence  of  Filter  Performance 

Even  with  zero-tilt  and  no  difference  in  mode  shapes,  modal  cross-talk  cannot  be 
avoided  completely  (see  Figure  3-1  at  zero-tilt).  This  minimum  level  of  cross-talk  is 
borne  of  under  sampling  and  partial  coverage  (i.e.  discretization  error).  The  exact 
amount  of  cross-talk  is  dependent  on  the  mode  shapes  in  question,  and,  as  such,  the 
zero-tilt  filter  performance  at  the  WVLA  is  time  varying.  Figure  3-2  plots  the  "Best 
Case"  (i.e.  zero-tilt  with  no  time  offset)  performance  of  the  mode  1-4  filters  calculated 
from  thermistors  on  the  WVLA  during  the  period  between  about  1300  on  yearday  216 
and  0130  on  yearday  217.  The  temporal  variation  is  significant,  but  the  cross-mode 
rejection  levels  never  drop  below  20  dB.  The  mean  "Best  Case"  mode  filter  rejection 
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Best  Case  Time  History 
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Figure  3-2:  Mode  rejection  levels  from  mode-filter  outputs,  Ai(*),  A2(x),  A3(o),  and 
Aj(-r-),  as  a  function  of  transmission  number.  The  four  plots  are  based  on  a  series 
of  untilted  inputs  of  pure  mode  1,  mode  2,  mode  3,  and  mode  4  arrivals,  where  the 
mode  filters  and  mode  arrivals  are  generated  from  the  same  series  of  SSPs. 
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levels  from  Figure  3-2  are 


^oiogfvM2'-')^2'"')^ 


p(zh) 


0  47  46  40 

49  0  34  31 

46  34  0  26 

40  31  26  1 


(3.10) 


in  matrix  form,  where  levels  are  in  dB. 

The  above  performance  levels,  as  stated,  are  optimistic;  they  may  be  taken  as 
upper  bounds.  The  "Worst  Case"  performance  levels  for  this  same  period,  shown  in 
Figure  3-3,  can  be  thought  of  as  lower  bounds.  They  assume  a  constant  tilt  of  one 
degree  and  a  30  second  offset  between  the  arrivals  and  the  mode  filters.  As  stated 
before,  tilts  may  exceed  1  m  from  time  to  time  and  time  offsets  of  up  to  42  seconds 
are  possible,  but  the  "true  worst  case"  is  a  fairly  elusive  concept;  mode  shapes  are 
only  tabulated  in  30  second  intervals,  and  the  tilt  distribution  tail  function  is  difficult 
to  quantify.  Given  all  this,  the  term  "nearly  worst  case"  might  be  better;  still,  the 
tilt  and  time  offset  parameters  chosen  for  Figure  3-3  will  bound  at  least  90  percent  of 
the  arrivals,  and  parameters  sufficient  to  bound  the  performance  100  percent  of  the 
time  would  provide  little  in  terms  of  a  useful  comparison. 

The  "Worst  Case"  filter  performance  is  significantly  reduced,  relative  to  the  "Best 
Case"  scenario.  Most  of  this  reduction  can  be  attributed  to  the  assumed  one  degree 
of  tilt,  but,  occasionally,  the  passage  of  a  soliton  at  the  WVLA  causes  a  large  rapid 
change  in  mode  shapes  during  a  mode  arrival.  The  assumed  30  second  time  difference 
between  arrivals  and  filters  can  then  cause  an  additional  marked  drop  in  performance. 
The  drop  seen  at  transmission  number  27  is  one  such  occurrence.  Overall,  the  "Worst 
Case"  mode  filter  rejection  levels  from  Figure  3-3  are  still  quite  acceptable,  but  there 
are  indications  that  the  quality  of  individual  arrivals  may  be  relatively  poor  from 
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Worst  Case  Time  History 
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Figure  3-3:  Mode  rejection  levels  from  mode-filter  outputs,  Ai(*),  ^(x),  ^3(0),  and 
A4(+),  as  a  function  of  transmission  number.  The  four  plots  are  based  on  1  degree 
tilted  inputs  of  pure  mode  1,  mode  2,  mode  3,  and  mode  4  arrivals,  using  mode  filters 
and  mode  arrivals  generated  from  the  SSPs  separated  by  30  seconds. 


50 


time  to  time.  The  mean  "Worst  Case"  values  in  matrix  form  are  (again  in  dB  levels) 


2Qlo  (™<}>i{zh,t)<t>*{zh,t-W)&z 


P(zk) 


0  13  24  32 

13  1  11  20 

24  11  1  11 

32  21  11  2 


(3.11) 


Where  (f>*(zh,t  —  30)  is  a  one  degree  tilted  version  of  <f>j(zh,t)  that  is  also  offset  in 
time  by  30  seconds. 

A  reasonable  estimate  of  the  expected  performance  for  the  four  mode  filters  might 
be  found  by  taking  the  average  of  the  mean  of  both  cases.  The  result  is 

0  31  35  36 

31  1  22  26 

34  22  1  18 

36  26  18  1 

in  matrix  form.  These  levels  are  quite  good,  but  they  must  be  considered  in  light  of 
the  potential  for  large  differences  between  modal  intensities  at  the  receiver  (e.g.  the 
average  filtered  peak  height  differential  between  mode  1  and  mode  4  is  about  10  dB.) 


3.3.3      Mode  Filter  Outputs 

As  discussed  previously,  we  do  not  expect  to  achieve  perfect  resolution  of  a  given 
mode  arrival,  but  we  may  reasonably  expect  the  lower  mode-filter  outputs  to  generally 
reflect  the  true  mode  arrival  structures.  Figure  3-4  shows  mode  1-4  filtered  outputs 
for  ten  sequences,  each  separated  by  6  minutes,  over  a  54  minute  period.  Most  of 
the  arrival  structures  do  not  seem  to  follow  any  clearly  recognizable  patterns,  and 
the  temporal  variability  over  a  six  minute  period  is  quite  dramatic.  The  patterns  are 
basically  uncorrected. 
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Figure  3-4:  Mode  Coefficient  outputs,  Ai(solid),  ^(dashed),  ^(dot-dash),  and 
Aj(dotted),  as  a  function  of  time.  The  ten  plots  depict  the  arrival  structure  for 
the  first  of  22  captured  sequences  out  of  ten  sequential  transmissions. 
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Figure  3-5  illustrates  the  temporal  variability  of  mode  arrivals  over  a  shorter  time 
scale  by  following  the  evolution  of  the  time  zero  transmission  in  Figure  3-4  over 
a  period  of  92  seconds.  This  figure  shows  some  fairly  strong  correlations  between 
sequences  separated  by  upwards  of  one  minute;  a  more  quantitative  measure  of  the 
temporal  coherence  can  be  reached  through  computing  the  correlation  coefficient 
between  the  first  sequence  of  a  transmission  and  the  twenty  subsequent  sequences 
in  that  transmission.  Figure  3-6  displays  the  results  of  these  computations  when 
averaged  over  ten  transmissions.  For  the  analyzed  time  period,  the  average  3  dB 
down  points  for  the  correlation  functions  fall  between  70  and  90  seconds.  A  good 
question  to  try  to  answer  is  why  the  lowest  mode  decorrelates  least  quickly  and  the 
highest  mode  decorrelates  most  quickly.  It  is  probably  related  to  the  percentage  of 
scattered  energy  in  each  arrival  pattern,  and  it  should  be  investigated  further  in  the 
future.  However,  it  will  not  be  simple,  as  the  broken  sampling  schedule  at  the  WVLA 
is  not  very  conducive  to  correlation  analysis.  Significant  levels  of  decorrelation  often 
do  not  appear  over  the  course  of  a  transmission. 

The  arrival  structures  we  have  observed  are  clearly  indicative  of  a  large  amount 
of  mode  scattering  between  the  source  and  the  receiver.  The  source  of  most  of  the 
scattering  is  very  probably  the  numerous  internal  wave  solitons  observed  by  the  ther- 
mistors on  the  WVLA  and  other  moorings  deployed  during  SWARM.  A  numerical 
model  for  mode  dispersion  and  coupling  will  be  a  useful  tool  in  identifying  the  relevant 
statistics  of  these  complicated  arrival  patterns. 
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Figure  3-5:  Mode  coefficient  outputs,  Ai(solid),  j42(dashed),  A3 (dot- dash),  and 
Ai(dotted),  as  a  function  of  time.  The  ten  plots  depict  the  arrival  structure  for 
every  other  sequence  from  the  first  to  the  19th  of  22  captured  sequences  from  the 
time  zero  transmission  in  Figure  3-4. 
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Mode  1-4  Average  Correlation  Coefficient 


1 
0.95  - 
0.9 

0.85 


c 
_® 

£    0.8- 
o 

c0. 
o 

JS 

£    0.7 

o 

O 

0.65 
0.6 

0.55 
0.5 


'•.■*.—"■! 

1 1 1 1 1 1 1 

1 

i 

■>».-.,..  ^.  •  -.•'«s. .- 

■s.-^  :       \^  : 

: : : :...;-....  \.s j/s^. 

:               :              :               :               :            '''.:  ^  n 

s 

\ 

•s. 
•  \           —  - 

:• .  n 

*    X 

X 

'■■.">, 

\  '"■-.. 

1              1              1              1              1              1 

1 

10  20  30  40  50  60 

Elapsed  Time  (s) 


70 


80 


90 


100 


Figure  3-6:  Average  correlation  coefficients  for  outputs,  Ai(solid),  ^(dashed), 
A3(  dot- dash),  and  A4(dotted),  as  a  function  of  elapsed  time.  The  average  is  taken 
over  the  ten  transmissions  of  Figure  3-4. 
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Chapter  4 


A  Simple  Scattering  Model 


Several  avenues  are  available  for  developing  a  model  that  exhibits  the  main  features 
of  the  scattering  behavior  seen  in  the  acoustic  data.  We  will  start  with  the  simplest 
possible  time  domain  model  and  examine  some  of  its  behavior;  we  will  then  transition 
to  a  broadband  model  that  includes  added  complexity. 

4.1      Normal  Mode  Pulse  Propagation 

In  order  to  perform  analysis  of  the  acoustic  arrivals  at  the  WVLA,  the  acoustic  signal 
transmitted  from  the  400  Hz  source  is  treated  as  a  sum  of  vertical  modes 

p8{z,t)  =  An<l>n{z)p{t),  (4.1) 

where  repeated  indices  indicate  summation,  <j)n{z)  is  the  mode  shape,  An  is  the  mode 
coefficient,  and  p(t)  is  the  pulse-shape  of  the  transmitted  sequence  after  baseband 
replica  correlation  (see  Figure  2-2).  Each  of  the  modes  transmitted  can  couple  into 
other  modes  when  large  horizontal  sound  speed  gradients  are  encountered  between 
the  source  and  receiver.  The  large  gradients  could  be  caused  by  temperature  fronts, 
salinity  intrusions,  or  the  internal  wave  solitons  which  are  of  particular  interest  here. 
To  estimate  the  coupling  of  modes  across  the  soliton  gradients,  an  approximate 
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single-scatter  one-way  coupled  mode  approach  from  Jensen  et  al.  [13]  is  employed 
that  uses  an  arithmetic  mean  of  pressure  and  velocity  matching  conditions  across  the 
interface  to  produce  the  elements  of  the  coupling  matrix 

where  k3rm  is  the  horizontal  wave  number  of  mode  m  in  region  j,  and  the  coupling 
among  mode  coefficients  is  given  by 

Ai+1  =  CJ^Al,  (4.3) 

or  in  matrix  notation, 

Aj+i  =  c>+1AJ.  (4.4) 

Parabolic  Equation  (PE)  simulations  by  Presig  and  Duda  [2]  indicate  that  the 
mode  coupling  caused  by  a  single  soliton  (whose  width  is  less  than  about  400  m)  can 
be  reasonably  modeled  by  a  sharp  transition  into  and  out  of  the  soliton's  sound  speed 
profile  (SSP).  This  sharp  interface  approximation  (SIA),  illustrated  in  Figure  4-1, 
provides  an  avenue  for  treating  acoustic  transmissions  as  an  ever-increasing  sum  of 
mode  pulses  propagating  through  the  waveguide.  Figure  4-1  also  shows  the  idealized 
SSPs  for  the  regions  inside  and  outside  the  inverted  top  hat.  The  nearly  isovelocity 
background  SSP  is  unrealistic,  but  its  sinusoidal,  frequency-independent,  mode  shapes 
are  useful  in  this  initial,  simplistic,  look  at  the  problem. 

4.1.1      Total  Field  Summation 

If  there  are  J  trapped  modes  in  both  the  background  and  "soliton  filled"  wave  guides, 
each  mode  pulse  will  couple  into  J  further  pulses  as  it  transitions  from  one  "environ- 
ment" to  the  other,  creating  an  J2  increase  in  the  number  of  pulses  for  each  of  the 
double  transitions  seen  in  individual  soliton  interactions.   Extrapolating  this  line  of 
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Figure  4-1:  An  illustration  of  the  SIA  approximation  for  solitons  (left),  and  the 
idealized  SSPs  (right)  used  to  model  the  interior  and  exterior  regions  of  the  SIA 
soliton.  The  dotted  curve  represents  the  exterior  or  background  region,  and  the  solid 
curve  represents  the  interior. 


thought  to  a  number  of  soliton  interactions  yields  the  following  time-domain  pulse 
synthesis  representation  for  the  signal  at  the  WVLA 


Pwvla{Z,t)  —  (pa{z)Lab    Cbc         Ccd         ...LjkAk 

x  exp[»A£(rWa  -  rN)  +  ik%-l{rN  -  rN~l)  +  •  •  •  +  ikQrk(rl  -  r0)] 


(4.5) 


xp(t 


(rwvta_rN)    ^    (rAT        N-l)    _ 


77=r 


), 


where  repeated  indices  indicate  summation,  and  u^  is  the  mode  a  group  velocity  in 
region  N. 

Figure  4-2  illustrates  the  effects  of  propagating  four  modes  through  32  km  of  this 
idealized  wave  guide  with  one  200  meter  width  soliton  at  the  20  km  point.  The  solid 
curves  represent  the  expected  mode  coefficients  or  mode-filter  outputs,  A\(t),  Azit), 
As(t),  and  A^{i).  The  dotted  curves  are  the  individually  scattered  subarrivals.  These 
individual  subarrival  pulses  are  not  detectable  in  actual  acoustic  receptions.  However, 
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Figure  4-2:  Propagation  of  four  modes  through  32  km  of  an  idealized  wave  guide 
with  one  200  meter  soliton  at  the  20  km  point.  The  solid  curves  are  the  absolute 
value  of  coherent  sums  of  the  J2  =  16  dotted  individually  scattered  arrivals.  Starting 
amplitudes  for  mode  1  through  4  were  1.0,  0.8,  0.6,  and  0.4  respectively. 
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as  seen  in  Figure  4-2,  we  can  sometimes  infer,  from  the  coherent  sum,  the  relative 
levels  of  influence  from  some  "general"  paths  through  range- modenumber  space. 

Each  of  the  mode  1-4  arrivals  in  Figure  4-2  shows  influence  from  four  different 
"general"  paths  through  range-modenumber  space.  We  call  the  first  hump  in  the  mode 

1  arrival  envelope  a  pseudo-adiabatic  mode  1  (PAM1)  arrival,  as  its  peak  arrives  on 
or  very  near  the  time  predicted  for  an  adiabatically  coupled  mode  1  arrival.  The 
general  path  for  this  PAM1  arrival  hump  could  be  referred  to  as  "the  1:1  path"  in 
this  simple  case,  because  all  the  energy  traveled  in  mode  1  before  reaching  the  soliton, 
took  various  paths  through  the  soliton,  and  transited  the  remaining  11.8  km  in  mode 
1.  The  energy  in  the  second  hump  of  the  mode  1  arrival  envelope  (2:1  path)  traveled 
the  first  20  km  in  mode  2,  took  various  paths  through  the  soliton,  and  transited  the 
remaining  11.8  km  in  mode  1.    Likewise,  the  energy  in  the  first  hump  of  the  mode 

2  arrival  envelope  (1:2  path)  traveled  the  first  20  km  in  mode  1,  took  various  paths 
through  the  soliton,  and  transited  the  remaining  11.8  km  in  mode  2.  Similar  "general" 
paths  could  be  inferred  from  the  remaining  humps  in  the  arrival  envelopes. 

The  addition  of  a  second  soliton  in  the  path  increases  the  complication  significantly 
(Figure  4-3);  each  mode  arrival  pattern  now  consists  of  the  coherent  sum  of  256 
(=  J2  x  J2,  where  J  =  4)  pulses  that  are  shifted  in  phase,  amplitude,  and  time.  The 
arrival  envelopes  have  more  humps  and  inferences  about  their  general  paths  are  more 
difficult  to  make.  The  first  humps  in  the  mode  1-4  envelopes  might  be  reasonably 
labeled  1:1:1,  1:1:2,  1:1:3,  and  1:1:4  respectively,  and  their  last  humps  4:4:1,  4:4:2, 
4:4:3,  and  4:4:4  respectively,  but  with  the  interior  it  becomes  more  difficult.  One 
could  consider  16  possible  general  paths  for  each  envelope,  but  the  near-symmetry  of 
the  soliton  placement  in  the  waveguide  would  create  very  little  separation  between 
some  of  them  (e.g.  the  1:2:1  and  2:1:1  paths  in  the  mode  1  arrival  envelope). 

Most  of  the  humps  do  not  correspond  to  any  particular  general  path.  Instead,  they 
are  the  result  of  interference  patterns  between  two  or  more  general  paths.  This  inter- 
ference actually  occurs  to  a  greater  or  lesser  extent  in  all  the  humps.  The  4:4:4  path 
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Figure  4-3:  Propagation  of  four  modes  through  32  km  of  an  idealized  wave  guide  with 
200  meter  solitons  at  the  10  and  20  km  points.  The  solid  curves  are  the  coherent  sums 
of  the  TV2  x  JV2  =  256  dotted  individually  scattered  arrivals.  Initial  mode  1  through 
4  amplitudes  were  1.0,  0.8,  0.6,  and  0.4  respectively. 
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(i.e.  the  PAM4  arrival),  for  example,  shows  little  interference,  while  the  interference 
in  the  neighborhood  of  the  1:1:3  path  is  significant. 

The  total  summation  approach  can  provide  some  valuable  insight  when  dealing 
with  one  or  two  solitons,  but  it  is  quickly  swamped  by  the  exponential  growth  in 
subarrival  pulses  for  multiple  solitons  in  the  path  between  the  source  and  receiver 
(e.g.  65,536  pulses  per  mode  with  four  modes  and  four  solitons). 

4.1.2      Step  Summation 

A  computationally  tractable  alternative  to  the  above  method  is  to  sum  the  pulses 
for  each  mode  after  each  coupling  event  and  create  a  new  pulse  shape  that  can  be 
propagated  to  the  next  step.  The  new  pulse  shape  for  step  n  is  given  by 

aw  =  cab  x  exp[i*srv  -  rn_1)]  x pr\t  +  {rn~ZT1]  -  ^'ZVh-  (4.6) 


U->  u 


b 

The  time  shifts  are  again  tracked  as  lag  with  respect  to  the  PAM1  arrival.  These 
time-shifts  are  subject  to  round-off  error,  if  the  step  size  is  too  small. 

The  step  summation  efficiently  handles  numerous  sequential  coupling  events.  Fig- 
ure 4-4  illustrates  temporal  variations  in  mode  propagation  through  7  solitons  as  the 
solitons  progress  in  lock-step  (0.8  m/s)  up  the  waveguide.  Even  though  lock-step 
propagation  is  unrealistic,  it  is  shown  as  a  type  of  lower  bound  on  the  variation 
in  mode  arrivals  caused  by  soliton  propagation.  Actual  solitons  have  phase  speeds 
that  are  inherently  variable,  and  they  can  propagate  in  offset  directions  that  cause 
dramatic  increases  in  apparent  phase  speed.  A  pattern  of  propagation  with  varying 
phase  speeds  causes  even  more  variation  in  the  arrival  patterns,  but  it  is  not  clear  as 
to  how  much  variation  in  phase  speed  would  constitute  a  "more  realistic"  environ- 
ment. Even  with  7  solitons,  the  PAM1  arrival  is  still  a  dominant  feature  in  many  of 
the  realizations. 

These  simulated  arrivals  show  time-spreading  and  distortion  levels  that  are  quite 

63 


Scattering  from  Seven  Solitons 


Time:  0  min 


Time:  36  min 


0.5 


/\ 

jr     V 

£>?«- -*;C"  "■"■'■■ 

T-i  ■         *»■     .    -  .. 

Time:  48  min 


0.5 


/ 

^vsT"-...-, 

50  100 

Time  (ms) 


0.5 


Time:  6  min 


Time:  54  min 


* "  \  — 

A    x-^ 

•&S---J   'lv 

0  50  100 

Time  (ms) 


Figure  4-4:  Simulated  mode  coefficient  outputs,  A\ (solid),  ^(dashed),  A3(dot-dash), 
and  A4 (dotted),  for  propagation  through  seven  idealized  solitons.  The  solitons  prop- 
agated toward  the  receiver  in  lock  step  at  0.8  m/s.  The  ten  plots  depict  the  arrival 
structures  for  sequential  realizations  at  six  minute  intervals. 
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Figure  4-5:  Time-series  of  temperature  at  a  depth  of  12.5  meters  on  the  WVLA. 


similar  to  the  mode-filtered  outputs  of  Figure  3-4. 


4.2      Incorporation  of  Field  Measured  SSPs 

The  WVLA  was  instrumented  with  six  thermistors  that  recorded  temperature  samples 
every  30  seconds.  This  record  is  principally  used  to  create  time-dependent  mode  filters 
that  are  used  to  analyze  the  acoustic  arrivals  at  the  array,  but  it  is  also  a  convenient 
source  of  realistic  SSPs  for  use  in  the  pulse  propagation  model. 

Solitons  influence  the  SSP  through  almost  the  entire  water  column,  but  the  largest 
temperature  fluctuation  signals  occur  at  depths  within  the  main  thermocline.  Figure 
4-5  shows  50  minutes  of  temperature  data  from  the  shallowest  (12.5  m)  thermistor 
on  the  array.  The  main  thermocline  for  SWARM  usually  has  a  top  edge  at  about 
8  to  10  m  depth,  so  the  location  of  this  thermistor  is  well  suited  for  the  detection 
of  solitons;  the  SSPs  from  minute  777  and  minute  781  were  chosen  as  representative 
of  profiles  with  and  without  a  soliton.  These  SSPs,  including  extrapolations  above 
12.5  m  and  below  60.5  m,  are  shown  in  Figure  4-6.   The  iso-velocity  model  utilized 
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Figure  4-6:  Field  measured  SSPs  for  use  in  the  pulse-propagation  model.  The  dashed 
line  is  representative  of  a  standard  SSP  with  no  soliton.  The  solid  line  represents  the 
soliton  SSP. 


a  pressure  release  bottom  for  simplicity.  The  field  measured  SSPs  were  coupled  to  a 
bottom  that  is  based  on  a  sedimentary  analysis  in  the  vicinity  of  the  WVLA  [1].  The 
bottom  has  a  sound  speed  that  increases  linearly  from  1500  to  1800  m/s  in  a  30  m, 
1.8  g/cm3  density,  sub  layer  over  a  rigid  bottom.  The  best  estimate  of  compressional 
wave  attenuation  in  this  sublayer  is  0.092  dB/m  at  400  Hz  (Altan  Turgut,  NRL, 
Washington,  D.C.,  private  communication),  but  as  later  results  will  indicate,  this 
number  is  either  too  large  by  at  least  a  factor  of  four  or,  more  likely,  is  not  applicable 
for  the  majority  of  the  actual  SWARM  waveguide. 

These  SSPs  should  provide  the  model  with  mode  dispersion  characteristics  that 
are  closer  to  reality  than  the  iso-velocity  SSP  used  previously,  but  they  also  bring  to 
light  another  question  about  the  narrowband  approximation  of  the  model. 
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4.2.1      Validity  of  400  Hz  Narrowband  Models  for  a  Broad- 
band Pulse 

Working  with  an  iso-velocity  water  profile  and  a  pressure  release  bottom  provides 
several  advantages,  including  the  fact  that  the  sinusoidal  mode  shapes  are  frequency 
independent.  As  a  result,  the  mode  group  speeds  are  also  frequency  independent. 
The  simple  pulse  propagation  model  previously  discussed  used  a  narrowband  approx- 
imation in  the  determination  of  coupling  matrices  and  group  velocities.  The  validity 
of  this  approximation  becomes  a  significant  issue  with  the  introduction  of  realistic 
SSPs  to  the  model. 

Figure  4-7  compares  soliton  present  and  soliton  absent  mode  shapes  at  the  ex- 
tremes of  the  pulse  frequency  domain  (355-445  Hz).  The  mode  shapes  for  these 
SWARM  SSPs  have  a  frequency  dependence  that  appears  minimal  in  modes  1  and  2, 
but  it  becomes  somewhat  more  visible  in  modes  3  and  4.  The  mode  shape  differences 
for  the  higher  modes  will  contribute  to  some  frequency  dependence  in  the  coupling 
matrices,  particularly  coupling  in  and  out  of  mode  4.  This  is  not  considered  in  the 
narrowband  pulse  propagation  model;  the  pulse  is  assumed  to  travel  and  couple  in 
accordance  with  the  400  Hz  mode  shapes. 

Mode  shape  frequency  dependence  also  leads  to  frequency  dependence  in  group 
velocity.  The  group  velocity  is  defined  as 

CLKrn 

Jensen  et  al.  [13],  using  perturbation  theory,  derive  an  equation  for  the  reciprocal  of 
the  group  velocity, 


au        kT  Jo 


*2{Z)   -dz,  (4.8) 


p{z)c2{z) 

that  illustrates  the  dependence  of  group  velocity  on  mode  shape.  Figure  4-8  shows 
the  degree  of  frequency  dependence  in  phase  and  group  speed  for  the  first  four  modes 
of  the  standard  field  measured  SSP.  Pulse  distortion  from  frequency  dispersion  in  the 
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Figure  4-7:  Mode  shapes  for  field  measured  (SWARM)  SSPs.  The  dashed  lines  are 
the  355  Hz  (thin  lines)  and  445  Hz  (thick  lines)  mode  shapes  corresponding  to  the 
standard  SSP  with  no  soliton.  The  solid  lines  are  the  same  modes  for  the  soliton  SSP. 
The  dotted  line  at  68  m  indicates  the  sediment  interface. 
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Figure  4-8:  Phase  and  group  speeds  for  the  standard  field  measured  SSP.  Mode 
1  has  the  lowest  phase  speed  and  the  highest  group  speed.  The  remaining  phase 
speeds  increase  monotonically  with  mode  number,  and  the  group  speeds  decrease 
monotonically  with  mode  number. 
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higher  modes  will  clearly  occur;  the  present  time-domain  model  does  not  replicate 
this. 

The  switch  from  a  time-domain  propagation  model  to  the  frequency  domain  is 
fairly  straight  forward.  Instead  of  simply  propagating  the  replica  correlated  pulse 
shape,  p(t),  at  the  400  Hz  modal  phase  and  group  speeds,  allow  its  individual  Fourier 
components  P{oj),  to  propagate  prior  to  synthesis,  the  resulting  transmitted  pulse 
shape, 

Pn(r,  t)  =  ~  J  PCwJe^e-*'^  (4.9) 

will  show  the  spreading  and  intensity  loss  that  is  due  to  frequency  dispersion.  The 
synthesis  results  for  different  modes  will  reflect  the  differences  in  group  speed  between 
modes  through  the  frequency  dependence  of  the  horizontal  wavenumbers,  krn(u>). 
The  narrowband  SIA  coupling  matrices  used  in  the  time  domain  model  are  linear 
operators  that  apply  equally  in  the  frequency  domain.  Figure  4-9  illustrates  both 
the  separating  effects  of  mode  dispersion  and  the  pulse  distortion  effects  of  frequency 
dispersion  that  might  be  expected  for  individual  mode  pulses  that  travel  the  entire  32 
km  waveguide  without  any  soliton  scattering.  An  ideal  frequency  domain  propagation 
model  would  also  include  frequency  dependence  in  the  coupling  matrices,  and  in  fairly 
simple  waveguides,  this  is  a  possibility  (e.g.  the  KRAKEN  [12]  field  synthesis  method, 
described  in  Appendix  B),  but,  as  also  shown  in  Appendix  B,  the  effects  are  minimal 
and  the  added  complexity  and  computation  time  is  not  warranted. 

4.2.2      Adding  Modal  Attenuation 

As  mentioned  earlier,  one  of  the  estimated  properties  for  the  field  measured  SSPs  was  a 
compressional  wave  attenuation  of  0.092  dB/m  in  the  bottom.  This  causes  significant 
attenuation  of  the  higher  modes  that  interact  with  the  bottom.  The  attenuation  can 
be  modeled  through  an  imaginary  term,  ia,  added  to  the  horizontal  wave  number,  kT. 
The  KRAKEN  [12]  calculated  400  Hz  mode  attenuation  levels  for  the  SWARM  SSPs 
are  given  in  Table  4-1.    The  presence  of  the  soliton  depresses  the  SSP  significantly, 
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Figure  4-9:  Effects  of  mode  and  frequency  dispersion  on  single  mode  transmissions 
over  32  km.  The  transmissions  include  no  soliton  scattering  or  intrinsic  mode  atten- 
uation. The  attenuation  in  the  higher  modes  is  due  solely  to  frequency  dispersion. 


Mode 

Attenuation  (km  *) 

Standard 

Soliton 

1 

0.0094 

0.0235 

2 

0.0706 

0.1782 

3 

0.2042 

0.3323 

4 

0.3060 

0.4090 

5 

0.3751 

0.5204 

6 

0.4639 

0.5547 

7 

0.5465 

0.5742 

8 

0.7132 

0.7457 

Table  4.1:  Intrinsic  mode  attenuation  parameters,  7,  for  various  mode  numbers. 
Values  for  both  a  soliton  SSP  and  the  standard  field  measured  SSP  without  a  soliton 
present  are  given. 
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Figure  4-10:  Effects  of  mode  and  frequency  dispersion  combined  with  intrinsic  mode 
attenuation  on  single  mode  transmissions  over  10  km.  The  transmissions  include  no 
soliton  scattering. 


causing  more  bottom  interaction  and  greater  modal  attenuation.  Even  the  standard 
SSP  attenuation  levels  are  quite  dramatic;  figure  4-10  illustrates  the  effects  of  their 
application  to  a  10  km  modal  transmission.  A  field  with  an  average  number  of 
solitons  would  have  higher  attenuation  levels  applied  over  about  10  percent  of  the 
range,  causing  even  greater  attenuation.  These  attenuation  levels  are  probably  too 
high,  considering  observed  levels  of  time-spreading  in  actual  data.  If  left  at  this  level, 
the  attenuation  will  all  but  eliminate  modes  3  and  4  after  a  full  32  km  transmission. 
As  will  be  shown  in  Chapter  5,  an  "effective  range  averaged"  compressional  wave 
attenuation  of  0.0115  dB/m  is  more  appropriate. 
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4.2.3      The  Effects  of  Bottom  Topography 

Waveguides  with  little  horizontal  variation  in  water  column  properties  may  still  ex- 
perience mode  coupling,  if  there  is  significant  variability  in  bottom  topography.  For 
soliton  filled  waveguides,  the  coupling  caused  by  bottom  topography  produces  a  base- 
line level  of  scattering  upon  which  is  added  time-varying  volume  (e.g.  soliton)  scat- 
tering. Given  the  bottom  topography  of  the  SWARM  waveguide  (see  Figure  2-1), 
some  noticeable  level  of  baseline  scattering  is  expected. 

The  simulations  shown  in  Figure  4-11,  which  used  a  range  dependent  waveguide 
based  on  the  bottom  topography  of  Figure  2-1,  can  provide  a  good  estimate  of  this 
baseline.  The  range  independent  SSP  used  in  the  "No  Ducting"  case  is  a  modified 
version  the  "non-soliton"  SSP  of  Figure  4-6.  The  modification  consists  of  imposing 
the  50.5  m  temperature  from  50.5  m  to  the  bottom,  producing  a  "worst  case"  sce- 
nario for  mode  1  scattering.  The  "Ducting"  SSP  is  the  same,  except  the  60.5  m 
temperature  is  imposed  from  50.5  m  to  the  bottom  depth,  noticeably  reducing  the 
topographic  coupling  effects  for  the  lower  modes.  The  thinner  lines  in  this  figure 
represent  predicted  arrival  shapes  for  modes  that  propagate  adiabatically  (i.e.  no 
scattering  of  energy  between  modes). 

4.3      Incorporation  of  a  Background  Internal  Wave 
Spectrum 

Although  the  primary  modeling  concern  in  this  paper  is  internal  wave  solitons,  linear 
internal  wave  fields,  which  are  also  present  in  shallow  water,  should  not  be  neglected. 
They  are  usually  modeled  and  studied  as  a  stochastic  process  that  produces  random 
sound-speed  fluctuations,  6c(r,  z,t),  about  a  background  profile,  cq{z).  A  "standard 
method"  of  producing  these  random  sound  speed  fluctuations  and  resultant  coupling 
matrices  in  deep  water  was  developed  by  Dozier  and  Tappert  [14];  it  generally  con- 
sists of  using  a  GM  spectrum  as  the  basis  for  creating  a  random  realization  of  the 
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Figure  4-11:  Mode  coupling  effects  of  the  bottom  topography  shown  in  Figure  2-1 
when  the  SSP  is  a  modified  version  the  "non-soliton"  SSP  of  Figure  4-6.  The  "No 
Ducting"  modification  consists  of  imposing  the  50.5  m  temperature  from  50.5  m  to 
the  bottom  depth,  and  the  "Ducting"  modification  consists  of  imposing  the  60.5 
m  temperature  from  50.5  m  to  the  bottom  depth.  The  thinner  lines  in  this  figure 
represent  predicted  arrival  shapes  for  modes  that  propagate  adiabatically. 
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background  internal  wavefield,  £(r,  z,t),  which  is  proportional  to  8c(r,  z,t).  The  re- 
sulting range  and  depth  averaged  <  (Sc/co)2  >  value  of  the  field  controls  the  rate 
at  which  energy  is  transferred  between  modes.  The  amount  of  energy  transferred 
between  modes  is  estimated  by  sampling  6c(r,  z,  i)  and  computing  coupling  matrices 
at  a  fairly  coarse  interval  (Dozier  and  Tappert  [15]  use  1  km).  The  results  are  then 
averaged  over  numerous  realizations. 

It  is  unclear  as  to  how  close  a  GM  spectrum  approximates  shallow  water  internal 
waves.  Lynch  et  al.  [16]  present  a  spectral  analysis  of  data  from  a  thermistor  string 
in  400  m  of  water  that  is  consistent  with  a  GM  spectrum  at  frequencies  above  about  8 
cpd,  but  the  spectrum  in  the  3  to  8  cpd  region  is  significantly  less  energetic  relative  to 
GM.  This  single  string  could  not  provide  any  information  concerning  the  directional 
spectrum  in  the  region. 

In  the  absence  of  any  other  data  to  the  contrary,  the  GM  spectrum  assumption 
provides  a  method  of  developing  somewhat  realistic  levels  of  vertical  and  horizontal 
coherence.  The  vertical  coherence  is  critical  to  the  calculation  of  the  coupling  ma- 
trix at  a  given  range  and  the  horizontal  coherence  provides  some  realism  to  the  range 
dependence  of  these  matrices.  On  the  other  hand,  since  the  GM  spectrum  may  not  re- 
ally be  representative  of  shallow  water,  the  modeling  of  correlation  functions  through 
a  rather  complicated  GM  internal  wave  synthesis  may  not  be  the  best  approach. 

4.3.1      EOF  Soundspeed  Perturbation  Modes 

Empirical  orthogonal  function  (EOF)  modes  of  thermistor  string  temperature  pertur- 
bations are  a  good  alternative  to  internal  wave  synthesis.  Unlike  the  GM  synthesis 
approximations,  the  EOF  modes  will  include  influence  from  both  linear  and  nonlinear 
internal  waves,  so  the  EOF  modes  can  be  used  as  stand  alone  random  realizations. 
The  range  and  depth  averaged  <  (6c/ cq)2  >  value  of  an  EOF  realization  is  eas- 
ily scaled  up  or  down.  When  scaled  down  it  can  be  used  in  conjunction  with  SIA 
solitons  as  a  background  field,  providing  additional  scattering  between  SIA  solitons. 
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The  main  limitation  associated  with  this  or  any  other  model  that  includes  Dozier 
and  Tappert  [14]  coupling  is  that  the  results  are  only  valid  in  a  stochastic  sense. 
Dozier  and  Tappert  [15]  restrict  their  own  results  to  averages  over  a  hundred  random 
realizations.  The  Dozier  and  Tappert  [14]  computations  and  several  issues  concerning 
the  implementation  of  EOF  fields  are  discussed  in  more  detail  in  Appendix  C. 


4.4     Fine  Scale  Coupled  Mode  Computations 

The  large  step  size  employed  with  Dozier  and  Tappert  [14]  coupling  is  the  primary 
impediment  to  achieving  a  deterministic  mode  arrival  model,  but  as  shown  in  Ap- 
pendix C,  even  when  implemented  with  fine  scale  sampling,  they  are  still  not  reliable 
predictors  of  deterministic  mode  arrival  patterns. 

A  more  accurate  and  workable  fine  scale,  range  dependent,  frequency  domain, 
coupled  mode  model  can  be  developed  using  the  equation  (4.2)  narrowband  coupling 
matrix.  The  range  dependence  of  the  model  is  determined  by  "propagating",  with- 
out dispersion,  the  SSPs  measured  at  the  WVLA  towards  the  source  at  0.8  m/s. 
Sampling  the  WVLA  SSPs  every  minute  yields  a  48  m  step  size  for  the  "propagated 
thermistor  string"  waveguide;  a  series  of  coupling  matrices  can  then  be  calculated  for 
the  transitions  between  the  SSPs. 

Equation  (4.2),  as  it  stands,  works  well  enough  for  the  limited  numbers  of  coupling 
events  required  with  SIA  solitons,  but  it  requires  modification  if  this  finer  scale  (48 
m)  range  dependence  is  to  be  employed.  Ideally  a  numerical  computation  of  equation 
(4.2)  using  the  same  SSP  on  both  sides  would  yield, 

an  identity  matrix,  but  the  six  decimal  place,  1001  point,  ASCII  mode  shapes  provided 
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by  the  KRAKEN  [12]  plotting  routine  "plotmodeM"  yields 

(j>n{z)<f>n{z) 


I 


M 


^1.  (4.11) 


The  integrals  are  very  close  to  one,  typically  around  0.9996  for  the  higher  modes,  but 
the  cumulative  effect  over  hundreds  of  coupling  events  is  noticeable.  The  numerical 
implementation  of  equation  (4.2)  also  produces  small  (4th  decimal  place)  erroneous 
off  diagonal  terms. 

All  these  small  error  sources  can  be  readily  eliminated  for  the  trivial  case  of 
coupling  into  an  identical  environment  by  multiplying  the  equation  (4.2)  coupling 
matrix  by  its  inverse,  creating, 

C'  =  C~1C  =  I,  (4.12) 

the  desired  identity  matrix.  In  a  fine  scale  implementation,  the  coupling  matrices 
are  non-trivial,  but  the  differences  between  SSPs  in  adjacent  regions  are  small.  This 
suggests  a  corrected  coupling  matrix, 

C[2  =  C^CU.  (4.13) 

Kuperman  et  al.  [17]  deal  with  a  truncated  coupling  coefficient  integral  in  a  similar 
fashion.  The  correction  matrix,  C^1,  is  always  close  to  an  identity  matrix,  insignifi- 
cant over  a  small  number  coupling  calculations,  but  it  can  be  important  when  small 
range  steps  are  used  over  large  distances.  Consider  a  32  km  waveguide  with  a  48  m 
step  size.  The  670  coupling  matrices  it  requires  will  generate  significant  cumulative 
errors  from  small  individual  errors  in  coupling  coefficients  (e.g.  .9996670  =  .765). 

The  effectiveness  of  this  correction  factor  can  be  seen  in  the  close  fit  of  the  mode 
arrival  patterns  in  Figure  4-12.  Both  sets  of  patterns  were  generated  using  the  same 
32  km  waveguide,  but  one  used  a  48  m  step  size  (dotted),  and  the  other  employed 
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Figure  4-12:  Comparison  of  fine  scale  coupled  mode  arrivals  using  24  m  (solid)  and 
48  m  (dotted)  step  sizes.  The  32  km  wave  guide  under  evaluation  was  generated  from 
range  propagated  WVLA  thermistor  records. 
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a  24  m  step  size  (solid),  requiring  twice  as  many  coupling  calculations.  The  32  km 
waveguide  is  derived  directly  from  range  propagated  WVLA  thermistor  records.  The 
transformation  from  time  dependence  to  range  dependence  is  made  by  assuming  all 
the  temperature  disturbances  are  propagating  without  dispersion  at  0.8  m/s  along  the 
same  line  as  the  acoustic  transmission  path.  The  reasoning  behind  this  assumption 
is  discussed  in  the  horizontal  EOF  analysis  section  of  Appendix  C.  Had  the  coupling 
matrices  not  been  corrected  properly,  the  smaller  step  size  model,  with  1340  coupling 
events,  would  have  diverged  noticeably  from  the  48  m  step  size  model.  The  closeness 
of  fit  also  indicates  that  the  model  is  at  or  near  solution  convergence  in  terms  of  step 
size  at  the  48  m  point. 

4.5     Putting  it  All  Together 

Starting  with  a  simple  near-isovelocity  narrowband  SIA  soliton  coupling  model,  we 
have  explored  and  considered  various  improvements  on  the  basic  model.  We  will 
now  sharpen  our  focus  and  limit  further  analysis  to  two  basic  model  types.  The  first 
model,  emphasizing  speed,  has  the  following  steps: 

1.  Create  waveguides  by  alternating  between  two  different  SSPs,  making  solitons 
in  the  Preisig  and  Duda  [2]  SIA  fashion.  The  SSPs,  shown  in  Figure  4-6,  are 
from  the  WVLA  at  periods  with  and  without  a  soliton  present. 

2.  Conduct  mode  propagation  in  the  frequency  domain,  using  interpolated,  KRAKEN 
calculated,  modal  wavenumber  and  attenuation  spectra  over  the  bandwidth  of 
the  SWARM  400  Hz  signal. 


3.  Use  the  narrowband  approximation  in  coupling  matrices. 


This  model  will  be  used  primarily  for  the  purpose  of  analyzing  the  effects  of  solitons 
on  mode  1  arrivals.  The  addition  of  bottom  topography  is  feasible,  but  complicated, 
and,  as  shown  in  Figures  4-11  and  4-12,  it  would  have  a  very  limited  impact  on  mode 
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1  scattering.  As  discussed  above,  consideration  was  given  to  adding  background 
scattering,  in  the  form  of  EOF  sound  speed  perturbations  and  Dozier  and  Tappert 
[14]  coupling,  to  this  SIA  model.  This  would  slow  the  model  considerably  and  add 
an  extra  analysis  dimension  in  the  form  of  the  background  field  strength.  This  added 
complexity  will  be  rejected  in  favor  of  an  even  slower,  but  more  straightforward, 
fine  scale  coupling  method  using  range-propagated  WVLA  thermistor  records.  This 
second  model  will: 

1.  Use  32  km  waveguides  derived  directly  from  range  propagated  WVLA  thermis- 
tor records.  The  transformation  from  time  dependence  to  range  dependence  is 
made  by  assuming  all  the  temperature  disturbances  are  propagating  without 
dispersion  at  0.8  m/s  along  the  same  line  as  the  acoustic  transmission  path.  We 
will  use  a  48  m  step  size  that  corresponds  conveniently  to  one  minute  of  the 
thermistor  record. 

2.  Conduct  mode  propagation  in  the  frequency  domain,  using  interpolated  modal 
wavenumber  and  attenuation  spectra  over  the  bandwidth  of  the  SWARM  400 
Hz  signal. 

3.  Use  the  narrowband  approximation  in  coupling  matrices. 

This  method  implicitly  includes  the  soundspeed  perturbations  from  both  linear  and 
non-linear  internal  waves.  The  addition  of  bottom  topography  in  this  model  is  not 
feasible,  because  the  modeled  temporal  variations  in  the  waveguide  rely  on  a  range 
independent  depth. 
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Chapter  5 


Model  Analysis 


5.1      Monte  Carlo  Analysis  of  the  SI  A  Model 

The  mode  arrival  structures  of  actual  data  (Figures  3-4  and  3-5)  are  very  complicated, 
so  choosing  relevant  statistics  to  describe  the  shape  or  distortion  of  arrivals  is  not 
straightforward.  The  SIA  numerical  model  described  in  the  previous  chapter  is  used 
to  analyze  the  sensitivity  of  various  measures  of  mode  arrival  distortion  (e.g.  spread 
and  bias)  to  changes  in  the  soliton  field  between  the  source  and  receiver. 

5.1.1      Relevant  Statistics 

The  main  purpose  of  these  simulation  runs  is  to  study  the  effects  of  varying  the 
numbers  and  locations  of  solitons  on  the  mode  1  arrival  spread  and  bias.  The  spread 
is  a  measure  of  the  width  of  the  arrival  pulse  and  the  bias  is  a  measure  of  the  difference 
between  the  average  arrival  time  of  mode  1  energy  and  the  PAM1  arrival  time. 

Measurement  of  spread  and  bias  in  an  arrival  is  somewhat  subjective;  there  are 
many  similar  mathematical  measures  that  give  slightly  different  results;  the  sub- 
jectivity comes  from  choosing  which  is  "best".  Some  reasonable  bias  measurement 
candidates  are  the  displacements  of  the  peak,  mean,  and  median  arrival  times  of  the 
mode  1  arrival  envelopes  relative  to  the  PAM1  arrival  time.  There  are  several  can- 
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Figure  5-1:  The  mean  (-),  median  (•••),  standard  deviation  (-•),  and  interquartile 
range  (-  -)  of  an  unscattered  mode  1  arrival  plotted  as  a  function  of  noise  cutoff  level. 
The  median  reaches  its  correct  value  of  zero  at  a  noise  cutoff  level  of  -38  dB  relative 
to  the  peak  of  the  arrival. 


didates  for  gauging  the  spread:  the  standard  deviation  of  the  arrival  envelope;  the 
interquartile  range  (the  difference  between  the  the  75th  and  25th  percentile  of  the 
arrival  envelope);  and  the  difference  between  the  mean  and  peak  arrival  time. 

Unfortunately,  several  of  these  commonly  used  statistical  quantities  can  be  cor- 
rupted by  the  trailing  noisy  side  lobes  of  the  replica  correlation  pulse  (Figure  2-2).  A 
clear  picture  of  this  corruption  can  be  seen  by  computing  the  mean,  median,  standard 
deviation,  and  interquartile  range  (IQR)  of  the  pulse  shape  in  figure  2-2  as  a  function 
of  noise  cutoff  level  (Figure  5-1).  At  a  noise  cutoff  of  level  of  40  dB  below  the  peak, 
the  sidelobe  induced  corruptions  in  the  mean  and  standard  deviation  are  about  4 
and  13  ms  respectively.  The  differences  between  initial  (-40  dB)  level  and  asymptotic 
levels  for  the  median  and  IQR  are  quite  small  in  comparison.  The  median  reaches 
its  correct  value  of  zero  at  a  noise  cutoff  level  of  -38  dB,  and  the  IQR  reaches  its 
asymptotic  value  at  a  cutoff  of  -32  dB.  As  the  cutoff  threshold  is  raised,  all  the  errors 
go  to  zero,  with  the  asymptotic  levels  in  the  standard  deviation  and  IQR  consistent 
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Figure  5-2:  The  solid  line  shows  the  mode  number  space  history  of  one  of  the  620 
subarrivals  that  would  make  up  the  mode  1  arrival  from  a  six  mode  transmission 
through  the  10  SIA  soliton  waveguide  that  is  indicated  by  the  dotted  lines. 

with  a  gaussian  pulse  (IQR  is  about  1.35  times  the  standard  deviation  for  a  gaussian 
pulse).  A  noise  cutoff  at  -25  dB  re  peak  was  chosen  as  appropriate  whenever  the 
mean  and  standard  deviations  are  to  be  used  for  spread  and  bias  quantification;  for 
consistency,  it  was  used  in  calculations  of  the  median  and  IQR  as  well. 

The  travel  time  statistics  analyzed  in  this  paper  will,  for  the  most  part,  be  limited 
to  the  above  described  statistics,  but  this  list  is  by  no  means  exhaustive.  Travel 
time  statistics  are  not  even  limited  to  computations  involving  the  time-domain  mode 
arrival  envelope.  Colosi  and  Flatte  [10],  for  example,  estimate  modal  travel  time  by 
determining  the  average  slope  of  modal  phase  as  a  function  of  frequency. 

To  some  degree,  the  spread  and  bias  of  a  mode  1  arrival  depends  on  the  wavenum- 
ber  space  (mode  vs.  range)  history  of  each  subarrival  that  makes  up  a  given  arrival. 
The  solid  line  in  Figure  5-2  shows  the  mode  number  space  history  of  one  of  the  620 
subarrivals  that  would  make  up  a  mode  one  arrival  in  the  10  SIA  soliton  waveguide 
that  is  indicated  by  the  dotted  lines  in  the  same  figure.  The  phase  relationships 
between  the  subarrivals  are  very  important,  but,  on  average,  higher  percentages  of 
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higher  modes  in  the  histories  cause  greater  bias,  and  greater  diversity  in  mode  number 
histories  of  the  subarrivals  cause  greater  spread.  In  other  words,  on  average,  the  bias 
in  an  arrival  is  determined  by  the  mean  modenumber  history  of  the  subarrivals  and 
the  spread  of  the  arrival  is  determined  by  the  variance  of  the  modenumber  histories 
of  the  subarrivals. 

5.1.2      Simulation  Parameters 

Eight  modes  were  propagated  and  scattered  through  the  simulated  waveguides  that 
will  be  discussed  in  this  section.  The  modeled  waveguides  support  upwards  of  40 
propagating  modes  for  the  400  Hz  source,  but  the  modes  above  8  are  of  little  conse- 
quence after  about  5  km.  The  most  serious  drawback  of  this  truncated  sum  approach 
is  exposed  when  one  considers  the  problem  of  adequately  modeling  coupling  events 
that  occur  within  the  first  few  kilometers  of  the  source.  The  "near-source"  coupling 
of  energy  from  modes  above  8  might  reasonably  affect  the  mode  1-8  arrivals  at  the 
receiver.  However,  trials  using  6,  8,  and  10  modes  in  the  model  show  convergence  for 
the  mode  1-4  arrivals  at  the  receiver  when  8  modes  are  used  in  "near-source"  coupling 
scenarios.  The  simple  nature  of  the  SIA  model,  as  implemented,  limits  the  number 
of  other  parameters  open  for  adjustment.  The  "soliton"  and  "no  soliton"  SSPs  are 
fixed  by  design;  other  parameters,  like  source  depth  and  receiver  range,  are  set  to  try 
to  match  the  conditions  in  the  actual  experiment. 

Another  parameter  we  would  like  to  match  to  the  actual  conditions  at  the  SWARM 
site  is  bottom  attenuation.  Chirp  sonar  inversions  for  bottom  parameters  near  the 
WVLA  infer  a  sediment  layer  (specific  gravity  -  1.8)  about  30  m  deep,  whose  sound 
speed  ranges  from  1500  to  1800  m/s  [1].  Given  these  characteristics,  Altan  Turgut 
(NRL,  Washington  D.C.,  private  communication)  estimates  the  compressional  wave 
attenuation  in  the  bottom  to  be  in  the  neighborhood  of  .092  dB/m  at  400  Hz.  Al- 
though this  level  is  on  the  low  end  of  published  attenuation  level  estimates  for  sed- 
iments with  soundspeeds  in  the  1500  to  1800  m/s  range  (e.g.    clay,  silt,  sand,  and 
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gravel)  [13],  the  resulting  modal  attenuation  levels  are  still  quite  significant  (see  Figure 
4-10). 

Figure  5-3  shows  six  random  realizations  of  the  type  of  SIA  soliton  waveguides  that 
will  be  referred  to  in  the  analysis  that  follows.  Using  a  compressional  wave  attenuation 
of  .092  dB/m  in  the  modeled  bottom  of  these  waveguides  yields  predictions  for  mode  1 
spread  that  are  significantly  below  the  values  observed  in  actual  arrivals  (not  shown), 
and  while  a  compressional  wave  attenuation  of  .046  dB/m  (half  of  the  predicted 
value)  provides  simulation  results  that  are  somewhat  better  (Figure  5-4),  they  are 
still  off  the  mark.  A  comparison  of  the  simulated  arrivals  in  Figure  5-4  with  actual 
arrivals  in  Figure  3-4  illustrates  this  point.  The  spread  is  still  not  large  enough, 
and  the  PAM1  arrival  is  much  stronger  than  the  average  real  data  arrival.  Reducing 
the  compressional  wave  attenuation  in  the  bottom  to  0.023  dB/m  provides  even  more 
spread  and  dampens  the  PAM1  arrival  (not  shown),  and  finally,  as  shown  in  Figure  5- 
5,  using  a  bottom  attenuation  of  0.0115  dB/m  produces  arrivals  that  are  very  similar 
to  the  outputs  shown  in  Figures  3-4  and  3-5.  Using  these  comparisons,  an  attenuation 
level  of  0.0115  dB/m  appears  to  be  the  best  choice. 

The  differences  between  results  using  0.023  dB/m  and  0.0115  dB/m  are  not  large; 
limited  Monte  Carlo  simulations  using  both  0.023  dB/m  and  0.0115  dB/m  were  con- 
ducted, and  the  results  were  similar,  but  the  average  mode  1  arrival  bias  levels  were 
about  20  percent  with  0.023  dB/m.  Part  of  the  "required  reduction"  in  the  level  of 
bottom  attenuation  stems  from  the  range  dependence  of  bottom  characteristics;  the 
sediment  layer  is  reduced  to  less  than  1  m  in  thickness  over  much  of  the  path  [1].  A 
reduction  in  the  thickness  of  our  simulated  bottom  to  1  m  reduces  modal  attenuation 
levels  by  a  factor  of  three.  However,  no  attempt  has  been  made  to  introduce  range 
dependence  into  the  bottom  characteristics;  the  simulations  shown  in  the  next  sec- 
tion were  all  calculated  using  the  30  m  sediment  with  an  "effective  range  averaged" 
compressional  wave  attenuation  level  of  0.0115  dB/m. 
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Random  SIA  Waveguides 
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Figure  5-3:  Random  realizations  of  test  waveguides  with  from  10  to  35  200  m  SIA 
solitons  over  32  km.  The  first  soliton  is  placed  randomly  about  2.5  km  from  the 
source,  and  the  remainder  are  added  randomly  from  there,  with  an  average  of  450  m 
between  leading  edges. 
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Figure  5-4:  Single  realizations  of  mode  1  (solid),  mode  2  (dashed),  mode  3  (dot- 
dash),  and  mode  4  (dotted)  arrivals  for  0  to  35  200  m  SI  A  solitons  over  32  km.  The 
simulated  waveguides  used  are  similar  to  those  shown  in  Figure  5-3.  The  simulated 
source  depth  is  27  m  in  68  m  of  water,  and  the  compressional  wave  attenuation  in  the 
bottom  is  .046  dB/m  at  400  Hz.  Note:  Mode  4  arrivals  are  not  visible  at  this  level  of 
attenuation. 
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Figure  5-5:  The  same  scenario  as  Figure  5-4,  except  the  compressional  wave  attenu- 
ation in  the  bottom  has  been  reduced  from  .046  to  0.0115  dB/m. 
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Figure  5-6:  Histograms  for  peak  position  of  mode  1  arrivals  for  5  through  40  random 
solitons.  The  solitons  start  about  2.5  km  from  the  source  with  an  average  separation 
of  450  m  between  soliton  leading  edges.  Time  zero  corresponds  to  the  PAM1  arrival. 

Soliton  Number  Effects 

Ten  computer  runs  of  100  realizations  each  were  conducted  to  evaluate  the  effects  of 
adding  a  number  of  solitons  between  the  source  and  receiver.  The  200  m  solitons  are 
randomly  placed  with  an  average  separation  of  250  m  between  solitons,  starting  at 
an  average  distance  of  2.5  km  from  the  source.  Figure  5-6  shows  that  the  addition  of 
solitons  is  clearly  accompanied  by  an  increase  in  the  probability  that  dispersed  mode 
1  energy  will  dominate  the  arrival  pattern.  The  averages  of  spread  and  bias  statistics 
taken  over  the  100  trials  are  displayed  as  function  of  the  number  of  solitons  in  Figure 
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Figure  5-7:  The  average  mean  (-),  median  (-  -),  and  peak  (•  •  •)  lag  with  respect  to 
the  pseudo-adiabatic  arrival,  and  standard  deviation  (-•)  of  scattered  mode  1  arrivals 
plotted  as  a  function  of  the  number  of  solitons  placed  between  the  source  and  the 
receiver,  where  the  solitons  are  added  starting  from  the  source  end.  The  upper  three 
lines  are  all  bias  measures,  and  the  lower  line  is  a  measure  of  spread. 


5-7.  Note  the  close  agreement  between  all  three  measures  of  signal  bias;  the  small  lags 
between  the  three  are  indicative  of  a  fairly  regular  mean  pulse  shape  with  a  longer 
tail  on  the  trailing  edge. 

Figure  5-8  shows  the  situation  to  be  entirely  different  if  the  solitons  are  added 
starting  from  the  receiver  end.  It  contains  a  re-plot  of  the  data  in  Figure  5-7  plus 
additional  curves  for  averages  of  spread  and  bias  statistics  when  the  solitons  are  added 
starting  near  the  receiver. 

Adding  a  few  solitons  near  the  receiver  dramatically  affects  the  spread  of  mode  1 
arrivals,  building  quickly  to  a  peak  spread  for  about  20  solitons.  The  effect  on  bias  is 
also  dramatic  when  measured  by  the  mean  arrival,  but  it  is  less  so  when  the  median  is 
used.  Both  of  these  measures  of  bias  seem  to  be  approaching  asymptotic  levels  at  the 
25  soliton  point.  When  peak  position  for  less  than  15  solitons  is  considered,  no  real 
difference  can  be  seen  between  the  "adding  near  the  source"  and  the  "adding  near 
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Figure  5-8:  The  average  mean  (-),  median  (-  -),  and  peak  (•  •  •)  lag  with  respect  to 
the  pseudo-adiabatic  arrival,  and  standard  deviation  (-•)  of  scattered  mode  1  arrivals 
plotted  as  a  function  of  the  number  of  solitons  placed  between  the  source  and  the 
receiver.  The  upper  curves,  averages  over  50  realizations,  represent  the  case  where 
the  solitons  are  added  starting  at  the  receiver  end.  The  lower  set  of  curves  is  replotted 
from  figure  5-5.  60  solitons  is  essentially  saturation,  and  the  same  points  are  used  in 
both  sets  of  curves. 
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the  receiver"  results.  The  large  differences  between  mean,  median,  and  peak  position 
for  moderate  numbers  of  near  receiver  solitons  are  indicative  of  a  mean  pulse  shape 
with  a  substantial  tail  on  the  trailing  edge. 

All  these  effects  can  be  understood,  when  one  considers  the  amount  of  modal 
dispersion  produced  by  this  waveguide  over  a  distance  of  32  km  (see  Figure  4-9).  If 
scattering  interactions  are  constrained  to  occur  only  near  the  receiver,  the  peak  of 
the  mode  1  arrival  is  difficult  to  affect,  but  the  tail  of  that  arrival  will  be  dramatically 
altered.  The  mean  arrival  time  and  the  standard  deviation  of  the  arrival  envelope  are 
very  sensitive  to  changes  in  the  tail.  The  median  and  IQR  of  the  arrival  envelope  are 
less  sensitive,  and  the  peak  is  by  definition  insensitive  to  the  tail,  unless,  of  course, 
the  tail  becomes  the  peak. 

The  asymptotic  level  of  bias  reached  in  the  "adding  near  the  receiver"  case  is  to  be 
expected.  With  no  scattering,  the  history  of  the  mode  1  arrival  (single  subarrival)  is 
mode  1  only.  As  the  number  of  scattering  events  is  increased  the  mean  modenumber 
history  is  weighted  more  and  more  by  higher  modes,  but  there  is  a  limiting  value  for 
the  average  percentage  of  mode  n  history  in  a  mode  1  arrival. 

Dozier  and  Tappert  [15]  have  shown  in  a  simplified  model  of  scattering  with  no 
modal  attenuation  that  the  equilibrium  energy  state  is  an  equipartition  of  energy 
between  the  modes.  The  asymptotic  limit  for  mode  1  modenumber  history  will, 
therefore,  approach  equal  percentages,  with  some  additional  weighting  (that  depends 
upon  the  average  distance  required  to  randomly  redistribute  the  energy  in  mode  1, 
relative  to  the  total  distance  traveled)  given  to  mode  1.  Similar  arguments  hold  for 
the  other  modes.  In  more  realistic  scenarios,  the  asymptotic  mode  1  modenumber 
history  percentages  will  still  approach  the  equilibrium  energy  partitions,  but,  as  shown 
by  Creamer  [18],  with  attenuation  these  will  not  be  equal  percentages.  The  energy 
partitions  and  modenumber  histories  will  be  weighted  towards  the  modes  with  the 
least  attenuation. 

The  soliton  effects  on  spread  can  be  more  easily  visualized  by  considering  the 
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most  simplistic  case  of  two  modes  that  sequentially  scatter  with  all  four  coupling 
coefficients  equal  to  0.5.  For  equal  mode  1  and  mode  2  intensities  and  equally  spaced 
scattering  events,  with  the  last  occurring  just  before  the  receiver,  this  scenario  will 
yield  an  average  mode  1  subarrival  modenumber  history  of  50  percent  travel  in  mode 
1  and  50  percent  travel  in  mode  2.  There  will  always  be  some  variance  between  the 
individual  subarrival  histories  with  some  subarrivals  spending  more  time  in  mode 
1  and  others  more  time  in  mode  2,  but  the  spread  will  decrease  as  the  number  of 
scattering  events  increases.  For  narrow  pulses  with  adiabatic  mode  1  and  mode  2 
travel-times  separated  by  AT,  the  spread  of  a  mode  1  arrival,  as  measured  by  the  N 
normalized  arrival  standard  deviation  is 

AT 

where  n  is  the  number  of  scattering  events  and  N  =  2n  is  the  number  of  subarrivals. 
Of  the  N  mode  1  subarrivals,  one  is  the  PAM1  arrival,  a  second  is  coincident  with 
PAM2  arrival  (scattered  off  right  before  the  receiver),  and  the  remaining  2n  —  2  arrive 
in  A;-groups,  where  k  is  number  of  increments  traveled  in  mode  2.  The  number  of 
subarrivals  in  a  fc-group  is  given  by  the  familiar  formula, 

(5.2) 


fc!(n-fc)!' 

for  finding  the  number  of  ^-combinations  of  n  distinguishable  items  [19].  The  PAM1 
and  "coincident  with  PAM2"  arrivals  can  also  be  considered  as  fc-groups,  with  k  equal 
to  0  and  n,  respectively. 

Table  5.1  summarizes  the  subarrival  statistics  for  the  case  where  n  =  4  and  AT  = 
1.  Clearly,  the  same  scenario  with  only  one  increment  (i.e.  one  scattering  at  the  very 
end)  will  still  yield  an  average  bias  value  of  .50,  but  the  spread  will  be  maximized 
(also  equal  to  .50).  Any  additional  scattering  prior  to  this  late  event  decreases  the 
spread,  with  no  change  in  bias.    Conversely,  the  two-mode  case  with  only  a  single 
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Simple  Two  Mode  Scattering  Scenario 
with  Four  Evenly  Spaced  Increments 

Subarrival  k-group 

Total 

0 

1 

2 

3 

4 

Number  in  Group 

1 

4 

6 

4 

1 

16 

Group  Bias 

0 

.25 

.50 

.75 

1.00 

.50 

Group  Spread 

0 

0 

0 

0 

0 

.25 

Table  5.1:  Subarrival  statistics  for  the  simple  two  mode  50/50  scattering  scenario 
where  there  are  four  equally  spaced  increments  and  the  difference  between  the  PAM1 
and  PAM2  arrival  times  is  AT  =  1. 

scatter  occurring  early  in  the  path  will  produce  a  small  amount  of  spread  and  bias; 
additional  scattering  events  occurring  later  will  cause  an  increases  in  both. 

The  Figure  5-8  simulation  runs  are  more  complicated  than  the  above  two-mode 
example,  but  the  same  effects  can  be  seen.  The  eight-mode  cross-modal  coupling 
coefficients  are  much  less  than  0.5,  so  the  "adding  near  the  receiver"  peak  spread 
is  not  reached  until  around  20  solitons.  We  can  also  see  the  bias  stabilize  at  about 
the  same  number  of  solitons,  indicating  saturation  in  mean  modenumber  history  of 
subarrivals.  The  two-mode  "adding  near  the  receiver"  example  was  saturated  from 
the  outset. 


Soliton  Location  Effects 

The  above  section  has  illustrated,  via  Figure  5-8  in  particular,  that  the  general  lo- 
cation of  solitons  in  a  waveguide  has  more  impact  on  spread  and  bias  levels  than 
does  the  number  of  solitons.  In  this  section  we  will  show  that,  in  addition  to  the 
general  location,  the  specific  location  of  solitons  contributes  significantly  to  the  na- 
ture of  mode  1  arrival  patterns.  Figure  5-9  demonstrates  this  with  two  trains  of  six 
SIA  solitons  propagating  in  lock-step  down  a  waveguide  at  0.8  m/s.  As  the  soliton 
trains  progress  towards  the  source,  the  location  of  the  last  scattering  event  moves 
towards  the  source,  reducing  the  potential  for  mode  1  arrival  bias.  As  expected,  this 
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general  trend  in  mean  and  peak  arrival  bias  is  evident  in  Figure  5-9,  but,  even  more 
interesting  is  the  fact  that  the  mean  and  peak  arrival  time  bias  is  very  sensitive  to 
small  shifts  in  soliton  train  positions. 

Figure  5-10,  a  companion  to  Figure  5-9,  provides  a  good  qualitative  sense  of  the 
mode  1  arrival  pattern's  dependence  on  the  position  of  the  soliton  trains.  The  outline 
of  the  first  arrival  at  the  bottom  of  the  figure  shows  four  distinct  humps  that  roughly 
correspond  in  location  to  the  mode  1-4  arrivals  shown  in  Figure  4-9.  For  the  majority 
of  the  arrivals  (see  Figure  5-9),  the  peak  arrival  resides  on  the  mode  1  hump,  but  it 
often,  almost  periodically,  shifts  to  the  mode  3  hump,  and  it  occasionally  resides  on  the 
mode  2  hump.  The  mean  arrival  fluctuations  are  equally  dramatic.  The  sensitivity  of 
the  arrival  shapes  and  bias  measurements  to  small  changes  in  the  location  of  solitons 
is  chiefly  due  to  the  constructive  and  destructive  interference  patterns  of  scattered 
subarrivals.  The  undulations  evident  in  the  "mode  1"  hump  on  Figure  5-10  contribute 
significantly  to  the  peak  position  jumps  in  Figure  5-9.  The  source  of  these  undulations 
is  alternating  constructive  and  destructive  interference  between  mode  1  and  2  in  the 
first  soliton  train.  The  mode  cycle  distance  between  mode  1  and  2  is  about  1420 
m;  this  combined  with  a  train  speed  of  48  m  per  minute,  yields  the  half  hour  cycles 
evident  in  both  figures. 

Figure  5-11,  a  companion  to  Figure  5-10,  shows  two  snapshots  of  the  arrival  shapes 
for  all  eight  modes  propagated  in  this  simulation.  The  left  hand  column  shows  the 
minute  0  arrival  shapes,  and  the  right  column  shows  the  minute  300  arrival  shapes.  A 
comparison  of  the  two  columns  shows  the  movement  of  the  two  soliton  trains  towards 
the  source  markedly  reduces  the  "spread  level"  for  all  eight  modes.  The  thinner  lines  in 
this  figure  represent  predicted  arrival  shapes  for  modes  that  propagate  adiabatically 
(i.e.  no  scattering  of  energy  between  modes).  The  eighth  mode  is  only  minimally 
excited  by  the  source  at  27  m  depth,  and  it  undergoes  the  most  attenuation,  so  it 
is  not  even  visible  as  an  unscattered  arrival.  There  is  little  modal  dispersion  beyond 
mode  6,  and  as  a  result,  the  scattered  arrivals  for  modes  6-8  show  only  negative  spread 
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SIA  Soliton  Waveguide  (0  Minutes) 
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Figure  5-9:  Time  dependent  waveguide  with  two  sets  of  six  SIA  solitons.  The  upper 
frame  shows  the  starting  positions,  the  middle  frame  shows  the  positions  five  hours 
later,  and  the  lower  frame  is  plot  of  peak  position  bias  (re  PAM1)  as  a  function 
of  time.  The  solitons  propagate  in  lock-step,  traversing  48  m  per  minute,  the  time 
between  each  peak  position  data  point  (o).  The  mean  arrival  data  is  also  in  one 
minute  increments,  but  the  points  are  connected  by  a  solid  line. 
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Mode  1  Arrival  Pattern 
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Figure  5-10:  Waterfall  of  mode  1  arrival  patterns  from  the  time  dependent  waveguide 
with  two  sets  of  six  SIA  solitons.  The  300  peak  positions  in  the  lower  frame  of  Figure 
5-9  were  picked  from  these  arrival  patterns. 
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Figure  5-11:  Snapshots  of  the  arrival  shapes  for  all  eight  modes  propagated  in  the 
Figure  5-10  simulation.  The  left  hand  column  shows  the  minute  1  arrival  shapes 
and  the  right  column  shows  the  minute  300  arrival  shapes.  The  thicker  lines  are  the 
scattered  arrival  shapes,  and  the  thinner  lines  show  the  corresponding  unscattered 
mode  arrivals.  The  eighth  mode  is  only  minimally  excited  by  the  source  at  27  m  depth, 
and  it  undergoes  the  most  attenuation,  so  it  is  not  even  visible  as  an  unscattered 
arrival. 
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Figure  5-12:  A  comparison  of  mean  (-)  and  peak  (o)  arrival  time  as  function  of  the 
number  of  evenly  spaced  solitons  in  a  32  km  waveguide.  The  first  soliton  is  always 
200  m  from  the  source  and  the  last  soliton  is  roughly  32  km  divided  by  the  number 
of  solitons  away  from  the  receiver. 


in  travel  time. 

Another  example  of  the  sensitivity  of  mode  1  bias  to  specific  soliton  location  can 
be  observed  by  considering  mode  1  bias  level  as  a  function  of  the  number  of  evenly 
spaced  200  m  SIA  solitons  in  a  32  km  waveguide.  The  first  soliton  is  set  at  200  m  from 
the  source,  the  rest  are  evenly  spaced  between  there  and  the  receiver.  An  additional 
soliton  in  the  arrangement  creates  a  series  of  small  changes  in  soliton  location  and 
spacing;  their  cascading  effects  cause  the  large  changes  in  the  mode  1  arrival  bias 
plotted  in  Figure  5-12. 
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Source  Excitation  Effects 

We  will  now  look  at  one  final  contributor  to  mode  arrival  shape  fluctuations.  In 
a  conceptual  sense,  the  mode  excitation  at  the  source  determines  the  size  of  the 
individual  mode  pulses  leaving  the  receiver;  these  pulse  sizes  in  turn  control  the  size 
of  all  the  descendent  subarrivals,  which  will  in  turn  combine  to  form  the  mode  arrivals 
at  the  receiver. 

Solitons  routinely  cause  upwards  of  10  m  depressions  in  the  thermocline;  the  effects 
of  these  depressions  can  be  seen  in  the  upper  frame  of  Figure  5-13  which  shows  a  2 
hr  record  from  the  shallowest  thermistor  on  the  WVLA.  The  source  excitation  levels 
(source  depth  of  29  m)  in  the  middle  frame  were  obtained  from  a  series  of  mode 
shapes  derived  from  the  same  2  hrs  of  WVLA  thermistor  data.  A  depth  of  68  m  was 
used  for  the  mode  shape  calculations.  The  waveguide  for  this  simulation  is  similar 
to  the  ones  shown  in  Figure  5-9,  except  it  is  frozen  (i.e  it  does  not  move);  the  source 
excitation  is  the  only  parameter  that  changes.  The  leading  edge  of  the  first  train  was 
fixed  at  5  km  from  the  source  and  the  trailing  edge  in  the  second  train  was  fixed  at 
5  km  from  the  receiver.  The  resulting  shifts  in  the  mean  arrival  are  quite  dramatic, 
but  the  mode  excitation  changes  are  insufficient  to  shift  the  peak  significantly. 

5.2      Analysis  of  the  "Propagated  Thermistor  String 
Model" 

5.2.1      An  Example  Waveguide 

The  fine  scale  "propagated  thermistor  string"  model  described  in  Chapter  4  is  much 
slower  than  the  SIA  model,  but  the  use  of  thermistor  string  data  brings  more  realistic 
size  and  spacing  characteristics  for  solitons,  as  well  as  added  scattering  from  other 
sources  (e.g.  linear  internal  waves).  An  example  of  this  type  of  waveguide  overlaid 
with  isopycnal  surface  displacements  taken  from  a  high  frequency  backscatter  record 
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Figure  5-13:  The  top  graph  shows  12.5  m  depth  temperature  readings  from  the  2  hrs 
of  thermistor  data  used  to  simulate  modal  excitation  levels.  The  middle  graph  shows 
the  resulting  modal  excitation  levels  for  mode  1  (solid),  mode  2  (dashed),  mode  3 
(dot-dash),  and  mode  4  (dotted).  Mode  5  is  similar  in  magnitude  and  used  as  an 
input  to  the  model,  but  it  is  omitted  from  this  figure  to  avoid  clutter.  The  lower 
graph  shows  the  output  bias  levels  corresponding  to  the  excitation  inputs  shown  in 
the  middle  graph.  The  solid  line  is  the  mode  1  arrival  envelope  peak  arrival  time  in 
excess  of  the  PAM1  arrival  time  and  the  dotted  line  is  the  mode  1  arrival  envelope 
mean  arrival  time  in  excess  of  the  PAM1  arrival  time. 
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provided  by  Marshall  Orr  (NRL,  Washington  D.C.,  private  communication)  is  shown 
in  Figure  5-14.  The  waveguide  in  this  case  is  produced  from  records  from  a  dedicated 
thermistor  string  (not  the  WVLA  thermistors)  used  during  SWARM.  This  "T-string 
598"  data  was  used  in  this  comparison  with  the  Orr  backscatter  record,  because  the 
WVLA  was  not  in  the  water  early  enough  to  record  these  events.  Various  comparisons 
between  different  instruments  point  to  non-linear  internal  wave  propagation  speeds 
in  the  range  of  0.6  to  0.8  m/s.  The  propagation  speed  of  0.62  m/s  used  in  Figure 
5-14  was  chosen  for  best  overall  agreement  between  the  records.  The  agreement  is,  in 
fact,  quite  good  —  better  than  one  would  have  expected  using  the  "constant  velocity" 
assumption  for  soliton  and  internal  tide  propagation. 

5.2.2      Model  Results 

The  WVLA,  with  its  thermistors  spanning  12.5  m  to  60.5  m  in  depth,  provides  better 
waveguides  for  simulating  tomographic  transmissions.  The  sound  speeds  are  extrap- 
olated for  depths  above  12.5m  and  below  60.5  m,  and  by  assuming  a  WVLA  temper- 
ature signal  propagation  speed  of  .8  m/s  along  the  acoustic  path,  each  670  minutes 
of  thermistor  string  data  is  transformed  into  a  32160  by  68  m  waveguide.  The  profile 
at  the  source  determines  the  "eight  mode"  starting  field  for  the  model,  and  an  "eight 
mode"  transmission  is  then  propagated  through  the  waveguide  and  relevant  statistics 
are  obtained  from  the  mode  arrivals. 

Peak  Arrivals 

Figure  5-15  shows  several  days  of  mode  1  peak  arrival  times  in  both  absolute  and 
wander-corrected  terms.  There  are  six  minutes  between  each  acoustic  transmission, 
and,  in  order  to  provide  some  qualitative  sense  of  the  variability  in  the  strength  of 
the  arriving  signals,  the  size  of  each  arrival  dot  is  scaled  by  the  amplitude  of  the  peak 
arrival.  The  peak  arrivals  consistently  display  a  tendency  towards  positive  travel 
time  bias  relative  to  the  PAM1  arrival,  but  there  is  a  notable  exception  for  the  period 
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T-string  598  and  Orr  Backscatter  Record  (July  30,  0330-0728Z) 
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Figure  5-14:  Pycnocline  displacements  from  an  Orr  backscatter  record  compared  to 
a  thermistor  string  record  using  a  propagation  speed  of  .62  m/s  over  the  distance 
between  the  back  scatter  instrument  and  the  thermistor  string  at  each  backscatter 
reading. 
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Figure  5-15:  Distribution  of  peaks  for  the  mode  1  arrivals  of  simulated  transmis- 
sions using  sequential  WVLA  thermistor  string  generated  waveguides.  The  sequen- 
tial 32160  by  68  m  waveguides  are  separated  by  six  minutes.  The  size  of  each  arrival 
dot  is  scaled  by  the  amplitude  of  the  peak  arrival.  The  upper  frame  includes  a  plot 
of  the  PAM1  arrival  time  in  conjunction  with  the  individual  peak  arrival  times;  the 
zero  reference  on  the  arrival  time  axis  corresponds  to  21.6  seconds.  The  lower  frame 
plots  the  bias  or  offset  of  the  peak  arrivals  relative  to  the  PAM1  arrival  time  for  each 
realization. 


104 


around  day  212.5.  During  this  period,  the  PAM2  arrival  is  actually  earlier  than  the 
PAM1  arrival,  allowing  for  negative  bias  in  some  mode  1  arrival  times. 

For  periods  where  the  PAM1  arrival  is  the  earliest  arrival,  the  distributions  of 
mode  1  peak  arrivals  tend  to  exhibit  a  fairly  clean  leading  edge  that  coincides  with 
the  PAM1  arrival  time.  The  arrivals  with  the  largest  peak  amplitudes  tend  to  fall 
within  10  to  20  ms  of  this  leading  edge.  Arrivals  with  peak  bias  levels  in  the  40 
to  60  ms  range  are  partially  the  result  of  destructive  interference  among  the  early 
subarrivals,  and  as  a  result,  the  peak  heights  are  noticeably  smaller. 

The  story  for  mode  2  arrivals  is  somewhat  different.  Figure  5-16  displays  the 
peak  mode  2  arrival  distributions  obtained  from  the  simulated  transmissions.  The 
temporal  variations  in  peak  amplitude  are  much  stronger  for  mode  2,  while  the  cor- 
relation between  bias  level  and  peak  amplitude  appears  weaker  relative  to  the  mode 
1  relationship.  Also,  save  for  the  period  where  the  PAM2  arrival  is  earlier  than  the 
PAM1  arrival,  we  generally  see  significant  peak  arrival  time  bias  in  both  the  negative 
and  positive  directions,  so,  unlike  the  PAM1  arrival  time,  the  PAM2  arrival  time 
is  not  readily  discernible  from  the  total  distribution  of  mode  2  peak  arrivals.  The 
PAM1  arrival  was  generally  discernible,  because  it  was  usually  "smeared"  in  only  one 
direction. 

A  simple  "thought  experiment"  may  prove  useful.  Consider  a  curve  thickly  drawn 
in  wet  paint  across  a  clean  white  surface.  Now  from  one  side  of  the  curve,  direct  a 
high  velocity  wind  across  its  surface.  The  paint  will  smear,  but  only  in  one  direction; 
the  original  shape  of  the  curve  is  still  discernible.  However,  if  the  winds  are  instead 
directed  in  a  random  fashion  from  both  sides  of  the  curve,  the  original  shape  will  be 
lost. 

Therein  lies  the  problem  with  the  mode  2  and  higher  arrivals.  They  are  "smeared" 
in  both  directions,  severely  limiting  their  utility  in  standard  tomographic  inversions 
for  range-average  soundspeed  profiles.  However,  the  smearing  levels  themselves  may 
prove  to  be  a  useful  input  to  an  inversion  for  something  else  in  the  future  (e.g.  location 
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Figure  5-16:  Distribution  of  peaks  for  the  mode  2  arrivals  corresponding  to  the  mode  1 
arrivals  of  Figure  5-15.  The  upper  frame  includes  a  plot  of  both  the  PAM1  arrival  time 
(dotted  line)  and  PAM2  arrival  time  (solid  line),  in  conjunction  with  the  individual 
peak  arrival  times;  the  zero  reference  on  the  arrival  time  axis  corresponds  to  21.6 
seconds.  The  lower  frame  plots  the  bias  or  offset  of  the  mode  2  peak  arrivals  relative 
to  the  PAM2  arrival  time  for  each  realization. 
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of  a  soliton  train). 

Spread  Statistics 

The  interquartile  range  (IQR)  of  a  mode  arrival  envelope  provides  a  robust  estimate 
for  that  mode's  time  spreading.  Figure  5-17  plots  these  spread  statistics  for  mode  1 
and  2,  in  30  minute  bin  averages,  as  an  overlay  with  the  WVLA  temperature  profiles 
for  the  simulation  test  period.  For  a  given  spread  reading  the  coincident  temperature 
profile  was  used  in  modeling  the  SSP  at  the  receiver.  The  profile  at  a  time  0.46  days 
to  the  left  was  used  at  the  source,  and  the  profiles  in  between  were  used  in  the  interior 
of  the  modeled  32,000  m  waveguide. 

On  average,  larger  amounts  of  spread  tend  to  be  accompanied  by  a  significant 
amount  of  soliton  activity  in  the  latter  half  of  the  waveguide  (i.e.  in  the  region  from 
the  time  in  question  to  about  0.25  days  to  left),  but  this  is  not  always  the  case.  The 
large  peaks  in  the  the  mode  1  spread  history  around  yeardays  213.4  and  214.8,  for 
example,  are  almost  inexplicable  in  terms  of  the  accompanying  soliton  activity  in  the 
waveguide  or  the  excitation  levels  at  the  source.  These  are  further  examples  of  the 
sensitivity  of  spread  and  bias  statistics  to  the  details  of  the  modeled  waveguide.  Mode 
2  spread,  on  the  other  hand,  seems  to  be  more  closely  tethered  to  the  "general  level" 
of  soliton  activity,  showing  a  near  tidal  periodicity  in  spread  fluctuations. 

One  other  notable  region  of  divergence  between  the  two  spread  histories  is  in 
the  region  around  day  212.5.  Mode  1  shows  moderate  spread,  while  the  spread  of 
mode  2  is  at  near  minimum,  and  the  corresponding  soliton  activity  is  low  to  moderate. 
Fluctuations  in  soliton  activity  are  probably  not  the  main  source  of  spread  fluctuations 
during  this  period,  rather  Figure  5-15  reveals  this  to  be  the  region  where  the  PAM2 
arrival  time  is  shifting  significantly  relative  to  the  PAM1  arrival  time. 

The  correlation  between  spread  levels  and  soliton  activity  is  difficult  to  quantify; 
there  is  no  single  statistical  measure  of  soliton  activity.  If  the  standard  deviation  of 
the  temperature  at  22.5  m  depth  over  the  previous  4.3  hrs  is  used  as  a  crude  measure 
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Figure  5-17:  30  minute  bin  averaged  IQR  measurements  and  the  WVLA  temperature 
profile  for  the  simulation  test  period.  The  profiled  temperatures  range  from  7.2  to 
25.9  degrees  C.  The  top  frame  is  for  mode  1  and  the  bottom  is  for  mode2. 
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of  soliton  activity  near  the  receiver,  then  an  attempt  at  measuring  the  correlation  can 
be  made.  Unlike  spread  levels,  this  statistic  cannot  change  appreciably  from  minute 
to  minute,  but  the  long  (~  12  km)  reach  is  necessary  to  encompass  the  most  critical 
of  region  soliton  scattering.  Figure  5-18  shows  1  hr  bin  averaged  spread  statistics 
plotted  against  1  hr  bin  averages  of  this  statistic  for  modes  2-4.  Here  we  see  similar 
rise  and  fall  patterns  in  all  the  statistics,  and  computed  zero-lag  cross-correlations 
between  the  temperature  statistic  and  the  mode  2-4  spread  levels  are  .43,  .14,  and 
.17  respectively  (all  relative  peaks  in  correlation,  except  mode  4  which  peaks  at  .34 
with  a  3  hr  lag).  Similarly,  the  zero-lag  cross-correlation  between  the  temperature 
statistic  and  mode  1  spread  is  .34. 

The  search  for  a  better  general  soliton  activity  statistic,  in  terms  of  correlation 
with  spread  levels,  will  be  difficult;  there  are  several  contributing  factors  to  consider. 
Small  changes  in  the  position  and  spacing  of  solitons,  as  well  as  variations  in  mode 
excitation  at  the  source,  can  have  a  significant  impact  on  spread  statistics,  and  almost 
any  general  measure  of  soliton  activity  imaginable  will  only  change  appreciably  over 
hourly  or  greater  time  scales.  Spread  levels  in  both  simulated  and  actual  data  can 
change  appreciably  in  less  than  one  minute. 

As  observed  in  Figure  5-15,  increased  bias  in  mode  1  peak  position  is  generally 
accompanied  by  a  reduction  in  the  mode  1  peak  height.  The  same  relation  holds 
between  spread  and  peak  height,  and  the  correlation  is  fairly  strong,  as  evidenced  by 
Figure  5-19.  This  figure  also  shows  that  peak  height  fluctuations  for  mode  2  have  a 
negative  a  zero-lag  correlation  with  mode  1  spread  statistics,  presumably  the  result 
of  a  corresponding  "smearing"  of  mode  2  energy.  However,  this  negative  correlation 
is  diminished  as  the  mode  number  increases.  This  sort  of  analysis  can  be  performed 
using  the  spread  statistics  from  higher  modes;  a  tabulations  of  the  zero-lag  results  are 
provided  in  Table  5.2.  Not  unexpectedly,  the  amounts  of  correlation  vary  considerably 
between  modes,  reflecting  the  complicated  relationship  between  the  degree  of  modal 
coupling  and  the  amount  of  spread  and  peak  intensity  in  higher  mode  number  arrivals. 
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Figure  5-18:  A  comparison  of  1  hour  bin  averaged  IQR  levels  for  the  simulated 
mode  2-4  data  (solid  lines)  and  a  1  hr  bin  average  of  the  Standard  Deviation  of 
the  Temperature  at  22.5  m  depth  over  the  previous  4.3  hrs  (dashed  lines).  The 
standard  deviation  of  the  temperature  is  in  tenths  of  a  degree  centigrade  to  facilitate 
comparison  with  the  IQR  levels  in  milliseconds. 
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Mode  1  -4  Peak  Height  vs.  Mode  1  Spread 
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Figure  5-19:  Cross-correlations  of  1.4  hr  averaged  peak  height  fluctuations  for  mode 
1  (solid),  2  (dashed),  3  (dot-dashed),  and  4  (dotted)  with  the  corresponding  1.4  hr 
averaged  mode  1  IQR  (spread)  levels. 


Zero- Lag  Cross- Correlations 

Spread 
Mode 

Peak  Height  Mode 

1 

2 

3 

4 

1 

-.67 

-.36 

-.07 

.06 

2 

-.27 

-.72 

-.35 

-.17 

3 

.10 

-.46 

-.62 

-.21 

4 

.30 

.04 

.23 

-.26 

Table  5.2:  Zero-lag  cross  correlations  between  1.4  hour  bin  averaged  mode  1-4  peak 
heights  and  mode  1-4  mean  spread  levels. 
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5.3      The  Transition  to  Real  Data 

A  problem  with  conversion  from  simulations  to  real  data  analysis  is  a  loss  of  the 
knowledge  of  the  exact  location  of  the  PAM1  arrival.  In  the  simulation  analysis  this 
location  was  easily  calculated,  because  the  environments  were  deterministic.  In  real 
data,  the  location  of  this  arrival  is  unknown.  Furthermore,  its  location  will  change 
from  sequence  to  sequence.  This  movement  is  referred  to  as  the  wander  of  the  mode  1 
arrival.  Some  small  amounts  of  wander  may  occur  over  short  time  scales,  but  wander 
signals  are  generally  dominated  by  large  scale  changes  in  temperature  or  current. 
These  carry  with  them  longer  time  scales  than  the  soliton  scattering  that  dominates 
the  spread  and  bias  of  the  mode  1  arrival.  In  our  shallow  waveguide,  mode  1  energy 
is  confined  to  the  cold  slow  region  below  the  pycnocline  (see  Figures  4-6  and  4-7), 
so  passage  of  a  soliton  tends  to  push  mode  1  deeper.  The  deepening  of  the  mode  1 
shape  almost  exactly  compensates  for  a  simultaneous  increase  in  mean  sound-speed, 
such  that  the  mode  1  group  velocity  is  relatively  insensitive  to  the  solitons  in  the 
waveguide.  This  translates  to  a  very  low  mode  1  soliton  wander  signal.  In  our 
SIA  model,  a  change  from  zero  to  30  solitons  caused  only  about  a  1  ms  drop  in  the 
PAM1  travel  time.  This  compensating  behavior  combined  with  the  fact  that  the  most 
energetic  linear  internal  waves  are  confined  to  the  upper  water  column,  should  keep 
most  sources  of  wander,  at  soliton  time-scales,  small  compared  to  the  bias  signals 
caused  by  soliton  induced  scattering.  The  difference  in  time  scales  between  the  major 
sources  of  wander  and  major  sources  of  bias  provides  the  potential  for  separating 
them! 

As  illustrated  by  the  relationship  between  the  peak  arrival  distribution  leading 
edge  and  the  PAM1  arrival  curve  plotted  in  Figure  5-15,  estimating  the  PAM1  arrival 
time  should  not  be  too  difficult.  The  histograms  of  SIA  model  peak  positions  (Figure 
5-6)  provide  some  further  insight;  the  PAM1  arrival  is  always  in  one  of  the  first  three 
bins,  usually  the  second,  so  for  a  group  of  100  random  realizations  of  a  simulated 
waveguide,  the  earliest  peak  sample  time  (plus  one  sample)  may  be  a  reasonably 
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good  estimate  for  the  PAM1  arrival. 

The  applicability  of  this  estimate  to  the  wander  of  real  mode  1  arrivals  does  not 
follow  directly.  The  SWARM  data  has  22  sequences  every  six  minutes,  so  approxi- 
mately 200  sequences  are  available  each  hour.  However,  many  of  the  sequences  will 
be  correlated.  The  measured  correlation  times  are  between  70  and  90  seconds,  so 
our  200  samples  of  actual  data  may  only  be  equivalent  15  or  20  random  realizations. 
However,  this  should  be  sufficient;  PAM1  estimation  is  clearly  possible  with  the  10 
samples  per  hour  data  rate  shown  in  Figure  5-15. 

Estimating  the  pseudo-adiabatic  arrival  times  (i.e.  wander)  for  higher  modes  will 
be  more  problematic.  This  makes  the  determination  of  bias  levels  for  higher  modes 
difficult  and  imprecise.  However,  the  measures  of  spread  used  in  the  simulated  data 
are  not  directly  dependent  on  knowing  either  the  wander  or  the  bias,  so  the  transition 
to  real  data  analysis  for  this  variable  is  straightforward. 
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Chapter  6 


Analysis  of  Real  Data 


6.1      400  Hz  Mode  1  Arrival  Analysis 

6.1.1      Peak  Arrival  Time 

The  nature  of  the  SWARM  waveguide  is  such  that  mode  1  generally  has  the  highest 
group  velocity,  making  the  PAM1  arrival  the  earliest  possible  arrival  for  most  trans- 
missions. There  are  exceptions.  During  the  early  hours  of  day  212,  at  the  WVLA, 
mode  2  has  a  higher  group  velocity  than  mode  1;  if  and  when  this  is  the  case  for  the 
majority  of  the  waveguide,  the  PAM2  arrival  will  precede  the  PAM1  arrival,  allowing 
for  some  positively  biased  mode  1  energy,  as  seen  in  Figure  5-15. 

The  distribution  of  SWARM  mode  1  peak  arrival  times  shown  in  Figure  6-1  has 
a  fairly  clean  leading  edge;  there  is  little  or  no  indication  of  negatively  biased  peaks. 
This  strongly  suggests  that  the  PAM1  arrival  was  the  earliest  arrival  over  most,  if 
not  all,  of  the  period  under  evaluation.  The  pattern  of  the  PAM1  arrival  time  (i.e. 
the  leading  edge)  shows  a  clear  semidiurnal  tidal  signal,  as  well  as  significant  wander 
at  even  lower  frequencies.  The  two  strongest  and  most  probable  sources  for  these 
fluctuations  are  barotropic  tidal  currents  along  the  path  between  the  source  and 
receiver  and  the  tidal  advection  of  deeper  warm  water  (e.g.  the  foot  of  the  shelfbreak 
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Mode  1  Peak  Arrival  Distribution 
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Figure  6-1:  Distribution  of  SWARM  Mode  1  Peak  arrival  times.  The  zero  reference 
is  at  an  arrival  time  of  21.985  seconds.  The  plotted  transmissions  occurred  every  six 
minutes,  with  each  transmission  producing  22  sequence  arrivals  over  a  two  minute  pe- 
riod. The  location  of  the  highest  peak  in  each  sequence  is  plotted.  Though  not  evident 
in  most  cases,  there  are  gaps  in  the  record  that  range  from  single  lost  transmissions 
to  occasional  larger  gaps  of  up  to  two  hours  in  duration. 
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front)  in  and  out  of  the  acoustic  path. 

Near  the  WVLA,  Mode  1  energy  is  for  the  most  part  confined  to  the  slower  sound 
speed  regions  in  the  40  to  60  m  depth  band  (see  Figure  4-7).  An  upward  looking 
Acoustic  Doppler  Current  Profiler  (ADCP)  moored  near  the  WVLA  in  75  m  of  water 
provides  a  measure  of  the  currents  in  the  region.  Taking  the  average  current  in  five  4  m 
depth  bins  between  39  and  59  m  depth,  low-pass  filtering,  and  integrating  the  effects 
on  travel  time  over  the  entire  32,000  m  acoustic  path  (assuming  range  independent 
currents)  yields  the  solid  curve  shown  in  Figure  6-2.  Using  14  bins  between  11  and  67 
m  depth  produces  the  dotted  line.  The  similarity  between  the  two  curves  shows  that 
the  vast  majority  of  the  current  signal  in  the  depth  band  of  interest  (solid  curve)  is 
barotropic  (baroclinic  currents  from  solitons  are  removed  by  the  lowpass  filter).  The 
barotropic  current  clearly  accounts  for  some  of  the  semidiurnal  PAM1  wander,  but  it 
contains  negligible  lower  frequency  signals. 

Subtracting  the  current  signal  from  the  mode  1  peak  arrival  times  yields  the 
residual  distribution  of  peak  arrival  times  plotted  in  Figure  6-3.  The  residual  mode 
1  wander  at  the  semidiurnal  tidal  frequency  is  substantial,  and  a  good  portion  of  it 
can  be  correlated  with  temperature  fluctuations  in  the  vicinity  of  the  WVLA.  These 
were  logged  in  deep  water  thermistor  records  from  the  WVLA  and  three  dedicated 
thermistor  strings  moored  in  a  cluster  4.5  km  shoreward.  The  solid  curve  in  Figure  6-3 
is  obtained  from  the  average  change  in  sound  speed  seen  at  seven  of  these  thermistors 
applied  over  a  distance  of  8000  m.  The  seven  records  used  (three  from  the  WVLA,  two 
from  one  dedicated  string,  and  one  each  from  the  other  two)  were  in  the  depth  band 
of  interest  for  mode  1  (40-60  m)  and  judged  to  be  reliable.  The  foot  of  the  shelfbreak 
front,  which  has  been  observed  (in  CTD  profiles)  to  advect  across  the  position  of  the 
WVLA,  seems  a  likely  source  of  these  temperature  fluctuations. 

The  actual  range  between  the  thermistors  is  only  4500  m,  but  the  effects  may  be 
reasonably  extended  further,  and  the  choice  of  8000  m  provides  a  level  of  predicted 
variance  close  to  that  observed  in  the  acoustic  data.    However,  the  PAM1  wander 
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Mode  1  Peaks  and  Predicted  Current  Signals 
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Figure  6-2:  Predicted  current  effects  on  PAM1  arrival  time  plotted  with  the  distribu- 
tion of  peak  arrival  times  shown  in  Figure  6-1.  The  solid  curve  is  based  on  the  mean 
current  between  39  and  59  m  depth  and  the  dotted  curve  is  based  on  the  average 
current  between  11  and  67  m  depth. 
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Current-Corrected  Mode  1  and  Predicted  Temperature  Signal 
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Figure  6-3:  Predicted  temperature  fluctuation  effects  on  PAM1  arrival  time  (solid 
curve)  plotted  with  the  distribution  of  "current-corrected"  peak  arrival  times.  The 
solid  curve  is  zero  mean  (by  definition)  and  offset  from  the  data  for  convenience. 
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Mode  1  Peak  Arrival  Bias 
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Figure  6-4:  Mode  1  peak  arrival  time  bias  histograms  for  the  data  shown  in  Figure  6-1 
and  the  model  distribution  shown  Figure  5-15.  The  bias  measurements  are  relative 
to  estimated  PAM1  arrival  times  (i.e.  the  leading  edge)  for  the  data  and  calculated 
PAM1  arrival  times  for  the  model. 

occurring  on  yearday  214  is  a  clear  indicator  that  large  temperature  fluctuations 
are  not  restricted  to  the  area  around  the  WVLA,  and  it  is  likely  that  changes  in 
temperature  along  the  entire  32  km  path  contribute  at  various  times  in  varying  degrees 
to  the  mode  1  wander  signal  at  the  WVLA. 

Regardless  of  its  source,  the  wander  signal  evident  in  Figure  6-1  can  be  estimated 
and  subtracted  from  the  distribution  of  peak  arrival  times.  A  histogram  of  the  re- 
sulting peak  arrival  "bias"  levels  is  shown  in  the  upper  frame  of  Figure  6-4.  A  mode 
1  peak  arrival  time  bias  histogram  obtained  from  the  "propagated  thermistor  model" 
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in  Chapter  5  is  shown  in  the  lower  frame  of  the  same  figure.  The  histograms  are  quite 
similar,  even  though  the  "Data"  is  based  on  over  32,000  arrivals  and  the  "Model" 
uses  only  816.  Not  surprisingly,  the  "Model"  histogram  is  not  quite  as  smooth;  it 
also  shows  a  negatively  biased  pocket  that  stems  from  the  period  when  the  modeled 
path  averaged  group  speed  for  mode  2  was  higher  than  mode  1  (see  Figure  5-15). 

6.1.2      Arrival  Spread 

The  range  and  variability  of  positive  bias  seen  in  the  Figure  6-1  peak  arrival  distribu- 
tion is  indicative  of  significant  time  spreading  in  the  mode  1  arrival  patterns.  Spread 
is  perhaps  the  simplest  quantifying  measure  of  signal  distortion,  and  as  such,  it  is  a 
clear  indicator  of  the  mode  scattering  produced  by  the  waveguide. 

Spread  is  conceptually  simple,  yet,  as  discussed  in  Chapter  5,  the  choice  of  which 
mathematical  definition  is  "best"  to  use  is  somewhat  subjective.  The  four  measures 
that  will  be  considered  here  are:  (1)  the  IQR  of  the  arrival  envelope;  (2)  the  standard 
deviation  of  the  arrival  envelope;  (3)  the  absolute  value  of  the  peak  minus  the  mean 
arrival  time;  and  (4)  the  standard  deviation  of  the  peak  times  in  a  22  sequence 
transmission.  This  list  is  obviously  not  exhaustive  but  is  a  good  selection  of  the 
possibilities. 

Using  method  one,  the  IQR,  several  different  bin  averages  of  the  WVLA  mode  1 
arrival  spread  are  shown  in  Figure  6-5.  The  upper  frame  is  a  plot  of  the  average  spread 
for  each  transmission,  while  lower  frame  includes  five  different  bin  averages  ranging 
from  1.4  to  4.1  hours  for  comparison.  The  upper  frame  of  Figure  6-5  reveals  significant 
variability  over  a  time  scales  ranging  from  several  minutes  up  to  semidiurnal  tidal 
oscillations.  The  tidal  oscillations,  among  some  other  low  frequency  signals,  are  also 
evident  in  the  longer  averages  shown  in  the  lower  frame. 

The  second  method,  standard  deviation,  is  very  similar  to  the  first;  the  main 
difference  being  in  the  sensitivity  to  outlying  energy.  The  third  method,  absolute 
value  of  peak  minus  mean,  strictly  speaking,  is  more  a  measure  of  "nongaussianess" 
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Figure  6-5:  Mode  1  arrival  spread,  as  measured  by  the  IQR  of  the  arrival  envelope 
(unscattered  IQR  is  about  7.2  ms).  The  upper  frame  is  plot  of  the  spread  when 
averaged  over  each  transmission  (22  sequences  per  transmission).  The  lower  frame 
shows  bin  averaged  spreads  for  five  bin  widths  ranging  from  1.4  to  4.1  hrs. 
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than  a  measure  of  spread,  but  as  we  will  see,  these  usually  go  hand-in-hand.  Whereas 
the  first  three  methods  require  evaluation  of  the  arrival  envelope,  the  fourth  method, 
standard  deviation  of  peak  times,  only  requires  the  peak  arrival  times  to  be  recorded. 
The  tradeoff  for  this  is  that  only  one  measurement  can  be  obtained  per  transmission, 
as  opposed  to  22  per  transmission  for  the  other  three. 

Applying  each  of  these  methods  to  the  mode  1  arrival  records  from  SWARM 
produces  varying  degrees  of  similarity  between  methods.  Using  1.4  hour  bin  averages 
(see  Figure  6-5)  as  benchmarks,  the  following  matrix  of  zero-lag  cross-correlation 
values  can  be  constructed: 

P  1.00      .97      .85      .57 

1.00      .82      .60 

1.00      .29 

1.00 

As  expected,  the  correlation  is  the  strongest  between  methods  one  and  two,  and  the 
other  combinations  show  limited  similarity.  The  IQR  will  remain  the  primary  method 
by  which  spread  is  measured. 


Tidal  Signals  in  Spread  Spectra 

Figure  6-6  contains  the  spectra  of  the  spread  time  histories  in  the  lower  frame  of  Figure 
6-5.  The  spectra  confirm  the  existence  of  significant  variability  near  the  semidiurnal 
tidal  frequency.  This  tidal  variability  is  consistent  with  our  understanding  of  the  most 
probable  source  of  the  spread  in  the  arrivals:  internal  wave  solitons. 

The  relationship  between  tides  and  soliton  formation  appears  to  be  strong,  but 
it  is  not  strict  and  exclusive,  and  it  is  certainly  not  linear;  the  numbers  and  sizes  of 
solitons,  as  well  as  the  extent  of  solibores  (internal  tide  depressions)  they  are  often 
associated  with,  can  change  dramatically  from  one  tidal  cycle  to  the  next.  This, 
coupled  with  multiple  soliton-generation  sites  and  the  nonlinear  relationships  between 
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Figure  6-6:   Power  Spectral  Density  of  the  bin  averaged  spreads  shown  in  the  lower 
frame  of  Figure  6-5. 
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Mode  1  IQR  vs  WVLA  Temps 
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Figure  6-7:  WVLA  temperature  as  function  of  time  and  depth  overlaid  by  the  1.4  hr 
bin  averaged  IQR  of  mode  1  arrivals  at  the  WVLA.  The  temperatures  in  this  figure 
range  from  7.2  to  26.5  degrees  C. 


spread  and  soliton  density,  accounts  for  the  broadband  nature  of  the  low  frequency 
spread  spectrum. 

Figure  6-7  illustrates  both  the  variability  of  the  soliton  packets  from  cycle  to  cycle 
and  the  variability  in  the  degree  to  which  mean  spread  levels  are  controlled  by  the 
density  of  solitons  in  the  near  receiver  half  of  the  waveguide. 

This  figure  is  similar  to  its  simulation  counterpart,  the  upper  frame  of  Figure  5-17, 
but  it  differs  in  several  important  considerations:  first,  the  mean  spread  levels  are  from 
actual  arrivals,  not  simulations;  second,  the  mean  spread  levels  are  calculated  using 
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upwards  of  100  sequences  per  hour  (depending  on  the  number  of  gaps  in  the  data), 
instead  of  just  10  per  hour  in  the  simulations;  third,  the  bin  width  for  averaging  is  1.4 
hrs,  vice  30  minutes;  and  finally,  the  temperature  profiles  in  the  area  670  minutes  to 
the  left  of  the  spread  reading  in  question  are  only  an  indication  of  what  the  waveguide 
may  have  looked  like  at  the  time  of  the  transmission. 

Both  figures  show  some  degree  of  correlation  between  spread  levels  and  soliton 
activity  near  the  receiver,  but  it  is  difficult  to  quantify.  If  the  standard  deviation  of 
the  temperature  at  22.5  m  depth  over  the  previous  4.3  hrs  is  again  used  as  a  crude 
measure  of  soliton  activity  near  the  receiver,  then  we  can  see  a  clearer  picture  of  the 
correlation,  at  least  for  the  first  half  of  the  record,  as  shown  in  Figure  6-8.  This  figure 
also  shows  that  there  is  some  correlation  between  the  mean  spread  levels  from  the 
simulations  and  the  mean  spread  levels  found  in  the  data.  Obviously  the  similarities 
could  be  coincidental,  as  the  correlation  is  not  strong  and  consistent,  but  given  the 
nonlinear  spread  response  to  soliton  density  and  the  unavoidable  differences  between 
the  modeled  and  actual  waveguides,  one  cannot  expect  extremely  close  agreement. 
The  mere  fact  that  the  behaviors  show  some  limited  agreement  may  not  be  solid 
proof  of  the  veracity  of  the  model,  but  neither  are  the  differences  evidence  against 
that  proposition! 

6.1.3      Peak  Arrival  Intensity 

Mode  1  arrivals  with  large  amounts  of  distortion,  gained  that  distortion  through 
mode-coupling  interactions.  Greater  amounts  of  spread  generally  indicate  greater 
percentages  of  subarrivals  having  spent  large  amounts  of  time  in  higher  modes.  The 
act  of  time-spreading  alone  will  reduce  the  peak  intensity  of  a  mode  1  arrival,  and 
this,  in  combination  with  travel  in  higher  modes  causing  more  attenuation,  should 
create  a  signal  in  the  mean  peak  intensity  of  the  arrivals  that  is  inversely  proportional 
to  mean  spread  in  the  arrivals.  A  comparison  between  average  levels  of  spread  and 
relative  intensity  loss  of  mode  1  arrivals  is  made  in  Figure  6-9.  The  degree  to  which 
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Mode  1  IQR  (ms)  and  T-String  STD  (centidegrees) 
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Figure  6-8:  A  comparison  of  1  hour  bin  averaged  IQR  levels  for  the  SWARM  data 
(solid  line),  the  corresponding  propagated  thermistor  string  model  (dotted  line),  and 
the  Standard  Deviation  of  the  Temperature  at  22.5  m  depth  over  the  previous  4.3  hrs 
(dashed  line).  The  standard  deviation  of  the  temperature  is  in  tenths  of  a  degree  C 
to  facilitate  comparison  with  the  IQR  levels  in  milliseconds. 
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1 .4  hr  Bin  Averages  of  Mode  1  Peak  Intensity 
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Figure  6-9:  A  comparison  of  1.4  hr  averaged  peak  intensity  fluctuations  and  the 
corresponding  1.4  hr  averaged  IQR  levels  for  the  SWARM  data.  The  peak  intensity 
fluctuations  are  plotted  in  terms  of  relative  loss  to  provide  a  signal  that  is  in  phase 
with  the  spread  fluctuations. 
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these  correlate  (the  peak  cross-correlation  is  .29  at  zero  lag)  is  an  indication  of  the 
strength  of  the  link  between  mode  1  spread  and  peak  intensity.  Much  of  the  near-tidal 
variation  in  both  statistics  appears  to  be  in  phase,  but  the  majority  of  the  remaining 
variance  does  not  appear  to  correlate.  While  scattering  is  probably  the  dominant 
influence  on  spread  and  a  major  source  of  intensity  fluctuations,  Figure  6-9  shows  the 
relationships  to  be  nonlinear  and/or  the  influence  of  other  factors  to  be  substantial. 
The  mode  excitation  at  the  source,  for  example,  can  have  variability  over  a  wide 
range  of  time  scales.  The  effect  of  this  on  peak  intensity  is  somewhat  predictable, 
but  the  influence  on  spread  is  situational.  Conversely,  as  observed  in  the  Chapter  5 
simulations,  shifts  in  relative  group  speed  can  influence  spread  without  necessarily 
affecting  peak  intensity  one  way  or  the  other. 

6.1.4     Higher  Mode  Statistics 

Peak  Height 

An  interesting  comparison  can  be  made  using  the  mean  peak-height  statistics,  shown 
in  Figure  6-10,  for  the  first  four  modes.  The  zero-lag  cross-correlation  coefficients 
between  pairs  of  these  peak  height  statistics  are  compared  with  the  same  statistics 
for  simulated  data  in  Table  6-1.  The  peaks  fluctuate  in  height  coherently  to  a  great 
degree;  roughly  half  the  variance  is  coherent  between  modes.  The  simulated  data 
is  the  same  shown  in  Figures  5-14  and  5-15,  and  the  cross-talk  contamination  is 
simulated  with  incoherent  mixing  based  on  the  "Average"  filter  performance  matrix 
shown  in  section  3.3.2.  The  increases  in  the  statistical  correlations  caused  by  this 
mixing  are  in  some  cases  significant,  and,  no  doubt,  some  part  of  the  correlation  seen 
in  the  SWARM  data  is  also  from  cross-talk.  The  zero-lag  cross  correlations  between 
the  "unmixed"  simulated  data  and  the  "cross-talk"  simulated  data  are  0.98,  0.98, 
0.92,  and  0.96  for  modes  1  to  4  respectively. 

Figure  6-10  shows  that  much  of  the  similarity  in  the  SWARM  data  is  from  what 
appear  to  be  coherent  tidal  fluctuations  in  mode  excitation  and/or  attenuation.  The 
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Figure  6-10:  A  comparison  of  1.4  hr  averaged  peak-height  fluctuations  for  the  first 
four  modes  of  the  SWARM  data.  Mode  1  is  solid,  mode  2  is  dot-dashed,  mode  3  is 
dashed,  and  mode  4  is  dotted. 
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Bin  Avg 

Peak  Height 

1 

Zero-Lag  C 

ross-  Correlations 

Mode 
Pair 

SWARM 
Data 

Simulated  Data 

Unmixed 

Cross- Talk 

1:2 

.61 

.48 

.50 

1:3 

.52 

.26 

.34 

1:4 

.44 

.24 

.29 

2:3 

.74 

.65 

.65 

2:4 

.63 

.39 

.49 

3:4 

.83 

.57 

.86 

Table  6.1:  Zero-lag  cross  correlations  of  1.4  hour  bin  averaged  peak  heights  for  the 
SWARM  data  and  the  simulated  data  of  Chapter  5.  The  unmixed  data  (modes  1  and 
2)  are  shown  in  Figures  5-14  and  5-15;  the  cross-talk  contaminated  is  from  the  same 
source,  with  incoherent  mixing  based  on  the  "Average"  filter  performance  matrix 
shown  in  section  3.3.2. 

tidal  signal  becomes  more  apparent  and  shifts  in  phase  as  the  mode  number  increases. 


Peak  Height  and  Spread  Correlations 

Experience  with  the  simulated  transmissions  in  Chapter  5  and  the  existence  of  tidal 
fluctuations  in  the  peak-height  statistics  suggests  a  probable  correlation  between 
spread  and  peak  height  statistics.  Figure  6-11  shows  the  result  of  cross-correlations 
between  mode  1  spread  and  mode  1-4  peak  heights.  Confirming  our  visual  analysis  of 
Figure  6-9,  we  see  that  peak  height  fluctuations  for  mode  1  have  a  negative  a  zero-lag 
correlation  with  mode  1  spread  statistics.  This  is  in  agreement  with  simulations  using 
the  "propagated  thermistor  string"  (see  Figure  5-18).  The  correlation  becomes  pos- 
itive and  develops  a  tidal  period  as  the  mode  number  is  increased;  this  is  consistent 
with  the  concept  of  scattering  being  an  energy  sink  for  mode  1  and  an  energy  source 
for  higher  modes. 

The  zero-lag  cross-correlations  between  spread  levels  for  modes  1-4  and  peak 
heights  for  modes  1-4,  tabulated  in  Table  6-2,  are  similar  in  magnitude  and  trend 
to  those  observed  with  the  Chapter  5  simulations.  The  spread  measures  of  this  table 
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Figure  6-11:  Cross-correlations  of  1.4  hr  averaged  peak  height  fluctuations  for  mode 
1  (solid),  2  (dashed),  3  (dot-dashed),  and  4  (dotted)  with  the  corresponding  1.4  hr 
averaged  mode  1  IQR  (spread)  levels  for  the  SWARM  data. 


Zero-Lag  Cross- Correlations 

Spread 
Mode 

Peak  Height  Mode 

1 

2 

3 

4 

1 

-.53 

-.31 

-.21 

-.08 

2 

-.20 

-.48 

-.36 

-.26 

3 

-.11 

-.23 

-.35 

-.22 

4 

-.05 

-.27 

-.40 

-.43 

Table  6.2:  Zero-lag  cross  correlations  between  1.4  hour  bin  averaged  mode  1-4  peak 
heights  and  mode  1-4  mean  spread  levels.  The  spread  measurement  employed  for  these 
cross-correlations  is  the  standard  deviation  of  peak  arrival  times  over  a  transmission. 
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are  computed  by  method  number  four  (as  labeled  above),  the  standard  deviation  of 
the  peak  arrival  times  taken  over  a  transmission.  This  method,  chosen  primarily  for 
its  speed,  produced  a  greater  magnitude  correlation  between  mode  1  spread  and  peak 
height  than  did  the  calculation  using  the  IQR,  shown  in  Figure  6-11. 

Peak  Arrival  Time 

Peak  arrival  scattergrams  for  the  higher  modes  look  similar  to  the  mode  1  scattergram 
shown  in  Figure  6-1.  The  scattergram  data  for  the  first  four  modes  are  summarized, 
by  way  of  their  leading  edges  and  bin  averaged  mean  arrival  times,  in  the  plots  shown 
in  Figure  6-12.  Although  somewhat  busy,  this  figure  shows  still  another  statistic  with 
a  strong  tidal  influence.  The  mean  mode  4  arrival  time  has  a  definite  tidal  influence 
that  is  over  and  above  the  tidal  modulation  of  the  PAM1  arrival  (i.e.  the  bumps 
in  the  latest  dotted  line  have  larger  amplitudes  than  the  bumps  in  the  earliest  solid 
line).  These  mean  mode  4  arrivals  are  also  negatively  biased  relative  to  the  predicted 
PAM4  arrival,  because  they  are  generally  within  80  ms  of  the  mode  1  leading  edge 
(80  ms  being  the  average  difference  between  computer  simulated  PAM1  and  PAM4 
arrivals). 

6.2     224  Hz  Arrival  Analysis 

The  224  Hz  tomography  source,  moored  about  1000  m  farther  up  the  shelf,  provides 
a  second  source  of  acoustic  arrivals  at  the  WVLA.  The  narrower  bandwidth  of  the 
224  Hz  source  yields  a  pulse  width  of  about  60  ms  versus  the  10  ms  pulse  width  of 
the  400  Hz  signal. 

6.2.1      Peak  Arrival  Time 

The  wider  pulse  width  of  the  224  Hz  source  no  doubt  contributes  to  the  somewhat 
fuzzy  distribution  of  224  Hz  mode  1  peak  arrival  times  seen  in  Figure  6-13.    It  is 
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Figure  6-12:  Mode  1  (solid),  mode  2  (dashed),  mode  3  (dot- dashed),  and  Mode  4 
(dotted)  peak  arrival  time  leading  edge  (connects  earliest  peak  arrivals  in  1.1  hr  bins) 
and  1.4  hr  bin  average  mean  peak  arrival  times.  The  leading  edges  for  modes  2-4  are 
bracketed  between  the  two  mode  1  solid  lines,  and  the  mean  arrivals  for  modes  2-4 
are  positioned  above  the  upper  solid  line. 
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Figure  6-13:  Distribution  of  SWARM  224  Hz  Mode  1  Peak  arrival  times.  The  zero 
reference  is  at  an  arrival  time  of  22.453  seconds.  The  plotted  transmissions  occurred 
every  five  minutes.  Only  8  sequences  (every  4th  from  1  to  29)  are  shown  of  the  29 
available.  An  upshifted  leading  edge  of  the  400  Hz  peak  arrival  distribution  is  also 
shown  for  comparison.  The  later  zero  reference  time,  relative  to  Figure  6-1,  stems 
from  the  difference  in  source  mooring  locations. 
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Daily  Average 

Range  Independent  Travel  Times 

in  Excess  of  21  Seconds  (ms) 

Mode 

Freq 

Yearday  1995 

211 

212 

213 

214 

215 

216 

1 
1 

224 
400 

522 

518 

576 
574 

586 
579 

589 

584 

550 
545 

483 
476 

Difference 

4 

2 

7 

5 

5 

7 

2 
2 

224 
400 

553 
533 

624 
584 

632 
603 

625 
608 

593 
562 

539 
504 

Difference 

20 

40 

29 

17 

29 

35 

3 
3 

224 
400 

596 
556 

677 
636 

668 
640 

672 
622 

637 
601 

597 
542 

Difference 

40 

41 

28 

50 

36 

55 

Table  6.3:  Hypothetical  range  independent  travel  times  corresponding  to  daily  av- 
eraged SSPs  at  the  WVLA.  The  travel  time  differences  between  224  Hz  and  400  Hz 
modes  are  also  listed. 


difficult  to  estimate  just  what  portion  of  the  variability  in  the  224  Hz  peak  arrival 
time  should  be  attributed  to  the  broader  pulse,  but  using  the  shape  of  the  400  Hz 
leading  edge  as  a  guide,  an  uncertainty  of  around  ±5  ms  appears  to  be  appropriate. 
A  similar  qualitative  analysis  of  Figure  6-2  leads  us  to  believe  the  uncertaintity  for 
the  400  Hz  peak  arrivals  is  less  than  about  ±2  ms. 

Even  with  the  added  uncertainty,  tidal  modulation  in  the  224  Hz  PAM1  arrival 
is  still  apparent,  and  the  separation  between  the  224  Hz  and  400  Hz  PAM1  arrival 
times  can  be  seen  to  increase  by  10  to  15  ms  over  the  course  of  the  week.  One  possible 
source  of  this  divergence  is  perhaps  the  frequency  dependence  of  the  temporal  changes 
in  modal  group  speed.  In  support  of  that  hypothesis,  Table  6-3  lists  hypothetical 
range  independent  travel  times  corresponding  to  daily  averaged  SSPs  at  the  WVLA. 
Unfortunately,  the  table  shows  that,  at  least  in  the  vicinity  of  the  WVLA,  temporal 
changes  in  mode  1  group  speed  differences  do  not  appear  to  be  large  enough  to  account 
for  the  divergence  seen  in  Figure  6-13.  The  trend  is  basically  correct,  but  the  range 
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in  values  is  only  5  ms.  There  may  be  greater  (unmeasured)  temporal  fluctuations  at 
other  locations  along  the  acoustic  path  that  might  account  for  the  divergence,  but  an 
equally  likely  explanation  might  involve  a  fading  of  the  PAM1  arrival  as  the  dominant 
source  of  the  leading  edge.  The  edge  is  quite  fuzzy  to  start  with,  and  the  loss  of  a 
strong  PAM1  signal  is  not  unprecedented;  PAM1  signal  fading  is  clearly  evident  in 
Figure  5-15.  A  gradual  fading  of  224  Hz  PAM1  dominance  would  account  for  a  trend 
towards  increased  separation  between  the  400  Hz  leading  edge  and  the  224  Hz  leading 
edge. 

If  one  assumes  the  224  Hz  leading  edge  as  never  being  exclusively  controlled  by  the 
PAM1  arrival  time,  but  rather  being  influenced  by  the  PAM1,  PAM2,  and  PAM3  ar- 
rival times  in  combination,  another  explanation  of  the  Figure  6-13  divergence  emerges. 
The  PAM2  and  PAM3  arrival  times  have  more  frequency  dependence  (224  Hz  vs. 
400  Hz)  than  the  PAM1  arrival  time,  and  an  average  of  the  three  at  224  Hz  pro- 
duces the  sizeable  ramp  shown  in  Figure  6-14  (when  the  400  Hz  PAM1  arrival  time 
is  subtracted).  This  ramp  has  the  correct  trend  and  shows  an  average  increase  in 
separation  of  about  15  ms  over  six  days;  but,  be  cautioned,  it  does  not  represent  any 
sort  of  physical  model;  it  is  simply  a  reasonable  combination  of  statistics  with  proper 
physics  in  the  problem  that  can  perhaps  support  the  divergence  observed  in  Figure 
6-13. 

The  rough  tidal  signal  observed  in  the  224  Hz  arrivals  of  Figure  6-13  is  not  incon- 
sistent with  the  idea  of  combined  influence  from  all  three  modes  at  the  leading  edge. 
The  leading  edges  of  the  higher  mode  400  Hz  arrivals  (see  Figure  6-12)  are  certainly 
influenced  by  the  higher  mode  pseudo-adiabatic  arrivals,  yet  the  tidal  signal  is  still 
evident;  the  same  holds  for  the  mean  400  Hz  arrival  times  of  Figure  6-12. 

The  corresponding  tidal  signal  fades  more  rapidly  in  the  224  Hz  arrivals.  Only 
very  rough  tidal  signals  are  evident  in  any  of  the  time  series  of  average  224  Hz  peak 
arrival  times  and  mean  peak  heights  plotted  in  Figure  6-15,  but  as  the  corresponding 
power  spectral  densities  (Figure  6-16)  illustrate,  M2  tidal  influences  are  still  present 
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Figure  6-14:  Average  arrival  time  difference  between  the  224  Hz  PAM1,  PAM2,  and 
PAM3  arrivals  and  the  400  Hz  PAM1  arrival.  The  arrival  times  are  based  on  12  hour 
average  SSPs  at  the  WVLA  employed  in  a  range  independent  fashion  over  a  distance 
of  32  km. 
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Figure  6-15:  224  Hz  1.4  hr  bin  averaged  peak  arrival  times  and  peak  heights.  Mode 
1  is  solid,  mode  2  is  dashed,  and  mode  3  is  dotted. 

in  all  six  curves.  Subtidal  oscillations,  on  the  other  hand,  which  are  evident  in  the 
Figure  6-12  time  series  of  400  Hz  mean  peak  arrival  time,  can  still  be  clearly  observed 
in  the  time  series  of  the  224  Hz  records. 

6.2.2      Spread  Statistics 

The  time  histories  for  224  Hz  and  400  Hz  bin  average  mode  1  spread  (IQR)  statistics 
are  compared  in  Figure  6-17.  Both  records  exhibit  fairly  strong  tidal  modulation 
of  spread  during  the  first  part  of  the  recorded  period.  The  latter  portions  of  both 
records  have  limited  tidal  signatures,  but  the  modulation  that  does  exist  is  somewhat 
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Figure  6-16:  Power  spectral  densities  for  the  224  Hz  1.4  hr  bin  averaged  peak  arrival 
time  and  peak  height  time  series  of  Figure  6-15.  Mode  1  is  solid,  mode  2  is  dashed, 
and  mode  3  is  dotted. 
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Figure  6-17:  Comparison  of  224  Hz  and  400  Hz  1.4  hr  bin  averaged  IQR  (spread). 
The  IQR  of  an  unscattered  224  Hz  pulse  would  be  about  37  ms  (7  ms  for  a  400  Hz 
pulse). 


141 


similar,  and  the  mean  levels  of  spread  exhibited  are  both  about  23  ms  above  their 
respective  unscattered  levels  of  7  and  37  ms.  The  zero-lag  cross  correlation  between 
the  two  entire  records  is  0.72. 

6.2.3      Utility  of  Further  224  Hz  Analysis 

The  analysis  of  224  Hz  arrivals  has  to  this  point  been  somewhat  limited  compared 
to  the  analysis  of  the  400  Hz  arrivals,  but,  given  the  reduced  resolution  of  the  224 
Hz  signals  and  the  similarity  they  have  exhibited  with  corresponding  400  Hz  signals, 
further  investigation  in  this  area  is  less  enticing.  This  is  not  to  say  that  good  infor- 
mation, possibly  of  better  quality  than  exists  at  present,  cannot  be  extracted  from 
the  224  Hz  arrivals.  Rather,  the  information  obtained  would  be,  for  the  most  part, 
very  similar,  but  of  a  poorer  quality,  to  that  obtained  from  the  400  Hz  signals. 

Figure  4-7  illustrates  the  variation  in  mode  shapes  from  445  to  355  Hz  (90  Hz 
spread).  The  differences  between  400  and  224  Hz,  a  176  Hz  spread,  are  similar,  but 
as  one  might  expect,  larger  by  a  factor  of  about  two.  These  differences  are  notable  but 
not  dramatic,  so  there  is  little  to  be  gained  in  the  way  of  sampling  a  different  portion 
of  the  water  column.  The  224  Hz  and  400  Hz  coupling  characteristics  are  predictably 
similar,  as  evidenced  by  Figure  6-17,  so  if  additional  effort  is  to  be  expended  on  data 
analysis,  the  majority  should  be  spent  on  additional  analysis  of  the  clearer  400  Hz 
signals. 
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Chapter  7 


Conclusions 


7.1      Conclusions 

The  conclusions  reached  in  this  investigation,  with  some  overlap,  can  be  placed  in  two 
broad  categories:  modeling  and  data  analysis.  The  overlap  was  particularly  important 
in  estimating  the  level  of  bottom  attenuation  and  verifying  the  importance  of  soliton 
location  to  the  level  of  signal  distortion.  A  detailed  listing  of  conclusions  follows. 

7.1.1      Modeling 

The  following  conclusions  can  be  drawn  from  the  modeling  efforts  of  Chapter  4  and 
Chapter  5: 

1.  The  average  400  Hz  compressional  wave  attenuation  in  the  bottom,  along  the 
SWARM  acoustic  path,  is  in  the  neighborhood  of  0.02  to  0.03  dB/m.  A  sed- 
imentary compressional  wave  attenuation  level  of  0.092  dB/m  can  be  inferred 
from  chirp  sonar  inversions  for  sound  speed  (linear  increase  from  1500  to  1800 
m/s  over  30  m)  and  density  (1.8  g/cm3)  in  sediment  near  the  WVLA.  Ex- 
tending these  bottom  parameters  in  a  range  independent  fashion  towards  the 
source  produces  almost  complete  attenuation  of  energy  above  mode  2  in  simu- 
lated transmissions,  but  the  actual  data  shows  significant  energy  transmission 
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up  through  and  including  mode  4.  An  attenuation  level  of  around  0.0115  dB/m 
proved,  through  numerous  trials,  to  be  appropriate  for  our  modeled  30  m  sub- 
bottom;  slightly  higher  levels  (0.02  to  0.03  dB/m)  would  be  more  appropriate 
if  observed  along  track  reductions  in  the  actual  SWARM  subbottom  thickness 
were  considered. 

2.  Mode  arrival  spread  and  bias  statistics  can  be  very  sensitive  to  the  specific 
locations  of  scattering  events,  but,  in  general,  the  more  soliton  activity  there  is 
in  the  portion  of  the  waveguide  closest  to  the  receiver,  the  greater  the  spread  and 
bias  will  be  in  the  mode  arrivals.  The  SIA  model  results  show  this  very  clearly, 
with  the  "propagated  thermistor  string"  and  actual  SWARM  data  providing 
good  support. 

3.  Under  most  shallow  water  circumstances  (i.e.  when  mode  1  has  the  highest 
group  speed),  the  PAM1  arrival  time  can  estimated  from  a  distribution  of  peak 
mode  1  arrival  times  by  finding  the  leading  edge  envelope  of  the  distribution. 

4.  Under  most  shallow  water  circumstances  (i.e.  when  significant  coupling  between 
modes  exists),  the  PAM2  and  higher  arrival  times  are  not  easily  discernible  from 
distributions  of  peak  arrival  times. 

5.  At  present,  the  predictive  modeling  of  mode  arrivals  from  soliton  filled  oceanic 
waveguides  is  limited  more  by  inadequate  sampling  of  the  environment  in  ques- 
tion than  it  is  by  the  limitations  of  the  physics  included  in  the  computer  mod- 
els. (It  would  be  extremely  difficult  to  obtain  enough  environmental  data  in  a 
SWARM  like  area  to  "ground  truth"  a  model,  but  temperature  strings  placed 
every  5  km  or  so  would  probably  be  enough  to  put  the  ball  back  in  the  "physics 
limited"  court.) 
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7.1.2      Data  Analysis 

The  following  conclusions  can  be  drawn  from  the  data  analysis  covered  in  Chapter  6: 

1.  Mode  1  has  the  highest  path-averaged  group  speed  for  most,  if  not  all,  of  the 
duration  of  the  experiment. 

2.  Large  fluctuations  in  the  path-averaged  mode  1  group  speed  occurred  over  tidal 
and  subtidal  periods  during  the  course  of  the  experiment. 

3.  A  significant  percentage  of  the  mode  1  group  speed  fluctuations  can  be  at- 
tributed to  warm  deep  water  advections  in  the  latter  5  to  10  km  of  the  waveg- 
uide. These  advections,  which  are  probably  associated  with  the  foot  of  a  warm 
(and  salty)  shelfbreak  front,  can  be  seen  in  deep  thermistor  records  near  the 
receiver. 

4.  Both  comparisons  with  model  simulations  and  correlations  with  physical  mea- 
surements point  to  solitons  in  the  waveguide  as  the  major  source  of  mode  cou- 
pling during  the  SWARM  experiment. 

5.  Soliton-filled  waveguides  will  produce  mode  arrival  distortions  with  fluctuation 
time  scales  ranging  from  minutes  to  days. 

6.  The  amount  of  mode  arrival  distortion  a  particular  waveguide  will  produce  can 
only  be  predicted  in  a  general  sense  (e.g.  solitons  near  the  receiver  increase  the 
probability  of  mode  1  spread).  Correlations  between  simple  measures  of  soliton 
activity  and  simple  measures  of  mode  arrival  distortion  are  weak,  but  can  be 
seen,  when  averaged  over  one  or  two  hours. 
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7.2      Proposed  Directions  of  Future  Investigation 

7.2.1  Modeling 

In  the  area  of  modeling,  the  following  ideas  should  be  explored: 

1.  Investigate  the  possibility  of  incorporating  depth  dependence  in  the  propagated 
thermistor  string  model  through  linear  perturbations  on  depth  independent  cou- 
pling matrix  coefficients  and  depth  independent  modal  wavenumber  functions. 

2.  Look  for  ways  of  increasing  computational  speed  through  the  reduction  of  re- 
peated calculations  of  modal  wavenumber  functions. 

7.2.2  Data  Analysis 

In  the  area  of  data  analysis: 

1.  The  analysis  of  the  WVLA  and  NRL  VLA  records  should  be  coordinated  to  the 
greatest  extent  possible.  The  goal  being  to  obtain  similar  sets  relevant  statistics 
from  each  instrument. 

2.  An  attempt  to  measure  temporal  variations  in  the  mode  arrival  decorrelation 
times  should  be  made.  The  effort  will  be  hampered  by  the  broken  nature  of  the 
sampling  scheme. 

7.2.3  Experiments 

An  attempt  should  be  made  to  conduct  a  follow-on  study.  Some  suggestions  for  that 
study  are  as  follows: 

1.  Conduct  it  in  a  region  where  there  are  fewer  soliton  crossing  patterns  than  were 
observed  at  the  SWARM  site.  The  ideal  being  a  location  that  produces  regular 
soliton  packets  with  a  substantial  quiescent  gap  between  them. 
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2.  A  linear  array  of  four  sources  at  15,  20,  25,  and  30  km  from  the  receiver  should 
be  used  to  study  range  dependence  effects. 

3.  Each  array,  as  well  as  the  source,  should  be  instrumented  with  at  least  six 
thermistors.  Two  other  strings,  placed  at  5  and  10  km  from  the  source,  would 
round  out  a  good  sampling  of  the  waveguide.  This  would  aid  both  modeling 
efforts  and  acoustic/environment  correlation  analysis. 

4.  Schedule  transmissions  with  two  sequences  from  each  source  every  minute,  with 
the  30  km  source  starting  on  the  minute,  the  25  km  at  15  seconds,  etc.  This 
would  provide  one  pulse  compressed  sequence  per  source  per  minute  for  analysis. 

7.3      Contributions 

The  significant  contributions  this  investigation  has  supplied  to  the  field,  again  with 
some  overlap,  can  be  placed  in  two  broad  categories:  modeling  and  data  analysis. 

7.3.1      Modeling 

The  following  contributions  are  drawn  from  the  modeling  efforts  of  Chapter  4  and 
Chapter  5: 

1.  A  simple  method,  termed  the  pulse  propagation  model  (PPM),  by  which  range 
dependent  broadband  acoustic  propagation  can  be  modeled  is  developed  and 
employed  to  evaluate  the  effects  of  solitons  on  acoustic  propagation.  This  PPM 
method,  as  shown  in  Appendix  B,  compares  favorably  to  synthesizing  broad- 
band pulses  from  KRAKEN  [12]  complex  pressure  fields,  but  it  is  orders  of 
magnitude  faster. 

2.  By  comparison  with  Marshall  Orr's  backscatter  record,  the  use  of  a  "propagated 
thermistor  string"  record  to  simulate  soliton  filled  waveguides  is  shown  to  be  a 
reasonable  first  order  approach. 
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3.  The  importance  of  soliton  location  to  a  waveguides  mode  scattering  potential 
is  clearly  demonstrated. 

7.3.2     Data  Analysis 

The  following  contributions  are  drawn  from  the  data  analysis  covered  in  Chapter  6: 

1.  The  first  detailed  analysis  of  actual  acoustic  propagation  in  a  soliton  filled 
waveguide  is  presented. 

2.  Experimental  evidence  supporting  the  concept  of  soliton  location  being  critical 
to  a  waveguides  mode  scattering  potential  is  presented. 

3.  Experimental  evidence  supporting  the  veracity  of  a  "propagated  thermistor 
string"  record  to  simulate  soliton  filled  waveguides  is  presented. 
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Appendix  A 


Array  Navigation 


A.l      Navigation  Set-Up 

Three  transponder  balls  (11.0,  11.5,  and  12.0  kHz)  were  placed  in  a  pattern  around 
the  WVLA  (Figure  A-l).  The  pattern  is  not  optimum,  as  the  deflections  we  are 
concerned  with  (in  the  direction  of  the  acoustic  path)  are  nearly  orthogonal  to  both 
the  11.0  and  11.5  transponder  paths.  A  10.5  kHz  pinger  at  the  base  of  the  array 
transmitted  every  four  minutes  and  triggered  the  transponders.  The  arrival  time  of 
each  frequency  at  the  top  hydrophone  of  the  WVLA  (channel  zero)  was  recorded. 
The  ranges  from  each  transponder  to  the  WVLA  were  varied  to  ensure  good  time 
separation  in  the  arrivals.  With  the  mean  and  ramping  removed  from  the  travel- 
times,  any  remaining  fluctuations  that  were  not  common  to  all  three  paths  would,  in 
theory,  be  due  to  motion  of  the  the  receiving  hydrophone. 


A. 2     Expected  Motions  of  the  Array 

The  main  mooring  consists  of  a  one  meter  diameter  flotation  sphere  supporting  a 
one  inch  diameter  combination  strength  member  and  16  channel  hydrophone  array. 
(Figure  A-2).    The  length  of  the  of  the  mooring  is  approximately  60  meters  and  it 
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Figure  A-l:  Relative  positions  of  transponder  balls  around  the  WHOI  Vertical  Line 
Array  (WVLA).  Transponders  were  interrogated  every  four  minutes  by  a  pinger  at 
the  base  of  the  array.  Round  trip  travel-time  to  the  top  hydrophone  of  the  array  was 
recorded. 
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Figure  A-2:  A  simple  form  drag  model  for  deflection  of  the  WVLA  where  the  drag 
force  on  the  cable  is  all  applied  at  the  midpoint.  The  20  cm/s  current  shown  would 
create  a  displacement  of  about  .5  m  at  the  top  hydrophone. 


has  net  a  buoyancy  of  3230  Newtons.  The  sum  of  the  added  mass  and  actual  mass 
of  the  mooring  is  approximately  equal  to  its  displacement,  so  the  undamped  natural 
frequency,  ta  =  (<?//) *,  is  about  3  or  4  cycles  per  minute.  However,  large  impulse- 
response  type  oscillations  are  not  expected,  because  the  damping  from  drag  is  quite 
significant.  It  is  proportional  the  square  of  the  flow  velocity,  with  the  equivalent  linear 
damping  reaching  criticality  (becoming  over-damped)  at  flows  above  about  1  cm/s. 

After  calculating  the  form  drag  of  the  mooring  sphere  and  cable  in  a  barotropic 
current  field,  the  tilt  of  the  WVLA  can  be  estimated  by  determining  the  equilibrium 
point  between  the  drag  forces  and  the  righting  moment  provided  by  the  spheres 


151 


buoyancy.  At  flows  around  20  cm/s  the  drag  on  the  sphere  can  be  approximated  by, 


Sphere  «  ^(^X^2)  =  6.6JV, 


(A.l) 


and  the  drag  on  the  cable  would  be  about, 


^cable  «  l-2(-^2)(.0254m)(60m)  =  39JV. 


(A.2) 


If  all  the  cable  drag  was  applied  at  the  midpoint  of  the  cable,  the  deflection  at 
that  point  would  be, 


Ax  =  (30m)sin(arctan(39/3230))  =  .36m, 


(A.3) 


and  the  drag  on  the  sphere  would  provide  an  additional, 


Ax  =  (60m)sin(arctan(6.6/3230))  =  .12m, 


(A.4) 


for  a  total  of  about  .5  meters  at  the  top  hydrophone,  which  is  about  57  m  above  the 
bottom. 

A  more  complicated  inelastic  cantenary  equation  solution  can  be  derived  by  fol- 
lowing Triantafyllou  [20],  replacing  his  vertical  force  "cable  weight  per  unit  length  in 
water"  with  our  horizontal  force  "drag  per  unit  length" (our  cable  is  roughly  neutrally 
buoyant).  The  displacement  equation  is 


x  = 


V_ 
Dr 


1  + 


H  +  DrL 


V 


\ 


1  + 


H  +  Dr{L  -  s) 
V 


(A.5) 


where  5  is  the  distance  along  the  cable  of  length  L  measured  from  the  bottom.  When 
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the  previously  computed  20  cm/s  force  quantities, 

V(sphere  vertical  force)  =  3230iV 

//(sphere  horizontal  force)  =  6.6./V  (A. 6) 

Dr(cable  drag)  =  .65N/m, 

are  inserted,  a  deflection  of  .48  m  is  again  obtained  for  the  60  m  point  on  the  cable. 

Assuming  continuous  laminar  flow  (the  sphere  probably  transitions  to  turbulence 
at  around  30  cm/s),  the  displacement  for  barotropic  flows  will  be  roughly  proportional 
to  the  square  of  the  current  speed,  giving  1  meter  at  30  cm/s;  2  meters  at  40  cm/s; 
and  3  meters  at  50  cm/s. 

Figure  A-3  shows  the  upper  and  lower  depth  bins  of  a  bottom-mounted  acoustic 
Doppler  current  meter  (ADCP)  moored  about  3  km  from  the  WVLA.  Most  of  the 
barotropic  currents  are  less  than  30  cm/s.  The  largest  currents  are  baroclinic  soli- 
ton  currents  which  have  a  westward  component  near  the  surface  and  an  eastward 
component  near  the  bottom.  The  effects  of  baroclinic  currents  on  mooring  motion 
are  hard  to  predict,  but  the  movement  at  the  top  hydrophone  is  certainly  less  than 
what  would  be  produced  by  a  barotropic  current  of  the  same  magnitude.  This  line 
of  thought  coupled  with  the  ADCP  data  predicts  mooring  motions  that  are  generally 
less  than  a  meter  in  amplitude  at  the  top  hydrophone  (equivalent  to  1  degree  of  tilt) 
with  occasional  excursions  that  might  approach  2  meters. 

A. 3     Navigation  Data  Analysis 
A. 3.1      Inversion  for  Mooring  Motion 

Initial  inspection  of  the  navigation  data  reveals  travel-time  fluctuations  ranging  from 
one  to  five  milliseconds  (Figure  A-4).  This  roughly  translates  to  deflections  in  hy- 
drophone position  of  from  1.5  to  7.5  meters.  A  fairly  clear  tidal  signal  is  apparent 
throughout  the  11.0  Khz  record;  the  tidal  signal  is  weaker  in  the  11.5  Khz  signal, 
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Figure  A-3:  The  first  (..)  and  fourteenth  (-)  depth  bins  of  a  16  channel  bottom- 
mounted  acoustic  Doppler  current  meter  (ADCP)  moored  about  3  km  from  the 
WVLA. 
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Figure  A-4:  Travel-time  fluctuation  data- vectors  associated  with  the  three  transpon- 
der ball  paths  used  to  navigate  the  WVLA. 


nonexistent  in  the  first  half  of  the  12.0  Khz  record,  but  very  strong  in  the  latter 
portions  of  the  12.0  Khz  record. 

The  travel-time  fluctuations  can  be  more  directly  converted  to  equivalent  hy- 
drophone displacement  using  a  least-squares  inversion  of  the  following  model, 


-l 

c 


cos  #!  sin  0\ 
cos  02  sin  02 
cos  03    sin  03 


Ax 

Ay 


Ah 
At2 
At3 


(A.7) 
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Figure  A-5:  Model  vectors  of  an  over-determined  least-squares  inversion  of  the  navi- 
gation data  for  north  and  east  displacement  of  the  hydrophone  array. 


where, 


and 


0n  =  relative  bearing  of  due  east  from  transponder  bearing, 

Ax  =  east  displacement, 

Ay  =  north  displacement, 

Atn  =  travel-time  fluctuation  at  transponder, 

c  =  1500m/s. 


(A-8) 


(A.9) 


Figure  A-5  illustrates  two  problems  with  the  inverted  displacements:  they  are 
much  larger  than  expected  from  the  ADCP  data  and  form  drag  calculations;  and  the 
source  of  the  high  frequency  oscillations  is  unclear. 

The  low  frequency  oscillations  in  the  latter  part  of  the  record  have  a  clear  tidal 
signature,  the  east- west  component  being  forced  almost  directly  by  the  12.0  Khz 
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transponder.  The  orientation  of  the  tidal  ellipse  in  the  area  is  319  deg  T  for  the 
major  axis  component  (~  14  cm/s)  and  231  deg  T  for  the  minor  axis  component  (~  9 
cm/s).  Given  this  tidal  ellipse  and  the  strong  12.0  Khz  tidal  signal  starting  on  day 
215,  we  would  expect  to  see  reasonably  strong  tidal  signals  from  the  11.0  and  11.5 
transponders,  but  we  don't. 

The  source  of  the  high  frequency  displacement  signal  is  even  more  troublesome. 
If  we  choose  to  believe  the  low  frequency  inverted  displacements  of  order  ±5  m  are 
caused  by  tidal  currents  (i.e.  gross  errors  in  hydrodynamic  calculations),  then  we 
must  also  accept  the  strong  possibility  that  the  higher  frequency  displacements  of 
order  ±2.5  m  are  real,  the  result  of  high  frequency  current  oscillations  (see  Figure 
A-3). 

Another  source  of  both  high  and  low  frequency  transponder  travel-time  fluctua- 
tions is  soundspeed  changes.  An  internal  tide  can  depress  the  thermocline  causing 
a  travel-time  change  that  is  linearly  proportional  to  the  affected  pathlength;  since 
the  distance  to  the  12.0  Khz  transponder  is  twice  as  far  as  the  the  11.0  and  11.5 
transponders,  our  previous  inversion  could  interpret  this  as  a  mooring  motion.  This 
does  not  appear  to  be  a  significant  part  of  the  tidal  signal  in  the  12.0  Khz  record.  If 
it  were,  we  would  expect  to  see  similar  fluctuations  at  about  half  the  amplitude  of 
12.0  fluctuations  in  the  11.0  and  11.5  records. 

The  level  of  influence  from  internal  tide  displacements  can  be  estimated  more 
directly  by  using  a  modified  form  of  the  previous  inversion  (no  longer  over-determined) 
that  solves  for  an  average  sound  speed  fluctuation  along  with  the  east  and  north 
hydrophone  displacements.  The  speed  deviation  in  this  case  is  assumed  to  apply  to 
the  entire  return  path  from  each  transponder  to  the  upper  hydrophone  of  the  WVLA. 
In  actuality  it  only  applies  to  a  portion  of  each  path,  but  for  long  waves,  the  affected 
pathlength  will  be  roughly  a  common  fraction  of  each  total  pathlength.  The  near 
bottom  paths  from  the  interrogator  to  the  transponders  are  assumed  to  have  a  stable 
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sound  speed.  This  inversion, 


-l 

c 


cos  B\  sin  $i  d\  /c 
cos  $2  sin  92  d2/c 
cos  63    sin  03    d3/c 


Ax 

Ati 

Ay 

= 

At2 

Ac 

At3 

(A.10) 


where, 


and 


Ac  =  path  average  change  in  sound-speed, 


dn  =  distance  from  WVLA  to  transponder, 


(A-11) 


(A.12) 


produces  smaller  amplitudes  (still  larger  than  expected)  for  the  low  frequency  mooring 
motions  (Figure  A-6),  and  the  added  model  vector  for  speed  deviation  shows  some 
similar  (and,  as  expected,  smaller  magnitude)  fluctuations  when  compared  to  a  record 
of  sound  speed  in  the  upper  water  column  (Figure  A-7).  The  sound  speed  record  is 
converted  from  a  temperature  record  at  a  depth  of  12.5  meters  on  the  WVLA.  This 
inversion  shows  no  diminution  in  high  frequency  mooring  motion;  we  would  not  expect 
high  frequency  soundspeed  deviations  to  fit  into  the  long  wave  assumption. 


A. 3. 2      Inversion  with  Temperature  Correction  Applied  to 
Data 

The  similarity  between  the  thermistor  record  and  the  average  sound  speed  model 
vector  opens  the  possibility  for  applying  a  sound  speed  correction  to  the  data  vectors 
prior  to  inversion  for  mooring  motion.  As  it  turns  out,  the  travel-time  fluctuation 
data  vectors  can  be  modified  slightly  by  applying  the  predicted  local  change  in  sound 
speed  from  the  thermistor  data  to  one  eighth  (a  reasonable  and  optimum  fraction  for 
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Figure  A-6:  Model  vectors  of  a  least-squares  inversion  of  the  navigation  data  for  north 
displacement,  east  displacement,  and  an  average  deviation  in  sound  speed  common 
to  all  three  transponder  return  path  integrals. 
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Figure  A-7:  A  least-squares  inversion  for  sound  speed  deviation  produces  fluctuations 
similar  those  observed  in  a  12.5  meter  depth  thermistor  record  that  has  been  converted 
to  sound  speed. 
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this  data  set)  of  the  return  path  for  all  three  transponders,  but  this  inversion, 


cos  0i  sin  61 
cos  $2  sin  62 
cos  03    sin  03 

is  not  a  significant  improvement  over 


cos  0\  sin  0i 
cos  02  sin  02 
cos  03    sin  03 


Ax 
Ay 


cAt,  +  A* 
cAt2  +  AAfi 
cA*3  +  ^ 


(A.13) 


Ax 
Ay 


cA*i 

cA*2 
cA*3 


(A.14) 


because  applying  the  sound  speed  corrections  to  the  data  vectors  only  slightly  de- 
creases their  variances  from  2.4,  2.3,  and  12.3  m2  to  2.2,  2.2,  and  12.0  m2. 

Utilizing  the  thermistor  sound  speed  deviation  over  larger  portions  of  the  return 
path  will  actually  increase  the  variance  of  the  data  vectors.  This  is  not  unexpected; 
the  return  paths  start  out  deep  (70  m)  where  the  effects  of  solitons  on  the  sound 
speed  are  not  so  pronounced  as  one  will  find  at  12.5  meters.  Also,  it  must  be  kept  in 
mind  that  the  thermistor  record  is  influenced  by  both  global  and  local  temperature 
fluctuations. 


A. 3. 3      Inversion  for  Sound  Speed  Only 

Internal  tides  have  20  to  30  km  wavelengths,  but  internal  wave  solitons  have  wave- 
lengths on  the  order  of  200  meters  and  could  be  an  important  source  of  the  high 
frequency  oscillations  seen  in  the  transponder  records.  An  over-determined  inversion 
looking  only  for  speed  fluctuations  can  be  constructed  to  solve  for  two  separate  sound 
speed  fluctuations.  The  11.0  and  11.5  paths  are  more  likely  than  any  other  pair  to 
have  speed  deviations  in  common.  The  predominant  direction  of  solitary  wave  prop- 
agation is  300°  by  compass  heading,  so  with  the  soliton  phase  fronts  perpendicular 
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to  the  propagation  direction,  the  11.0  and  11.5  kHz  locations  could  very  possibly  ex- 
perience the  simultaneous  passage  of  a  solitary  wave  (Figure  A-l).  An  inversion  that 
solves  for  a  sound  speed  fluctuation  common  to  the  11.0  and  11.5  kHz  transponder 
return  path  and  a  separate  fluctuation  for  the  12.0  kHz  return  path, 
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<2i     0 

Ac12 
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d2     0 

Ac3 

— 

At2 

0     d3 

A<3 

(A.15) 


produces  very  interesting  results.  Figure  A-8  shows  that  the  deviations  in  sound 
speed  from  this  inversion  are  again  similar  to  the  temperature  (converted  to  speed) 
changes  recorded  at  the  top  of  the  WVLA.  A  closer  comparison  between  inverted 
speed  deviations  for  the  11.0  and  11.5  kHz  transponders,  Aci2,  and  a  demeaned 
upper  thermistor  speed  record  shows  very  good  agreement  for  the  early  part  of  the 
record  (Figure  A-9). 

A  comparison  between  the  12.0  khz  travel-times  and  the  two-way  transponder 
range  divided  by  the  average  of  the  three  deepest  sound  speeds  on  the  WVLA  is  shown 
in  Figure  A-10.  While  Figure  A-9  indicates  some  sensitivity  for  shallow  temperature 
fluctuations  in  the  11.0/11.5  khz  records,  this  figure  shows  the  12.0  khz  to  be  very 
sensitive  to  deep  temperature  fluctuations.  The  size  and  path  dependence  of  these 
temperature  signals  make  any  inversion  for  mooring  motion  unreliable. 

A. 4      Conclusions 

The  following  conclusions  are  reached  in  regard  to  navigation  data  at  the  WVLA: 


1.  The  average  length  scale  of  temperature  fluctuations  in  the  vicinity  of  the 
WVLA  was  less  than  1  km  (i.e.  most  temperature  fluctuation  signals  were 
not  common  to  all  three  paths). 
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Figure  A-8:  A  least-squares  inversion  for  sound  speed  deviation  produces  fluctuations 
similar  those  observed  in  a  12.5  meter  depth  thermistor  record  that  has  been  converted 
to  sound  speed. 
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1 1 .0  and  1 1 .5  kHz  and  Upper  Thermistor  Record 
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Figure  A-9:  A  least-squares  inversion  for  sound  speed  deviations  yields  an  11.0/11.5 
kHz  path  sound  speed  fluctuation  (solid  line)  that  is  very  similar  those  observed  in  a 
12.5  meter  depth  thermistor  record  (dashed  line). 


164 


12  Khz  Travel-time  vs  Range/(Avg  3  Deep  Speeds) 
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Figure  A- 10:  12.0  kHz  path  travel-time  fluctuations  (upper  line)  are  very  similar 
those  calculated  by  taking  the  range  to  the  transponder  and  back  divided  the  average 
of  the  WVLA  soundspeeds  at  40,  50,  and  60  meters  depth. 


165 


2.  Mooring  motion  signals  in  the  navigation  data  were  too  small  to  be  detected 
above  the  background  of  temperature  fluctuations. 

3.  Hydrodynamic  estimates  limit  mooring  motions  to  less  than  one  degree  of  tilt 
for  most  periods,  with  occasional  excursions  approaching  two  degrees  during 
the  passage  of  large  solitons. 
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Appendix  B 

Comparison  of  the  Pulse 
Propagation  Model  with  Pulse 
Synthesis  using  KRAKEN 
Complex  Pressure  Fields 


The  pulse  propagation  model  (PPM)  described  in  Chapter  4  relies  on  the  KRAKEN[12] 
normal  mode  code  for  computation  of  mode  shapes,  wave  numbers,  and  attenuation 
levels.  The  acoustic  FIELD  program,  within  the  KRAKEN  model,  can  also  com- 
pute the  complex  pressure  field, C(r,  z,u)  =  a(r,z,u>)  +  ib(r,  z,u>),  at  a  number  of 
selected  ranges  and  depths  for  a  given  source  depth,  source  frequency,  and  acoustic 
environment. 

If  the  field  calculations  are  repeated  for  a  number  of  frequencies  (covering  the 
bandwidth  of  the  transmitted  signal),  the  resulting  complex  spectrum  of  the  field  at 
a  given  point,  when  multiplied  by  the  Fourier  components  of  the  source  pulse,  P{u>), 
can,  after  Fourier  synthesis,  yield  a  time  domain  pulse, 

P(r,z,t)  =  ±Jc{r,z,u)P(u)e**du.  (B.l) 
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Since  complex  pressure  fields  are  not  smooth  functions  of  frequency,  an  accurate 
solution  using  the  SWARM  400  Hz  source  requires  1201  KRAKEN  runs  for  each 
different  SSP  in  a  range  dependent  problem.  This  procedure,  which  we  will  refer  to 
as  KRAKEN  Field  Synthesis  (KFS),  computes  the  arrival  at  a  single  hydrophone;  it 
can  be  extended  to  predict  mode  arrivals  by  mode-filtering  the  predicted  arrivals  at 
vertical  array  of  hydrophones. 

Conversely,  the  mode  arrival  predictions  of  the  pulse  propagation  model  can  be 
synthesized  to  yield  the  arrival  at  a  given  single  hydrophone.  This  computation, 

P{r,z,t)  =  Y^Pn(r,t)(j)n(zs)<f>n(z),  (B.2) 

n 

where  the  pn(r,t)  are  the  mode  arrivals,  the  <f){z)  are  the  mode  shapes,  and  zs  is 
the  source  depth,  is  very  simple  and  provides  an  output  for  direct  comparison  with 
the  KFS  method.  Since  the  modal  wavenumber  spectra  it  uses  for  propagation  are 
fairly  smooth  functions  of  frequency,  the  PPM  method  requires  only  5  KRAKEN 
runs  for  each  SSP.  The  wavenumbers  and  attenuation  levels  for  each  modes  1201 
non-zero  frequency  bins  are  then  computed  through  cubic  spline  interpolation.  The 
5  KRAKEN  runs  with  interpolation  used  in  the  PPM  are  orders  of  magnitude  more 
efficient  than  the  1201  KRAKEN  runs  needed  in  KFS. 

Figure  B-l  shows  a  comparison  between  the  two  methods  for  simple  range  depen- 
dent case.  The  modeled  waveguide  uses  the  SSPs  shown  in  Figure  4-6,  with  coupling 
from  one  SSP  to  the  other  occurring  at  13,000  and  18,000  m.  The  agreement  between 
the  two  methods  is  very  good  at  all  ranges  for  this  case. 

A  more  complicated  scenario  (Figure  B-2)  yields  results  that  are  not  as  closely 
matched  but  still  quite  good,  even  at  the  20  Km  point.  There  are  several  potential 
contributors  to  the  deviation  between  the  PPM  and  KFS  results.  The  most  obvious 
is  the  narrowband  approximation  used  in  PPM  coupling  matrices.  The  mode  shapes, 
as  illustrated  in  Figure  4-7,  do  have  a  slight  dependence  on  frequency,  and  as  such, 
the  coupling  matrices  are  also  frequency  dependent.   This  dependence  is  accounted 
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Figure  B-l:  A  comparison  of  single  hydrophone  pulse  arrivals  computed  from  the 
PPM  (solid)  and  KFS  (dotted)  methods.  The  modeled  waveguide  uses  the  SWARM 
SSPs  shown  in  Figure  4-6,  with  coupling  from  one  SSP  to  the  other  occurring  at 
13,000  and  18,000  m.  The  source  depth  is  27  m  and  the  receiver  depth  at  all  five 
ranges  is  50m. 
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Figure  B-2:  An  eight  soliton  waveguide  (upper  panel)  and  two  corresponding  arrivals 
computed  from  the  PPM  (solid)  and  KFS  (dotted)  methods.  The  modeled  waveguide 
uses  the  SWARM  SSPs  shown  in  Figure  4-6  to  simulate  the  solitons.  The  source  (*) 
depth  is  27  m  and  the  receiver  (o)  depth  at  both  ranges  is  50m. 
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Figure  B-3:  Single  hydrophone  RIP  pulse  arrivals  computed  from  the  PPM  (solid) 
and  KFS  (dotted)  methods  using  various  source  depths  with  the  receiver  depth  fixed 
at  50m. 


for  in  the  KFS  method,  but  the  PPM  method  assumes  the  400  Hz  coupling  matrices 
apply  at  all  frequencies.  It  is  obviously  not  a  major  effect  for  one  or  two  coupling 
events  (see  Figure  B-l),  but  it  may  have  cumulative  effects  over  numerous  coupling 
events. 

The  coupling  matrices  are  not  the  only  sources  of  deviation.  Figure  B-3  is  evi- 
dence of  the  fact  that  deviations  can  exist  between  the  two  methods  even  in  range 
independent  propagation  (RIP),  where  the  aforementioned  PPM  narrowband  cou- 
pling approximation  is  not  a  factor.  The  source  depth  dependence  of  the  deviations 
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in  this  figure  point  to  the  source  excitation  factor,  ^(^s),  as  a  probable  cause  of  the 
most  significant  deviations,  but  the  modal  wavenumber  spectra  and  modal  attenua- 
tion spectra  used  in  the  PPM  cannot  be  ruled  out  as  additional  contributors  to  the 
deviations. 
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Appendix  C 

Incorporation  of  Background 
Internal  Wave  Spectrum  through 
Dozier  and  Tappert  Coupling 


Our  primary  modeling  concern  in  this  thesis  is  internal  wave  solitons;  their  strong 
horizontal  sound  speed  gradients  cause  significant  scattering  of  energy  between  acous- 
tic modes.  Linear  internal  wave  fields,  which  are  also  present  in  shallow  water,  cause 
additional  scattering  of  energy  between  acoustic  modes.  They  are  usually  modeled 
and  studied  as  a  stochastic  process  that  produces  random  sound-speed  fluctuations, 
8c(r,  z,<),  about  a  background  profile,  cq(z).  A  "standard  method"  of  producing  these 
random  sound  speed  fluctuations  and  resultant  coupling  matrices  in  deep  water  was 
developed  by  Dozier  and  Tappert  (DT)[14].  We  will  follow  their  method  here. 


C.l      Coupled  Mode  Formulation 

The  DT  method  generally  consists  of  using  a  GM  spectrum  as  the  basis  for  creat- 
ing a  random  realization  of  the  internal  wavefield,  £(r,  z,t),  which  from  Munk  and 
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Zachariasen  [21]  is  proportional  to  Sc(r,  z,t)  through, 


N2(z) 
6c(r,z,t)a:—^ar,z,t),  (C.l) 


where  g  is  the  gravitational  acceleration  and  N(z)  is  the  Brunt-Vaisala  (buoyancy) 
frequency.  Dozier  and  Tappert  [14]  use  the  following  normal  mode  formulation  (with 
cylindrical  spreading  removed)  for  the  pressure  field, 

p(r,  z,  t)  =  e-iut  f)  &(r,  t)(kTn)-^2<f>n(z),  (C.2) 

n=l 

where  £n(ri*)  =  An(r,  t)(krny^2  is  a  modified  mode  amplitude  that  yields  a  sym- 
metric coupling  matrix.  Using  the  quasi-static  and  narrow  angle  forward  scattering 
approximations,  they  developed  the  coupled  mode  equation 

-^p(r,  t)  -  ikrn£n(r,  t)  =  -i  ^2  Rnm{r,  t)£m(r,  t),  (C.3) 

where  the  coupling  matrix  is 

ft-*'*)  =  snTT  m/2  I dz6-^j^-K{z)K(z).  (C.4) 

Creamer  [18]  closely  follows  Dozier  and  Tappert  in  his  analysis  of  shallow- water 
waveguides;  he  retains  the  GM  spectrum  assumption  and  the  basic  form  of  the  coupled 
mode  equations,  but  he  adds  a  bottom  loss  via  the  addition  of  ian  to  the  horizontal 
wave  number,  fcrn,  and  introduces  a  simplification  to  the  determination  of  8c(r,  z,t). 
Although  the  procedure  is  still  complicated,  he  bypasses  the  actual  determination  of 
the  internal  wave  field  in  favor  of  computing  a  8c(r,  z,t)  field  that  is  (1)  zero-mean 
and  gaussian  and  (2)  has  a  correlation  function  consistent  with  that  due  to  a  GM 
internal  wave  spectrum. 

It  is  unclear  as  to  how  closely  a  GM  spectrum  approximates  shallow  water  inter- 
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nal  waves.  For  example,  Lynch  et  al.  [16]  present  a  spectral  analysis  of  data  from 
a  thermistor  string  in  400  m  of  water  that  is  consistent  with  a  GM  spectrum  at  fre- 
quencies above  about  8  cpd,  but  the  spectrum  in  the  3  to  8  cpd  region  is  significantly 
less  energetic  relative  to  GM.  This  single  string  could  not  provide  any  information 
concerning  the  directional  spectrum  in  the  region. 

Absent  any  other  data  to  the  contrary,  the  GM  spectrum  assumption  provides  a 
method  of  developing  somewhat  realistic  levels  of  vertical  and  horizontal  coherence. 
The  vertical  coherence  is  critical  to  the  calculation  of  the  coupling  matrix  at  a  given 
range  and  the  horizontal  coherence  provides  some  realism  to  the  range  dependence  of 
these  matrices.  Conversely,  since  the  GM  spectrum  may  not  really  be  representative 
of  shallow  water,  the  modeling  of  correlation  functions  through  a  rather  complicated 
GM  internal  wave  synthesis  may  also  not  be  the  best  approach. 

C.l.l      EOF  Sound-speed  Perturbation  Modes 

Empirical  orthogonal  function  (EOF)  modes  of  thermistor  string  temperature  pertur- 
bations are  a  good  alternative  to  internal  wave  synthesis.  Unlike  the  GM  synthesis 
approximations,  the  EOF  modes  will  include  influence  from  both  linear  and  nonlinear 
internal  waves,  so  the  EOF  modes  can  be  used  as  stand  alone  random  realizations 
(i.e.  no  SIA  solitons  are  required). 

For  small  perturbations  in  temperature,  the  associated  sound-speed  fluctuations 
can  be  approximated  by 

6c(r,z,t)  ~  (4.6  -  AIT)  x  6T{r,z,t),  (C.5) 

where  T  is  in  degrees  centigrade.  With  proper  scaling,  the  EOF  vertical  mode  shapes 
for  the  perturbations  recorded  by  a  vertical  array  of  thermistors  can  be  used  to 
synthesize  perturbations  in  sound  speed  as  a  function  of  depth.  For  a  single  realization 
of  sound-speed  perturbations  at  a  particular  range,  using  an  average  temperature  of 
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18  degrees  centigrade  in  the  thermocline,  we  let 

TV 

6c(r,  z)  =  2.6  £  Bn(r)Um(z),  (C.6) 

where  the  Bn(r)  are  zero-mean  gaussian  random  variables,  and  the  Uvn(z)  are  the 
vertical  EOF  modes  obtained  through  a  singular  value  decomposition  of  12  days 
(over  17,000  samples)  of  data  from  an  eleven  thermistor  string  located  between  the 
source  and  receiver.  The  procedure  determines  11  vertical  modes,  but  only  the  first 
six  are  used  in  the  numerical  trials  that  follow. 

The  determination  of  horizontal  correlations  through  the  range  dependence  of 
Bn(r)  is  not  as  straight  forward.  EOF  vertical  mode  coefficients,  bn(t),  are  determined 
from  modal  analysis, 

bn(t)  =  Uvn(z)-lT(z,t),  (C.7) 

where  T(z,  t)  is  the  thermistor  string  record.  The  single  time  series  bn(z,  t)  is  then 
broken  into  23  different  time  series  that  are  one  tidal  cycle  (12.42  hrs  or  746  minutes) 
long.  The  transformation  from  time  dependence  to  range  dependence  is  made  by 
assuming  all  the  temperature  disturbances  are  propagating  without  dispersion  at  0.8 
m/s  along  the  same  line  as  the  acoustic  transmission  path.  This  is  reasonable  for 
the  nonlinear  wave  perturbations,  as  the  majority  observed  during  SWARM  traveled 
in  that  direction.  The  solibore  internal  tide  appears  from  thermistor  records  to  also 
travel  at  a  speed  of  around  .8  m/s.  Propagation  directions  offset  from  the  acoustic 
path  would  weight  the  average  phase  speed  higher.  For  linear  internal  waves,  at 
frequencies  above  about  3  cpd,  phase  speeds  of  around  .35  m/s  are  appropriate.  An 
isotropic  directional  spectrum  for  these  linear  internal  waves  yields  a  median  phase 
speed  of  .5  m/s  relative  to  the  propagation  path.  The  mean  phase  speed  is  infinite, 
but  it  lowers  to  .8  m/s,  if  the  minimum  phase  angle  considered  is  set  at  .05  instead  of 
0  radians.  The  0  radian  limit  is  not  reasonable  if  wavefront  curvature  is  considered. 
All  things  considered,  0.8  m/s  is  reasonable  estimate  for  an  "all-wave  average  phase 
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speed"  relative  to  the  propagation  path.  The  resulting  range  dependent  temperature 
profiles  never  existed  in  reality,  but  they  are  approximations  of  range  dependent 
profiles  that  could  reasonably  exist.  By  this  reasoning,  each  746  minute  set  of  mode 
coefficient  data  is  assumed  to  represent  35.8  km  in  range  through 

Bn(r)  =  bn(t  =  zrf-r-).  (C.8) 

The  23  sets  of  vertical  mode  coefficients  can  then  be  analyzed  for  horizontal  EOF 
modes,  Uhm{r)-  A  single  realization  of  Bn(r)  is  computed  by 

M 

Bn{r)  =  £  PmUhm(r),  (C.9) 

771=1 

where  the  flm(r)  are  zero-mean  gaussian  random  variables.  This  makes  Bn(r)  and 
6c(r,  z)  in  turn  zero-mean  gaussian  random  variables.  The  23  sets  of  data  produce 
22  horizontal  EOF  modes,  but  only  the  first  15  are  used  to  synthesize  sound-speed 
perturbation  fields  in  the  numerical  trials  that  follow. 

The  standard  deviations,  <7$c,  of  6c(r,  z)  are  controlled  by  the  standard  deviations, 
<y-w.ni  of  the  {3m  for  each  f//im(r)  of  each  Uvn(z).  These  standard  deviations  are  simply 
that  statistic  as  obtained  from  the  23  EOF  analysis  outputs, 

Pm  =  Uhm(r)-1Bn(r),  (CIO) 

of  the  23  sets  of  horizontal  data.  The  amn  thus  produced  yield  a  model  with  a  range 
and  depth  averaged  <  (Sc/cq)2  >  of  4.4  x  10~6.  Creamer  [18]  considers  from  10"6  to 
10-5  to  be  typical  values,  and  Tielbuerger  et  al.  [9]  report  a  value  of  5  x  10~7  for 
their  model. 
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C.1.2     Numerical  Trials 

Numerical  trials  with  the  Dozier  and  Tappert  [14]  coupled  mode  equations  using 
EOF  mode  sound-speed  perturbation  fields  will  demonstrate  their  effects  on  modal 
propagation.  The  coupled  mode  equations  (4.14)  in  matrix  form  are 

f£  =  i{KN  -  R(r)]t,  (C.ll) 

or 

where  Kn  =  diag[fci,...,fcjv]  [15].  This  equation  has  a  local  solution  that  can  be 
represented  by 

£(r  +  Ar)  =  £(r)e^-*(r)]Ar .  (C.12) 

The  implementation  of  the  above  equation  in  discrete  steps  is  obtained  through  di- 
agonalization  of  the  coupling  matrix  [Ajv  —  R(r)]Ar.  R(r),  which  is  calculated  from 
(C.4),  derives  its  range  dependence  from  a  realization  of  8c(r,z)  computed  by  (C.6). 
The  step  size,  Ar,  used  in  these  trials  is  480  m.  This  undersamples  the  fields  actual 
perturbations  in  soundspeed,  as  did  Dozier  and  Tappert's  [15]  1  km  step  size,  but  the 
idea  is  that  the  application  of  a  8c  =  ^(8c(r)  +  8c(r  +  Ar)  to  a  region  Ar  is  equivalent 
to  using  a  <  8c  >,  spatially  averaged  over  Ar,  when  the  results  are  averaged  over 
numerous  trials. 

Propagation  without  Attenuation 

When  there  is  no  attenuation,  Dozier  and  Tappert  [14]  have  shown  that  these  cou- 
pled mode  equations  conserve  energy  and  tend  towards  an  equal  partition  of  energy 
among  the  modes  involved.  The  assumption  of  no  attenuation  is  not  reasonable  in 
our  case,  but  running  the  model  without  attenuation  provides  data  for  comparison 
with  Dozier  and  Tappert  and  demonstrates  the  coupling  effectiveness  of  the  EOF 
field.  Figure  C-l  shows  that  all  the  energy  is  initially  transmitted  in  mode  1,  and  an 
equal  partition  of  energy  is  achieved  at  about  120  km.  The  rate  at  which  the  system 
tends  towards  equilibrium  is  controlled  by  the  strength  of  the  coupling.  Figure  4-11 
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Mode  1-7  Propagation  (No  Attenuation) 
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Figure  C-l:  Effects  of  linear  internal  wave  scattering  on  mode  propagation.  The 
results  shown  are  the  average  of  50  trials  using  a  sound-speed  perturbation  model 
with  a  range  and  depth  averaged  <  (6c/ cq)2  >  of  1.4  x  10~5.  The  mode  shapes  used 
in  the  coupling  equations  were  from  the  standard  field  measured  SSP  at  400  Hz. 
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represents  the  average  of  50  trials  using  random  realizations  of  sound-speed  pertur- 
bation fields  produced  by  the  EOF  sound-speed  model  described  above.  The  range 
and  depth  averaged  <  (6c/ Co)2  >  was  scaled  up  to  1.4  x  10-5,  speeding  the  approach 
to  equilibrium.  The  field  strength  is  easily  scaled  up  or  down.  When  scaled  down  it 
can  be  used  in  conjunction  with  SIA  solitons,  similar  to  a  GM  random  realization,  as 
a  background  field. 

Propagation  with  Attenuation 

When  modal  attenuation  is  added  to  the  above  coupling  equations,  Creamer  [18] 
shows  both  analytically  and  numerically  that  in  the  limit  of  a  large  range,  both 
the  average  modal  intensity,  <  I  >,  and  fourth  moment  of  pressure,  <  I2  >,  decay 
exponentially  with  range.  The  exponential  decay  rate  of  <  I2  >  is  always  less  than  or 
equal  to  the  exponential  decay  rate  of  <  7  >2,  because  their  difference,  the  variance  of 
the  intensity,  must  be  greater  than  or  equal  to  zero.  Thus  in  the  limit,  the  scintillation 

index,  SI,  defined  as 

<  I2  > 
57  =  ^-1,  (C.13) 

grows  exponentially  in  range.  Repeating  our  numerical  trials  with  modal  attenuation 
applied  reveals  an  asymptotic  SI  +  1  exponential  growth  factor,  Xi  °f  about  0.01 
km-1  (Figure  C-2).  Creamer's  [18]  isovelocity  model  with  average  GM  <  (8c/ cq)2  > 
perturbations  of  10-5  yielded  an  SI  +  1  exponential  growth  rate  of  0.02  km-1.  The 
scattering  ability  of  the  EOF  field  clearly  mimics  that  seen  with  GM  fields,  both  with 
and  without  attenuation. 

C.2     Limitations  of  Dozier  and  Tappert  Coupling 

Regardless  of  the  form  of  internal  wave  field,  the  main  limitation  of  any  model  that 
includes  Dozier  and  Tappert  [14]  coupling  is  that  the  results  are  only  valid  in  a 
stochastic  sense.  Both  Dozier  and  Tappert  [15]  and  Creamer  [18]  restrict  their  results 
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Figure  C-2:  Asymptotic  behavior  of  SI+ 1  for  mode  1.  The  plotted  quantity  decreases 
to  the  asymptotic  value  of  the  growth  rate  factor  x-  The  mode  shapes  and  attenuation 
parameters  used  in  the  coupling  equations  were  from  the  standard  field  measured  SSP 
at  400  Hz. 
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to  averages  over  a  hundred  or  more  random  realizations. 


C.2.1      Fine  Scale  Application 

The  large  step  size  employed  with  Dozier  and  Tappert  [14]  coupling  is  the  primary 
impediment  to  achieving  a  deterministic  mode  arrival  model.  If  the  implementation 
step  size  is  reduced  significantly,  one  conceivably  moves  from  predicting  the  probabil- 
ity of  energy  transfer  between  modes  over  a  large  step  to  predicting  actual  levels  of 
coupling  over  small  steps,  where  6c(r)  at  a  point  in  the  step  is  representative  of  8c(r) 
over  the  entire  step.  The  feasibility  of  this  sort  of  fine  scale  application  is  evaluated 
using  a  test  waveguide. 

Test  Case 

In  lieu  of  an  arbitrary  random  realization  of  a  GM  or  EOF  internal  wave  field,  we 
will  attempt  to  apply  these  coupled-mode  equations  to  waveguides  derived  directly 
from  range  propagated  WVLA  thermistor  records.  Using  the  arguments  previously 
discussed  for  horizontal  EOF  analysis,  the  transformation  from  time  dependence  to 
range  dependence  is  made  by  assuming  all  the  temperature  disturbances  are  propa- 
gating without  dispersion  at  0.8  m/s  along  the  same  line  as  the  acoustic  transmission 
path. 

Following  the  propagation  step,  the  temperature  records  are  interpolated  in  depth 
and  converted  to  sound  speed.  The  resulting  range  dependent  sound  speed  profiles 

^'^-^SS^-  (ai4) 

never  existed  in  reality,  but,  again,  they  are  approximations  of  range  dependent  pro- 
files that  might  reasonably  exist. 

The  next  step  involves  subtracting  a  mean  profile  to  create  a  sound  speed  pertur- 
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Figure  C-3:  Mean  sound  speed  profiles  obtained  from  the  WVLA  thermistor  records. 
The  solid  line  is  an  average  over  yearday  211-217,  and  the  dashed  line  is  averaged 
over  yearday  211  only.  The  dotted  SSP  is  taken  from  early  in  yearday  211  when  no 
solitons  were  at  the  WVLA.  All  three  are  extrapolated  for  depths  above  12.5  m  and 
below  60.5  m. 


bation  field, 


8c(r,  z,  t)  =  c(r,  z,  t)  -  Cq(z). 


(C.15) 


Herein  lies  the  difficulty,  namely  the  choice  of  the  mean  profile  to  be  used  in  computing 
the  sound  speed  perturbation  fields.  Three  possible  choices,  a  one  day  average,  a  six 
day  average,  and  an  SSP  chosen  from  a  period  without  any  solitons  are  shown  in 
Figure  C-3. 

The  choice  of  a  mean  profile  determines  the  modeshapes,  which,  as  shown  in  Figure 
C-4,  affect  the  elements  of  the  coupling  matrices.  Each  mean  profile  also  has  slightly 
different  dispersion  curves;  the  phase  and  group  speeds  control  the  phasing  between 
modes  as  a  function  of  range,  and  the  cumulative  effect  of  the  differences  in  dispersion 
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Figure  C-4:  Magnitudes  of  discrete  coupling  matrix  coefficients  using  a  48  m  step 
size  with  three  different  background  profiles.  The  solid  lines  are  computed  using  the 
yearday  211-217  mean  profile,  the  dashed  lines  are  from  the  yearday  211  mean  profile, 
and  the  dotted  lines  are  computed  using  the  yearday  211  "no  soliton"  SSP. 
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Figure  C-5:  Predictions  of  mode  1  arrivals  made  with  three  different  background 
profiles.  The  solid  lines  are  computed  using  the  yearday  211-217  mean  profile,  the 
dashed  lines  are  from  the  yearday  211  mean  profile,  and  the  dotted  lines  are  computed 
using  the  yearday  211  "no  soliton"  SSP. 


curves  and  coupling  matrices  can  yield  dissimilar  outputs  over  long  transmissions. 

The  outputs  shown  in  Figure  C-5  are  from  32  km  simulated  transmissions  that 
employ  670  propagation  and  coupling  steps.  These  arrivals  show  dramatic  differ- 
ences between  models  for  the  minute  0  and  minute  6  predictions,  but  the  relatively 
unscattered  minute  12  arrivals  are  surprisingly  similar. 

The  statistical  significance  of  the  predicted  arrivals  using  this  method  can  be 
judged  from  Figure  C-6,  where  12  hours  of  arrivals  are  compared  using  two  different 
mean  profiles.  Both  show  transitions  from  low  to  high  scattering  levels  at  about  the 
same  point  in  time,  but  the  one  day  average  mean,  as  might  be  expected,  predicts  less 
scattering  in  general.  A  near  ideal  (but  impractical)  choice  for  the  background  profile 
would  be  to  use  the  range  average  of  the  soundspeed  profiles  used  in  each  realization. 
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Figure  C-6:  Twelve  hours  of  predicted  mode  1  arrivals  using  the  yearday  211  mean 
profile  and  the  yearday  211-217  mean  profile. 
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This  would  amount  to  a  series  of  1 1  hr  averages  and  minimize  the  average  magnitude 
of  the  sound  speed  perturbations  for  each  transmission's  waveguide.  As  it  is,  the 
yearday  211  mean  profile  yields  a  model  with  a  <  (6c/ cq)2  >  of  6.4  x  10~6,  and  the 
yearday  211-217  mean  profile  yields  a  model  with  a  <  (6c/ cq)2  >  of  1.0  x  10~5. 

An  argument  can  be  made  that  in  each  case  the  soundspeed  perturbations  ob- 
tained from  WVLA  thermistor  records  are  too  large  for  the  Dozier  and  Tappert 
coupled  mode  equations  to  be  applicable.  Clearly,  as  both  long  and  short  term  fluc- 
tuations in  sound  speed  decrease,  the  choice  of  the  background  profile  becomes  more 
deterministic,  and  the  "mean  profile"  uncertainties  on  predicted  mode  1  arrivals  are 
diminished.  Dozier  and  Tappert's  [14]  deep  water  model  has  sound  speed  pertur- 
bations that  are  at  least  an  order  of  magnitude  less.  They  have  a  maximum  range 
averaged  <  (6c/ Cq)2  >  of  about  2.5  x  10~7,  and  the  perturbations  decay  exponentially 
in  depth.  Their  model  has  no  uncertainty  in  the  mean  profile,  as  it  is  assumed  from 
the  outset.  Similarly,  Creamer  [18],  who  uses  the  Dozier  and  Tappert  coupling  equa- 
tions in  his  shallow  water  model,  has  no  uncertainty  associated  with  his  isovelocity 
mean  profile,  but  his  depth  and  range  averaged  <  (6c/cq)2  >  values  10-6  and  10~5 
are  closer  to  the  perturbations  employed  in  our  model. 

Given  the  size  of  our  soundspeed  perturbations  and  the  uncertainties  associated 
with  the  background  profile,  we  cannot  expect  to  accurately  predict  a  particular 
mode  1  arrival  using  equation  (C.4).  Therefore,  regardless  of  the  step  size  used,  the 
predicted  mode  1  arrival  structures  generated  with  Dozier  and  Tappert  [14]  coupling 
(Figure  C-5  and  C-6)  are  only  significant  in  a  statistical  sense. 


187 


188 


Appendix  D 


Modal  Reciprocity 


The  simulations  discussed  in  Chapter  5  illustrate  the  sensitivity  of  mode  spread  and 
bias  statistics  to  the  position  of  solitons.  Figures  5-9,  5-10,  and  5-11,  in  particular, 
show  how  the  movement  of  two  soliton  trains  (away  from  the  receiver  and  towards 
the  source)  reduces  the  potential  spread  in  mode  arrival  envelopes. 

This  result  predicts  significant  differences  between  the  mode  arrivals  for  reciprocal 
transmissions  when  the  distribution  of  scattering  regions  between  the  transducers  is 
highly  asymmetric.  In  a  test  of  this  concept,  a  highly  asymmetric  SIA  soliton  waveg- 
uide (from  Figure  5-9),  is  used  for  simulated  reciprocal  transmissions.  The  results 
of  this  test  are  shown  in  Figure  D-l.  The  mode  arrivals  clearly  show  the  expected 
decreases  in  modal  spread  for  the  reciprocal  transmission,  so  there  is  no  reciprocity 
in  the  individual  mode  arrivals.  However,  when  the  mode  arrival  predictions  are 
synthesized  to  yield  time  domain  pulses  at  the  individual  transducers, 

p(32A;m,27m,2)  =  ^pn(32fcm,*)<£n(27m)  (D.l) 

n 

(where  thepn(r,  i)  are  the  mode  arrivals  and  the  <f)(z)  are  the  mode  shapes),  reciprocity 
between  the  source  and  receiver  is  revealed. 

A  somewhat  simpler  scenario  is  illustrated  in  Figure  D-2.    In  this  case  a  front, 
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Figure  D-l:  A  comparison  of  reciprocal  transmissions  through  an  asymmetric  SIA 
soliton  waveguide  (upper  frame).  Mode  1  through  4  arrivals  as  well  the  time  domain 
pulse  at  a  single  transducer  (27  m)  are  shown  for  both  the  forward  and  reciprocal 
transmissions.  The  source/receiver  depth  is  27  m  for  both  transducers. 
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Figure  D-2:  A  comparison  of  reciprocal  transmissions  through  a  strong  front  (upper 
frame).  Mode  1  through  4  arrivals  as  well  the  time  domain  pulse  at  the  receiver  are 
shown  for  both  the  forward  and  reciprocal  transmissions.  The  source/receiver  depth 
is  50  m  at  0  km  and  27  m  at  32  km. 
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requiring  a  single  coupling  event,  is  simulated  at  the  28  km  point  of  32  km  waveguide. 
The  standard  "non-soliton"  SSP  is  used  from  0  to  28  km  and  the  standard  "soliton" 
SSP  is  used  from  28  to  32  km.  Again  there  are  significant  differences  in  the  individual 
mode  arrivals,  but  the  sums  at  the  transducers  show  the  reciprocity  between  the 
source  and  receiver. 
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